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In order to decide whether the multiplicative theory of 
showers offers an adequate explanation of the experimental 
data on large showers, or whether it is necessary to suppose 
that they are produced by some other process, such as has 
been outlined by Heisenberg, the following calculations are 
undertaken. The energy spectrum of the electrons and 
photons incident upon a small thickness of material is 
calculated from the observed frequency distribution in size 
of bursts of ionization produced by large showers, the effect 
of fluctuations being taken into account. The energy 
distribution obtained is of the form E~*, where a is 2.6 for 
energies of the order of 10° volts, and decreases slowly with 
increasing energy, in agreement with the energy distribu- 
tion calculated by Heitler to explain the variation of 
cosmic-ray intensity with altitude. The maximum number 
of electrons necessary to give the observed frequency of 
bursts is less than 0.5 percent of the total number of 


INTRODUCTION 


HE theory of the growth and decay of the 

secondary radiations which accompany a 
high energy electron in its passage through 
matter has been worked out in some detail by 
Carlson and Oppenheimer! and by Bhabha and 
Heitler.2 This theory has been applied, with 
considerable success, to explain the variation of 
cosmic radiation with latitude and altitude by 


*A preliminary report of this paper was presented at 
the New York meeting of the American Physical Society, 
February, 1938. 

1J. F. Carlson and J. R. Oppenheimer, Phys. Rev. 51, 
220 (1937). 

?H. J. Bhabha and W. Heitler, Proc. Roy. Soc. A159, 
432 (1937). 


particles observed in a cloud chamber at sea level in the 
energy range between 10° and 10!° volts, which is in good 
accord with the cloud chamber observations. This calculated 
incident energy spectrum is utilized to calculate the number 
and frequency distribution of large bursts for large 
thicknesses of material. These calculated values differ 
considerably from the experimental ones, but this difference 
is probably to be ascribed to the effect of the penetrating 
cosmic rays, and is not to be regarded as evidence of a 
breakdown of the cascade theory. The experiments on the 
absorption of a shower are also shown to be in harmony 
with the theoretical estimates. It is concluded that no 
mechanism involving the production of many shower 
particles in a single act need be invoked to explain the 
occurrence of large showers, but that the ordinary multipli- 
cative processes are entirely adequate when proper account 
is taken of the fluctuations. 


Heitler® and Nordheim.‘ The properties of small 
showers of rays are also well accounted for in this 
way. It remains to show whether the large 
showers, such as produce bursts of ionization, or 
Hoffmann Stésse, are consistent with this 
mechanism of production or not. It has been 
realized for some time that the formation of a 
large shower is not the result of a single act, but 
that some sort of a multiplicative process must be 
involved.> To what extent this ‘‘production of 
~ 3 W. Heitler, Proc. Roy. Soc. A161, 261 (1937). 

*L. W. Nordheim, Phys. Rev. 51, 1110 (1937). 

°C. G. Montgomery, Phys. Rev. 45, 62 (1934); G. L. 
Locher, Phys. Rev. 44, 779 (1933); C. G. way my 


D. D. Montgomery and W. F. G. Swann, Phys. Rev. 4 
512 (1935). 
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956 C. G. MONTGOMERY 
one shower by another”’ is important has not 
been clear. On the other hand, Heisenberg® has 
suggested a mechanism whereby a large number 
of electrons would be produced in the interaction 
with a single nucleus. It is the purpose of this 
paper to inquire whether some such mechanism 
as is postulated by Heisenberg is necessary to 
explain the experimental properties of large 
showers, or whether successive pair formation by 
the radiations resulting from the slowing down 
of high energy electrons is sufficient. 


CALCULATION OF THE ENERGY DISTRIBUTION 
OF INCIDENT RAys 

Suppose we have a piece of heavy material, say 
lead, placed above an ionization chamber, and we 
make observations of the rates of occurrence of 
the bursts of ionization as a function of the size 
of the burst. We suppose that the ionization 
produced is proportional to the number of 
electrons which make up the shower. Then it has 
been pointed out’ that the frequency of showers, 
R(N)dN, containing numbers of rays between NV 
and N+dN can be well represented by the 
empirical expression 

R(N) =A/N’. 


This expression is valid over quite a large range 
of N, and its form is almost independent of the 
thickness of the shower producing material. The 
value of the exponent necessary to obtain a best 
fit of the data varies somewhat with the obser- 
vations considered, but we may take it as 
integral for simplicity. 

The cascade theory of showers is given in the 
most convenient form for calculation by Bhabha 
and Heitler. They give, in graphical form, the 
expected number, ¢, of electrons, positive and 
negative, of energy greater than a critical value, 
E., (E.=10' volts for lead) which emerge from a 
plate of a given thickness when an incident 
electron of energy £ falls upon the plate. Let us 
first consider the case of a lead plate 1.2 cm 
thick. The actual number of rays, N, which 
emerge from the plate will not, of course, in 
general be equal to this expected number, e. 
Furry* has shown that, for thin pieces of material, 
"6 W. Heisenberg, Zeits. f. Physik 101, 533 (1936). 

7C. G. Montgomery and D. D. Montgomery, Phys. Rev. 


48, 786, 969 (1935). 
8 W. H. Furry, Phys. Rev. 52, 569 (1937). 
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Fic 1. The energy distribution of electrons and photons 
necessary to produce the observed distribution in size of 
large showers. 


the probability of finding N rays when « rays are 
expected is given by the expression (1 —1/e)*~'/e. 
Let us suppose that the number of ion pairs 
observed in our chamber is proportional to the 
number of rays above the critical energy. Then if 
we have a number of electrons F(£)dE of energy 
between E and E+dE falling upon the lead plate, 
the number of showers containing N rays will be 
given by 


oo 
R(N) =4/N= { e1(1—e")4"F(E)dE. 
Ec 
We thus have an integral equation which can be 
solved for F(Z). An approximate solution of this 
can easily be obtained. We note that, for a given 
thickness, € is a function of E only, so that a 
change of variable gives us 


R(N) -{ e"(1—e")*f(e)de. 
1 


If we assume f(e) is of the form B/e*, the integral 
may be evaluated directly, and we have 
R(N) =BI(s)P(N)/T(s+N). 
Now if we equate BI'(s) =A and s=3, we obtain 
R(N) =A/N(N+1)(N+2), 
which, for large N, is indistinguishable from the 
experimental A/N?*, and hence f(e) =A/2é is the 
desired solution. 
We may then immediately apply Bhabha and 


Heitler’s calculations, and find which F(£) 
corresponds to f(e)=A/2é. Fig. 1 gives the 
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values of F(E) so calculated. This energy 
distribution may be considered to be proportional 
to E-* where a is 2.6 for E=10° volts, and a 
decreases with increasing E to a value of 2.1 for 
E=10" volts. This decrease in a cannot continue 
indefinitely, of course, since the total energy of 
the distribution would diverge. It may be shown, 
however, that if the function R= A/N* is a good 
representation of the experimental data from 
N=10 to N=1000, then the function F(£) is 
determined as above from 10° to 10" volts. The 
F(E) calculated here has been spoken of as if it 
were a distribution of electrons which produced 
the observed showers. If, on the other hand, we 
had carried through the calculations assuming 
that all the observed showers were produced by 
photons, we would have been led to the same 
F(E), which, however, would then represent an 
energy distribution of photons, since the behavior 
of a high energy photon is almost identical with 
that of an electron. Actually we know that both 
electrons and photons produce the showers, and 
we must regard F(Z) as made up of both entities 
in amounts of comparable importance. Since the 
shower producing effects are independent of the 
proportion of electrons and photons, we cannot 
utilize the shower phenomena to give us any 
information regarding this proportion. 

We may compare our calculated form of F(£) 
with the energy distribution of electrons at sea 
level obtained in another way. Heitler,* in order 
to explain the variation of the soft component of 
the cosmic radiation with elevation, assumes a 
distribution of electrons entering the earth’s 
atmosphere. This results in an energy distribution 
of electrons at sea level which is of the form 
E--5 approximately.® Thus there is good agree- 
ment between the two methods of calculation, 
which are based on quite independent data. 

We have determined the form of the energy 
distribution of the shower producing rays; let us 
now determine the number of them necessary to 
produce the observed number of showers. The 
quantity measured as the size of a shower is the 
number of ion pairs produced. Since N is the 
number of rays above the critical energy, E., and 
since there are roughly an equal number of rays 
above and below the critical energy,'® we should 


® Reference 3, page 276. 
10H. J. Bhabha, Proc. Roy. Soc. A164, 257 (1938). 


choose an “‘effective specific ionization” approxi- 
mately twice the specific ionization of a single 
electron. Now there is considerable uncertainty in 
this quantity, but the maximum value of the 
constant A will be given by the minimum value 
of the specific ionization. We shall take this to be 
30 ion pairs per cm. We choose for illustration the 
rate of occurrence of bursts greater than 1.5 x 10° 
ion pairs produced in a magnesium chamber 
containing about 15 atmospheres of nitrogen and 
covered by one centimeter of lead.'' We are thus 
dealing with showers containing of the order of 
100 rays. These data result in an upper limit for 
A of 9.8 X10 cm~ sec... If we now consider the 
region of the energy spectrum from 10° to 10!° 
volts, we should need only 4.0X10-° electrons 
and photons per sq. cm per second to account for 
all the showers produced. The total number of 
charged particles observed in a cloud chamber in 
this energy range may be taken from the 
observations of Blackett," and is 950 X 10-* cm~? 
sec. Thus at most only 0.5 percent of the 
particles in this energy range observed in a cloud 
chamber need be electrons. That there are some 
electrons is shown by the fact that showers have 
been observed, in cloud chamber photographs, 
whose total energies are greater than 10° volts". 
Thus we have good consistency with cloud 
chamber observations. Similarly even the largest 
bursts which have been reported can be produced 
electrons or photons whose number is a negligibly 
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Fic. 2. The contribution to the number of 100 ray 
showers by the rays of various energies. The vertical line, 
Em, marks the median energy. 


1C. G. Montgomery and D. D. Montgomery, Phys. 
Rev. 48, 786 (1935). 

2 P. M.S. Blackett, Proc. Roy. Soc. A159, 1 (1937). 

SC. D. Anderson and S. H. Neddermeyer, Phys. Rev. 
50, 263 (1936). 














958 C. G. MONTGOMERY 
small proportion of the total number of particles. 
We see that it is unnecessary to invoke any 
special mechanism for the production of large 
showers such as suggested by Heisenberg. 

It is interesting to see what energies are 
effective in producing showers of a given size. 
Fig. 2 gives the contribution of each energy to the 
total number of showers for V=100. The lower 
scale of abscissae gives the expected number of 
rays, e, for each energy. We note the striking 
fact that the rays which are most effective in 
producing 100 ray would 
expected to produce about 20 rays, and that 
the rays for which we expect 100 shower particles 
give a wholly negligible contribution. Table I 
shows some characteristics of this distribution. 
Thus all the observed showers of 100 rays are 
fluctuations from expectation. 

It might seem from a cursory examination of 
the problem that since f(e) and R(N) are the 
same functions of their respective variables, the 
effect of the fluctuations is only to introduce a 
small constant factor. However, this is quite 
illusory since, in the absence of the fluctuation 
phenomena, the rays which would be effective in 
producing 100 ray showers would not lie in the 
energy range between 10° and 10!° volts, but 
would have considerably higher energies in the 
neighborhood of 5X10! volts. Thus in the 
absence of fluctuations we should need, in order 
to account for the observed number of bursts, 
considerably more electrons of energy around 10° 
volts than all the particles of this energy ob- 
served. It is the neglect of the fluctuations which 
has led many investigators," ourselves included,!® 
to the view that it is necessary to have some 
mechanism whereby many shower particles are 
produced in a single act to account for the large 
showers from small thicknesses. However, when 
the fluctuations are properly taken into account, 
the ordinary cascade picture is entirely adequate. 


showers only be 


TABLE I. Characteristics of distribution in Fig. 2. 











ENERGY 
Most probable 5.0 X 10° volts 20 
Median 1.2 101° 35 
Mean 4.3 X 101° 51 








4 Reference 10, page 265; H. Euler, Physik. Zeits. 38, 
943 (1937). 

1’ C. G. Montgomery and D. D. Montgomery, Phys. 
Rev. 50, 490 (1936). 








AND D. D. MONTGOMERY 
& 

6- 

LOG,RCN . 
447 4 
2- 
° Lj T 
10 100 1000 «ON 


Fic. 3. The variation of frequency distribution in size of 
showers with thickness of the producing material. The 
numbers indicate the thickness of lead in centimeters. 


THE VARIATION OF THE PRODUCTION OF LARGE 
SHOWERS WITH THICKNESS OF 
THE MATERIAL 


We may use the energy spectrum of the 
incident rays derived above to calculate how 
many showers would be expected at large 
thicknesses of material. However, we enter here 
on more uncertain ground. Although the effect of 
the fluctuations is all important at small thick- 
nesses, its importance should become less and less 
as the thickness increases. In the absence of a 
more accurate knowledge of the effect of fluctu- 
ations, the best approximation possible at the 
present time would appear to be to neglect the 
fluctuations entirely for large thicknesses. The 
calculation of the number of showers is then quite 
straightforward, and we give only the results. 
Fig. 3 shows the distributions in shower sizes to 
be expected at several thicknesses of lead. The 
lines in the figure have been drawn as straight, 
although they have a slight curvature concave 
upward. 

These results differ considerably from experi- 
ment.'® This disagreement may be expressed in 
this way: for showers of a given size, the pre- 
dicted decrease at large thicknesses is much too 
large, and the shapes of the distribution in size 
curves, that is, the slopes of the lines in Fig. 3, 
change too much with thickness. The experi- 
mental data do show some decrease in slope, but 
not nearly as large a one as is indicated here. 
This is usually expressed as a shift of the maxi- 


16 E.g., R. T. Young, Phys. Rev. 52, 559 (1937); J. K. 
Béggild, Diss. Copenhagen, 1937. 
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Fic. 4. Observations on the absorption of showers. Open 
circles: data of H. Nie recalculated by the authors; double 
circle: authors’ measurements, first method; filled circles: 
authors’ measurements, second method. 


mum to larger thicknesses in the relation between 
the number of showers and the thickness of the 
material, as the size of the showers under 
consideration increases. These points of dis- 
agreement are probably the result of neglecting 
the effects of the penetrating cosmic radiation, 
and probably do not represent any fault in the 
theory of showers produced by electrons or 
photons. As this point of view has been much 
discussed previously,'®:'7 we merely wish to 
emphasize that there need be no discrepancy in 
the cascade theory of showers here. 


THE ABSORPTION OF A SHOWER 


Another means of testing the multiplicative 
theory of showers is to utilize the experiments on 
the passage of a shower through lead. Such 


17 R. H. Woodward, Phys. Rev. 49, 711 (1936); H. Euler, 
PR R. T. Young and J. C. Street, Phys. Rev. 52, 552 
(1937). 
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experiments have been performed by Nie!* and 
the authors,'® and lead to results substantially in 
agreement. Fig. 4 shows the results of these 
experiments. Again the showers dealt with con- 
tain approximately 100 rays. Now, on account of 
the importance of the fluctuations, we cannot 
expect to be able to treat this question exactly, 
but can only obtain a rough estimate of the 
absorption from the theory. We may represent the 
experimental points by the equation N= Noe, 
where JN is the number of rays in the shower after 
it has passed through ¢ cm of lead, No the 
original number of rays (No = 100), and u«=} cm™ 
of lead. Nordheim‘ has shown that a particular 
solution of the diffusion equations of Carlson and 
Oppenheimer is given by 


P(t, E)=eE, 


for energies above the critical energy, £,, with 
the relation between m and uz: 


0.4u=4/3—1/n—[(3 —1/n)?+4/3n(n—1) }} 


when yu is measured in cm~ of lead. Now if we 
suppose that the No rays in the ‘average’ 
shower have an energy distribution of the form 
E-", we can compute what the value of m must be 
in order that the shower will be absorbed 
exponentially with the observed coefficient. This 
leads to the value n = 2.07. The total energy of the 
shower iseasily seen to be E= NoE.(n—1)/(n—2), 
or for this case, E=1.5 X10! volts. This ‘‘aver- 
age’ energy of the shower is equal to the energy 
of the incident ray, and may be compared with 
the energies listed in Table I for 100 ray showers. 
It is seen to be of the correct order of magnitude. 
Thus, although these calculations are only 
approximate, here again we obtain agreement 
with the cascade theory, and find no evidence 
that any other mechanism is necessary to 
account for the observed behavior of large 
showers. 

The authors wish to thank Professor W. F. G. 
Swann for much valuable discussion of these 
matters. 

18H. Nie, Zeits. f. Physik 99, 776 (1936). 


19C, G. Montgomery and D. D. Montgomery, Phys. 
Rev. 49, 705 (1936); Zeits. f. Physik 102, 534 (1936). 
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Transition Effects of Cosmic Rays in the Atmosphere 


H. SNYDER 
University of California, Berkeley, California 
(Received April 22, 1938) 


In this paper we give an improved treatment of the 
multiplication and absorption of cosmic-ray electrons and 
gamma-rays in the atmosphere. After a brief recapitula- 
tion, in I, of the principal processes involved in the cascade 
theory of showers and the resulting diffusion equations, we 
give, in II and III, a solution of these diffusion equations, 
making no further approximations. In IV, an expression 
for the probable number of charged particles is given; the 
values of the tunctions necessary for numerical calcula- 


tions are given in a tabular form. In V, it is shown how 
the Bhabha-Heitler cut-off energy method for estimating 
the effect of ionization losses may be consistently used. 
In VI, a comparison is made with the Carlson-Oppenheimer 
calculations. The differences near the top of the atmosphere 
are small and in a direction to improve the agreement with 
observations. At sea level the new calculations give an 
intensity roughly % of the Carlson-Oppenheimer results. 





ECENT experiments with ionization cham- 

bers and vertical coincidence counters have 
extended the study of cosmic-ray intensities 
nearly to the top of the atmosphere. In this paper 
we shall return to the question of a theoretical 
estimate of the intensity to be expected on the 
assumption that cosmic rays consist only of 
electrons, positrons, and gamma-rays. It is known 
that this assumption fails to give an adequate 
description of the observations, and that other 
phenomena, which have to do with the produc- 
tion and absorption of the ‘‘penetrating’’ com- 
ponent, must be essentially involved. We want 
here to improve the theoretical treatment of the 
multiplication and absorption of the “‘soft’’ 
radiation, in the hope that this will help, when 
further experimental results are available, to 
throw light on the behavior of the penetrating 
component. 

Carlson and Oppenheimer! have given an 
expression for the average number of charged 
particles that would be expected at a given depth 
below the top of the atmosphere, provided only 
that this distance is not too small. In their 
calculations they made two types of approxi- 
mations: (a) Analytic forms for the probabilities 
of pair production and gamma-radiation which 
closely approximate the quantum theoretical 
expressions in the high energy region were used; 
(b) Whereas, in calculating the energy distri- 
bution of emitted gamma-rays, the theoretical 
spectrum with constant intensities was used, a 


! Carlson and Oppenheimer, Phys. Rev. 51, 220 (1937). 
Referred to as C.O. 


distribution corresponding to constant proba- 
bility for various quanta was used to estimate the 
direct effect of radiation on the energy distri- 
bution of the charged particles. We have found it 
possible to calculate the probable number of 
charged particles at a given depth without 
making the second approximation, and to extend 
the calculation to the case of small thicknesses. 
The approximations (a) are very good, and these 
have been retained. 


I 


For convenience, a brief resume of the principal 
processes will be given. The three essential 
processes which will be considered are: (1) pro- 
duction of gamma-rays by charged particles, 
(2) pair production by gamma-rays, (3) ioniza- 
tion losses of the electrons and positrons. The 
Compton recoil of electrons under impact by 
photons was not explicitly introduced, because in 
air at the low energies for which pair production 
gives too few particles, the recoil electrons 
absorb enough gamma-radiation to keep the 
gamma-ray absorption coefficient nearly con- 
stant. The Compton effect produces one particle 
per photon, whereas pair production gives two. 
Either process contributes the same quota to the 
ionization, since such charged particles have 
energies so low that they are stopped by ionization 
losses so quickly that they radiate very little. No 
consideration will be given to other processes, as 
production of high energy secondaries or pairs in 
the nuclear fields by charged particles, because 
they are relatively much rarer, and because the 
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TRANSITION EFFEC 


resultant particles would be nearly as penetrating 
as the incident one. 

The analytically amenable expressions for the 
probabilities of pair production, emission of 
gamma-radiation, and ionization losses, which 
were chosen by Carlson and Oppenheimer will be 
used for these calculations. 

The probability that an electron or positron of 
energy greater than E shall radiate a photon of 
energy between E and E+AE, while passing 
through matter a distance Ax, will be taken to be 

PAEAx =(K/E)AEAx. (1) 


If the atomic number of the matter is Z and 
the density of atoms is N, 
4Z*e®N 


K =——— In (200/Z!). 
im*c® 





The probability that a gamma-ray of energy E 
shall produce a pair having energies between E’ 
and E’+AE’, and E—E’ and E—E’—-AE’, will 
be taken to be 

P’ AE’ Ax = (K’'/E)AE’ Ax’ (2) 
with E’ <E, and for which 
K'/K=o~}. 

The ionization losses will be assumed to be 
independent of the energy and will be evaluated 
at energies for which they are important. If one 
regauges the measure of distance by the trans- 
formation t= Kx, the energy loss by ionization 
per unit distance becomes 


= (4rNZe'/Kmc?) In (8/ZRh). (3) 


Let N(t, E) be the probable number of charged 
particles of energy greater than E to be found at a 
thickness ¢. Further, take y(t, E)AE to be the 
probable number of photons having energies 
between E and E+AE at t. Then, in accordance 
with (1), (2) and (3), diffusion equations take the 


form 
dy/dt= —ovy+ N/E, (4) 


aN aN y(E”) 
antl ef 7 aR" 
at 








E N(E+E’) an 
+ { = dE’ 
“9 E 
© N(E+E’)—N(E’) 
+f - -dE’. (5) 
E E’ 
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If one supposes that the incoming radiation is a 
charged particle of energy Eo, the boundary 
conditions which the solution to (4) and (5) must 
satisfy are 

1 E<Ep 
y(t=0, E)=0; N(t=0, E)= (6) 
0 E>Es. 
Il 


An exact solution of Eqs. (4) and (5) with the 
boundary conditions (6) may be found for 8=0. 
If one writes 


vyo=e~*'s0/E, (7) 
then it follows that 
No =e~*'2o. (8) 


In Eq. (8) and succeeding equations a dot will 
represent differentiation with respect to ¢. The 
substitution of (7) and (8) in (5) gives rise to the 
integral -differential equation 


“0 © 2o(E’’) 
K (tu) = 20—oiy—20 | ar’ [ : dk” 
JE Je (BE)? 


E 3(E+E’) — i0(E) 
-{ asics lisa 
0 EF’ 


and 2o( E+E’) —2(E’) 
-{ — —— dk’'=0. (9) 
E | 


Simple solutions of the form z9=e*'E~” exist, 
provided yu is one of the roots of the quadratic 
equation 


uw? —(o+L(y))u—20e/y(y+1) =0, (10) 
in which 
1 (1+x)-"-1 (1b x)- ¥— x 
L(y) = { ——_———dx + = ‘dx 
0 x ae | x 
=—[¥(y)+7] (11) 
¥(y) =(d/dy) InT(y+1); y=0.577---. 


The roots of (10) are 


Me _ oth (y) o+Liy)\? 20 ' 
(2 
v 2 2 y(y+1) 


u>v. (12) 
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The general solution of (9) is now given by 


1 dys,E\~» 
n=—| - (-) [A(v)e**+B(y)e"] (13) 
2ridc y \ Ey 


if the contour C and A(y) and B(y) are de- 
termined so that this integral exists. 
To satisfy the boundary conditions (6) one 
must have 
1 E<Eo 


z(¢t=0, E)=0; 2(t=0, E)= 
0 E>Ebo. 
The contour for the integration with respect to 
y was chosen to be a straight line parallel to the 
imaginary axis, at a distance 6 to the right, and 
from —i~ +6 toio +6. 
Equation (13) may then be written 


1 dy 
—e[ A (y)e“"+B(y)e’*] 


2rid ¢ ¥y 


20>-— 
with A=In (E/E). 


From 2(t=0, £)=0, it must be true that 


B(y) = —A(y). (14) 
The second boundary condition leads to 
1 dy 1 A>0 
— —eA(y)[u—v]= (15) 
2ridc y 0 A<0. 


Multiplying (15) by e~%* and integrating with 
respect to \ from 0 to ©, one gets 
1 dy uv 


— | —A(y)—-=-—— for 
2mid¢ y y—Yo Yo 


R(yo) >0. (16) 


The integral in (16) may be evaluated by 
taking the residue at y=yo, from which one 
obtains 


A(y)=1/(u—»). (17) 


The solution satisfying the boundary condi- 
tions (6) is: 


1 pdy_ evt—er 


Sa — ey 
2m1 cy 





(18) 


Z0 
un-v 

Equation (18) is an analytic expression of the 
solution of the cosmic-ray problem considered by 


Bhabha and Heitler.? 
? Bhabha and Heitler, Proc. Roy. Soc. A159, 432 (1937). 
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III 
If one takes as before 


y=oe"Z/E; N=e~*z, (19) 


a solution of (4) and (5) with 80 which satisfies 
the boundary conditions (6) to terms of order 
8/E 9 may be found. The equation which z must 
satisfy is: 

K(z) =602/dE (20) 
with the operator K defined in Eq. (9). 

In analogy with the work of Carlson and 
Oppenheimer, we look for simple solutions of 
(20) which depend upon a parameter, say y, and 
which may be combined to satisfy the boundary 
conditions. We set 


ext E\~-»7 E\-* 
2(y) -— fas(—) (—) 
2r1 8 Eo B 


~ T(—s)l(yts) 


xX Cy, s) (21) 





r'(y) 


The contour for the integration with respect to 
sis from —i toi. The contour is kept to the 
left of the poles at 0, 1, 2, etc. and to the right of 
the poles at —y, —y—1, —y—2, etc. We look for 
a C(y, s) which for R(y) >0 is free of singularities 
in the strip 0=R(s)=1, since one can see that 
this condition is necessary to make the total 
number of particles in the shower finite. Substi- 
tuting (21) into (20), one obtains: 


E\-*s,E\~- r'(—s)T(yt+s) 
fa(—) (—) C(y, Ss) 
s \Eo B r'(y) 


x {u?—(o +L(y+s))u—20/(y+s)(yts+1)} 


E\-*s/,E\7-*" r'(—s)(y+s) 
= — fas( ) (—) C(y, s) 
Eo B r'(y) 


Xu(yts). (22) 











In order that the integrals over E obtained by 
this substitution converge, it is necessary that 
R(y+s)>0, although the expressions thus ob- 
tained may be analytically extended into the left 
half-plane. One may now note that the coefficient 
of C(y, s) in the left-hand side of (22) differs from 
that of the right side by a first power in £. If one 














(19) 


itishes 
order 
> must 


(20) 


1 and 
ms of 
y, and 
ndary 


(21) 


ect to 
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+s) 
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then takes the C(y, s) to satisfy the functional 
relation 
{u?—(o+L(yt+s))u—20/(y+s)(v+s+1)} 


s(yv+s—1)uC(y, s—1) 
xX Cy, s) = ~ (23) 
(y+s) 





(21) is at least a formal solution of (20). To 
insure that (21) be a solution of (20), y will be 
restricted to R(y)>0, and yu will be taken to 
satisfy (10). When these conditions are satisfied, 
one may deform the contour for the integration 
with respect to s for the last term in (22) and 
show that (21) is a solution of (20). 

If uw? is eliminated from (23) by means of (10), 
it may be written: 


{ul v(vt+s) —¥(y) ]+20[1/y(y+1) 
—1/(y+s)(y+s+1)]}C(y, s) = 
s(y+s—1)y 


(y+s) 


Thus for R(y) >0, C(y, s) is free of singularities 
for —1— R(y) <R(s). 

As there are two roots to (10), there are two 
linearly independent solutions of (20) for each 
value of y. The C(y, s) corresponding to the two 
roots » and » will be denoted by C,(y, s) and 
C,(y, s). In addition, the C’s will be chosen so 
that C,(y, 0) =C,(y, 0)=1. A solution satisfying 
the boundary conditions (6) to terms of order 
B/Eo may now be written: 


1 dy E\-* l'(—s)'(y+s) 
$= —-— —e f'as(—) 
4r* cy 8 B r (y ) 


| Cay, s)e*—C,(y, s)e" 








C(y,s—1). (24) 








(25) 


uv 


for evaluating the integral over s in terms of the 
residues of the poles at 0, 1, 2, etc., we find 


1 dy ert—ert 1 B 
sa— | — oe ——_ - — dye(— 
2ridc y uy 2riv ¢ E 





Cy(y, 1)e**—C,(y, 1)e” 
x| | +008 E)?+--+. (26) 
uv 


From (26) it follows that (25) satisfies the 
boundary conditions (6) for E> 8, at least to 


order 6/£; this suggests that the total energy of 
the shower is Ey to terms of order 8. This is 
shown to be true in Eq. (41). Further, the total 
number of particles of energy greater than zero at 
t=0 determined from (25) is 1+0(8/E»), as 
8/E 0, which follows from Eqs. (28) and (34). 
It is only for large Eo/8 that (25) satisfies the 
boundary conditions, and it is only for this case 
that we shall calculate the total number of 
particles. 


IV 


For the purpose of calculating the number of 
charged particles in the shower, deform the 
contour to the left. There are no singularities 
between s=0 and —y. In the limit E=0, the 
only contribution to the integral comes from the 
simple pole at s= —y. One then has 

1 dy (Cul, —yert'—Cy, — ye") 
(E=0)=— | — —— : 


2ridcy uv J 


e=In (Eo/B). (27) 


From Eqs. (19) and (27), one can easily obtain 


e~*' edy 
ya [Pe 
2ridec y 


Cy(y, —v)me**—C,(y, —y)ver"| 
x|— Diasec, iil 
u-» 


It is convenient to divide (28) into two parts, 


1 edyfuC,(y, —y) 
nat f FP) prea cy 
2midc y uy 


1 dy vC,(y, —y) 
nat fF —Meeon Go 
2ridc y um 


The first part, (29), gives the principal 
contribution to the shower. In the second, » 
is less than zero for all y on the contour, so (30) 
falls off exponentially with increasing /. 

For large values of ¢ and ¢, (29) may be calcu- 
lated by the saddle point method. The quantity 





uC, (vy, —y)/(u—v) (31) 


is a slowly varying function of y, and may be 
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taken from under the integral sign and evaluated 
at the saddle point. One then obtains 





Ca lVer — Va) peet¥et (eset 
N,=- (32) 
[2ar(us’ye2t+1) }'(us— ve) 
with the saddle point determined by 
us t+e—1/y,=0, (33) 


in which p’ =dy/dy; uw” =dy' /dy. 

‘In fact, (32) gives satisfactory results even for 
small ¢. For small values of ¢ and large values of « 
the saddle point determined by (33) is at y very 
near zero, so it is possible to approximate yu by 
(20/y)!. When this is effected, (29) may be 
evaluated in a rapidly converging series, and one 
obtains 

Cu(Ve, —Vs)Ms oo [(2oe)*t]” 
= e~°! . (34) 
n=0 '(An-+1)P(n+1) 








Nu 
(Ms agit Vs) 
Numerical calculation of (34) shows that even 
for =0 the corrections to be applied to (32) are 
small; for t=0, they amount to a decrease of 
eight percent. 

The function C,(y, —y) is a very complicated 
but slowly varying analytic function of y for all y 
near the saddle points determined from (33). It is, 
however, easily calculated from (24) for integral 
values of y. Its values for these points were 
plotted, and a smooth interpolating curve drawn 
through them. 

If the substitutions 


H(y) =(27)C,(y, —y)u/(u—¥), (35) 


kR=yu-o, a=—yp'y, b=p"y? 


are made in Eqs. (32) and (33), they may be 








TABLE I. 

y k a | b | H 
0.2 1.8876 1.485 2.06 0.286 
0.4 0.9538 1.207 1.51 0.330 
0.6 0.4944 1.063 1.41 0.351 
0.8 0.2026 0.957 1.37 0.352 
1.0 0.0000 0.857 1.35 0.342 
1.2 —0.1479 0.762 1.33 0.320 
1.4 — 0.2588 0.669 1.30 0.294 
1.6 — 0.3418 0.581 1.26 0.267 
1.8 — 0.4055 0.504 1.18 0.243 
2.0 — 0.4542 0.428 1.08 0.221 
2.2 — 0.4920 0.366 0.99 0.201 
2.4 —0.5215 0.313 0.91 0.185 
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written in the simple forms 


N,= (1+ 01) Het (36) 


and 
t=(ey—1)/a. (37) 

In Table I, the values of a, 6, k, and H/ are 
given over a range of values sufficiently large to 
cover showers from the top of the atmosphere to 
sea level. The value of ¢ given by (37) increases 
uniformly as a function of y. 

For air the value of B is about 90 Mev. ¢ is in 
units of 0.4 m water equivalent. 

N, is completely negligible for values of ¢ 
greater than 0.5. For values of ¢ less than this, it 
is easy to estimate its contribution to the size of 
the shower: (a) VN, =3, N,=3 at t=0; (b) at /=0 
the slope of NV, is equal and opposite to that of 
N,; (c) N, falls exponentially to a half-value of 
} at about ¢=0.2. 

V 

We will now look at Eq. (18) and see if we can 
find a rational way in which to introduce a low 
energy cut-off by means of which we can estimate 
the total size of the shower. This is the procedure 
suggested by Bhabha and Heitler.? We introduce 
E., supposed independent of ¢, and assume that 
the energy distribution of particles of energy 
greater than £, is not disturbed by the ionization 
losses, and that because of the ionization losses 
there are no particles of energy lower than this. 
Since the initial energy of the shower Ey» must be 
absorbed by ionization, we must have 


hed 


| No(E,.)dt= Eo/B 


0 
1 a.) ~dy pele—7 they 
-—| a { ————--- 
2mido0 oe uv 


with A,=In (Eo/E.). 
These integrals are easy to evaluate if we 
neglect terms of order B/E; then (38) gives 


(38) 


E.=|— 
[u(1)—v(1) Ju’(1) 





u(1)(—1) 


1 
=—__—_—8=0.4678. (39) 


CAi, —1) 
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(37) 
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In accordance with (39), the number of 
particles in the shower becomes 
(0.467 )—9y,e°4 thes 08 
Ny=-- (40) 
[2r(us’’vst+1) }'(us—vs) 
with the saddle point determined as before by 
Eq. (37). This equation, (40), is of the same form 
as (32) with C(y, —y) replaced by 


el ( 1 ) 
u(1) ~ (0.4677 | 


For y=0 and 1, these give identical results, and 
for y=2 they give 3.28 and 4.58, respectively. It 
should be observed that the proper ‘‘cut-off 
energy” is less than half 8. 

Calculating the total energy of the shower, as 
given by (29), and again neglecting terms of 
order 8/Eo, we get for this energy 

Cy(1, —1)u(1) 
— —————— Ey) = E. (41) 
[u(1) —v(1) Ju’(1) 
VI 

In their work, Carlson and Oppenheimer 
replace the last two terms of Eq. (5) by an 
approximate expression which enabled them to 
find a differential equation for z. The approxi- 
mation was made in such a manner that the 
number of particles lost in one energy range 
appeared in some other, and so that the average 
energy lost by radiation by the charged particles 
was gained by the gamma-rays. The approxi- 
mation distorts the energy distribution among 
the charged particles so that there are too many 
very high energy particles, too many low energy 
particles, and too few of intermediate energy 
(~8). In Fig. 1 the effect of this distortion on the 
results for the average number of particles is 
shown for Eo=1008. The curve marked J was 
obtained from the C. O. calculations, and the 
other, marked J/, from Eqs. (36) and (37). Their 
approximation : (a) slightly lowers the height of 











Fic. 1. J is a plot of the number of particles against ¢ as 
determined from the C. O. calculations. JJ is a plot of the 
same as given by Eqs. (36) and (37). JJJ shows the 
contribution to the shower by N,. The calculations were 
made for Eo=1008. One unit of ¢ is about 0.4 m water 
equivalent. 


the maximum, (b) moves the maximum about 
0.2 m to greater thicknesses, (c) considerably 
overestimates the size of the shower at great 
depths. At ten meters depth, near sea level, the 
C. O. formula gives N=7.5X10-*, and we find 
N=4.5X10-%, for a 9X10* Mev particle with 
vertical incidence. 

The changes in the theoretical curve near the 
top of the atmosphere help to remove some of the 
discrepancies found in this region by Bowen, 
Millikan and Neher* between their experimental 
results, and the calculation of C. @. applied to an 
isotropic incident distribution. So detailed a 
comparison can hardly be significant, however, 
unless it takes into account the large deviations 
from isotropy to be expected from the effects of 
the earth’s field. Near sea level, on the other 
hand, our results are about half those of C. O. 
and thus accentuate still further the discrepancies 
with the experimental curves, discrepancies 
clearly associated with the penetrating com- 
ponent. 

The author wishes to thank Professor J. R. 
Oppenheimer and Dr. R. Serber for many 
valuable suggestions that have furthered the 
writing of this paper. 
~ 3 Bowen, Millikan and Neher, Phys. Rev. 52, 80 (1937). 
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Measurements of x-ray reflection intensities at two temperatures have been made on the 
2131 and 1015 lines from an alloy of copper and zinc. The idea is put forward that the atomic 
vibrations should be greatest along the c axis of the crystal in case the axial ratio is greater 
than 1.633, least when the axial ratio is less than 1.633. The values obtained for the relative 
amplitudes of vibration along and perpendicular to the c axis confirm this idea. The relative 
atomic structure factors are also found and show a surprising dependence on the angle which 


the normal to the reflecting plane makes with the ¢ axis. 





1. INTRODUCTION 


EVERAL workers! have shown that the 
atomic vibrations in zinc crystals are aniso- 
tropic. Brindley’s? work on cadmium crystals 
indicates that its atomic vibrations are also 
anisotropic. Now zinc and cadmium both crystal- 
lize in the so-called ‘hexagonal close-packed 
system,” and have axial ratios of 1.86 and 1.89 
respectively, whereas a system of close-packed 
geometrical spheres would have an axial ratio of 
1.633, so that Zn and Cd both have an abnormal 
separation of the atoms along the direction of 
the c axis of the crystal. The interesting thing is 
that the atomic vibrations are greatest along 
the direction in which the atoms are furthest 
apart. Brindley and Ridley,*® in fact, believe 
that the asymmetry of these crystals is due 
wholly or in part to these asymmetric atomic 
vibrations. However, it has been repeatedly 
observed in this laboratory that the axial ratio 
changes by only a very small amount when the 
crystals are cooled to the temperature of liquid 
air, at which temperature the atomic vibrations 
are only about 3 as large as they are at room 
temperature, so that the atomic vibrations may 
be said to have very little to do with the ani- 
sotropy of the crystal. It would seem that the 
asymmetry of the atomic vibrations is an effect 
rather than the cause of the lattice asymmetry. 
We may picture (rather crudely) the potential 
function along a given direction, such as the c 


1G. W. Brindley, Phil. Mag. 21, 790 (1936); G. E. M. 
Jauncey and W. A. Bruce, Phys. Rev. 50, 408 (1936); 
E. O. Wollan and G. G. Harvey, Phys. Rev. 51, 1054 
(1937). 

2G. W. Brindley, Proc. Leeds Phil. Soc. 3, 200 (1936). 

3G, W. Brindley and P. Ridley, Nature 140, 461 (1936). 


axis, which would exist if the atom under 
consideration were removed from the crystal 
without altering the position of any of the other 
atoms, as being a U shaped curve (see Fig. 1a) 
with the bottom of the U at the normal rest 
position of the atom that has been removed. 
The last requirement is necessary in order that 
the crystal be stable. If now, the crystal is 
stretched out in the direction of the c axis, as 
it is in Zn and Cd, the whole curve would be 
stretched horizontally (see Fig. 1b), the bottom 
of the U would become more nearly flat and for 
a given thermal energy the atom would execute 
larger vibrations than before. In the case of 
close-packed geometrical spheres (axial ratio 
c/a=1.633) the distance from any atom to 
each of its nearest neighbors is the same, and 
thus one would expect the atomic vibrations to 
be the same in all directions. According to 
Brindley‘ this seems to be very nearly true for 
Mg (c/a=1.624). 

The purpose of the present research was to 
see if this idea could be extended to crystals 
having an axial ratio less than 1.633. The 
elements beryllium, titanium, zirconium, haf- 
nium, yttrium, ruthenium, and osmium have 
axial ratios between 1.58 and 1.59. In addition, 
certain alloys have low axial ratios, and in 
particular the alloy of Zn and Cu which contains 
15 atomic percent of Cu has an axial ratio® of 
1.556. This alloy has the advantage that it 
differs only slightly from pure Zn since the 15 
percent of Zn atoms (atomic No. 30) are replaced 


*G. W. Brindley, Nature 138, 290 (1936). 
5 E. A. Owen and L. Pickup, Proc. Roy. Soc. 140, 179 


(1933). 
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by Cu atoms of atomic No. 29. Thus we have a 
crystal which is essentially a zinc crystal com- 
pressed along the c axis, and if our hypothesis is 
correct, the atomic vibrations should be greatest 
in a direction normal to the ¢ axis. 

For an isotropic crystal (in which case the 
amplitudes of atomic vibrations are equal in all 
directions) one gets a smooth curve when the 
atomic structure factors, fr, are plotted against 
(sin 6)/A, where @ is the Bragg angle for radiation 
of wave-length \. Since the fr are proportional 
to the square roots of the reflection intensities, 
and since the reflection intensities decrease as 
the amplitudes of atomic vibration normal to 
the reflecting plane increase, it appears possible 
to detect the relative amplitudes of vibration 
along and perpendicular to the c axis of the 
crystal by finding the fr values. Thus, if our 
hypothesis is correct, we should expect the fr 
values for reflections from planes making a 
small angle with the c axis to give lower values 
than those from planes making a large angle 
with the c axis. It will be seen in Fig. 2 that this 
is true, and hence the question of the relative 
amplitudes of atomic vibration might seem to 
have been answered. The theory of the tempera- 
ture effect on reflection intensities (see Eqs. 
(1), (3.0), (6)) shows, however, that the temper- 
ature effect should decrease with (sin @)/Xd. It 
will be noticed in Fig. 2 that just the reverse 
is true. This means that another effect is coming 
into play. What this effect is will be left for 
later discussion, but in any event fr values do 
not prove conclusively whether an asymmetry 
in the atomic vibrations exists or not. 











Fic. 1. V is the potential of the mth atom along a given 
line of atoms in a crystal. AV represents the thermal 
energy of the atom due to vibrations along this line, and 
u, u’ the amplitudes of vibration. d and d’ are the distances 
between atoms for crystals a and b, respectively. 
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Fic. 2. Atomic structure factor of Zn 85 percent, 
Cu 15 percent (atomic percent) vs. (sin @)/X. 


Since the fr values gave inconclusive evidence 
of atomic vibration asymmetry, it was decided 
to try the method used by Wollan and Harvey! 
to find the vibration asymmetry in zinc. This 
method consists in comparing the reflection 
intensities at two different temperatures. If it 
is assumed that whatever other causes may make 
the reflection intensities greater or smaller, their 
temperature dependence still follows the theories 
of Debye-Waller, and Zener, then it is possible 
to determine the quantities a7 and by in Eq. 
(6.0). Since a7 and br are proportional to the 
mean squares of the atomic vibrations along and 
perpendicular to the c axis, respectively, it is 
possible to test the atomic vibration asymmetry 
hypothesis independent of other effects. 


2. DESCRIPTION OF APPARATUS 


In determining the atomic structure factors, 
a metal, self-rectifying, x-ray tube designed by 
Dr. E. O. Wollan and the author was used in 
conjunction with a zinc crystal monochromator. 
Cu Ka rays were selected and allowed to fall on 
the powdered alloy at the center of a circular 
camera of radius 4.5 cm. The total length of the 
path of the x-rays from target to film was about 
14 cm. With this arrangement it was possible 
to get good photographs in from three to five 
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hours, and the use of a crystal monochromator 
has the advantage that there is almost no 
divergence of the incident x-rays and hence no 
focusing, which means that the lines are all of 
comparable width. This makes it possible to 
microphotometer all the lines on one film instead 
of having to microphotometer different sets of 
Debye-Scherrer lines on different films and then 
match the overlapped portions. The mono- 
chromatic rays from a crystal also give much 
less background fogging on the film, and the 
background is far more uniform than it is when 
filtered radiation is used. This latter fact makes 
it possible to determine the background with 
considerably more precision than is possible 
when filtered radiation is used. 

The apparatus used for the temperature effect 
work has been described in a paper by Wollan 
and Harvey.' In this case two nickel filters were 
used to get the Cu Ka rays. It was found 
advisable in all cases to clean the surface of the 
target frequently as this improves the relative 
amount of characteristic radiation given out by 
the tube. The background on the film is reduced 
and the exposure time is shortened. The x-ray 
outfit was supplied with current from a “constant 
voltage generator’ and the high voltage and 
emission current were manually adjusted to a 
set value every five minutes during the exposure 
in order to get constant radiation output. In 
half of the cases, the low temperature exposure 
was made first so that any diminution in x-ray 
intensity due to progressive coating of the target 
would average out when several films were used. 
The microphotometer used in measuring the 
line intensities on all films consisted of a caesium 
cell and a DuBridge balanced-bridge d.c. ampli- 
fier to operate a short period Moll galvanometer, 
together with a mechanical system for scanning 
the film and moving the sensitized paper by 
proportionate amounts. 


3. EXPERIMENTAL PROCEDURE 


The alloy was prepared by mixing powdered 
pure Cu and Zn in proper proportions to give an 
alloy containing fifteen atomic percent of Cu. 
The mixed filings were sealed in a quartz tube 
and heated to 600°C for a period of eleven hours, 
and then quenched from the molten state in 


ROBERT A. 


HOWARD 


cold water. This gives an alloy in the e phase of 
the Zn-Cu system with a melting point of 500°C.°® 
After quenching, the alloy was annealed for 
twenty-four hours, gradually heated to 400°C, 
slowly cooled to 200°C, and then allowed to cool 
in air. Filings made by a fine file were sealed in 
an evacuated Pyrex tube and annealed for 
thirteen hours at 300°C. The Debye-Scherrer 
lines obtained with this powder had too many 
spots on them for intensity measurements, so 
the filings were ground in an agate mortar with 
a small amount of powdered KCI. After one and 
a quarter hours of grinding the KCI was dissolved 
away with water and the alloy powder dried 
but not reannealed. The average particle size 
was about 2X10-* cm after grinding, which is 
only slightly smaller than before grinding. The 
chief function of the grinding seemed to be to 
cold-work the metal particles, as the spots on 
the Debye-Scherrer lines were now small enough 
that they did not interfere with intensity 
measurements. Owen and Pickup’ give the value 
1.556 as the axial ratio of this alloy. By plotting 
the experimentally determined values of (X/sin @) 
for the alloy used here on Hull and Davey’s’ 
logarithmic scale and then using their chart, it 
was confirmed that the crystal belonged to the 
hexagonal close-packed system and had an axial 
ratio of 1.56 as nearly as it could be determined. 

In finding the atomic structure factors, several 
Debye-Scherrer photographs were made. Each 
film had a set of intensity steps made on it so 
that any difference in the film or its development 
would not affect the determination of the relative 
intensities of the lines. The intensity steps were 
made by placing the film, held in a lead cassette 
with a rectangular aluminum window at one end, 
behind a rotating, step-cut wheel and allowing 
nickel filtered x-rays from a copper target to 
fall on it. The wheel was so cut that each suc- 
cessive strip on the film received the same 
increment of radiation energy. The relative in- 
tensities were measured by microphotometering 
the films and measuring the areas under the 
curves obtained for the various lines. After the 
intensities had been corrected for absorption in 


6 International Critical Tables, Vol. 2, p. 435. 
7A. W. Hull and W. P. Davey, Phys. Rev. 17, 549 
(1921). 
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ATOMIC V 


the briquet from Brindley’s® curve, the values 
obtained were multiplied by a suitable constant 
so that the corrected relative intensities from 
each film were on the same basis. The relative 
atomic structure factors, fr, were then calculated 
from the formula 


I/A =const. X pS*fr*o(), (1) 


where p is the multiplicity factor, S the crystal 
structure factor (=1, 2 or (3)! for hexagonal 
close-packed lattices), and 

1+ cos? 28 cos? 20 
$(0) = . 








———— 


sin @ sin 20 


where @ is the Bragg angle for the reflection and 
6 is the angle which the incident radiation makes 
with the monochromator crystal face. In this 
case cos? 28 = 0.649. 

In the temperature effect work, two exposures 
were made on the same film, one at room tem- 
perature and one at the temperature of liquid air. 
Intensity marks were made on each film as 
before. Nothing was changed between the two 
exposures except the position of the film. The 


8G. W. Brindley and F. W. Spiers, Proc. Phys. Soc. 

46, 841 (1934), give a set of curves for the quantity 
sin (20— 

y - 5b (which is the fraction by which the 
intensity is reduced by absorption in the briquet) as a 
function of the Bragg angle @ fcr several values of a, the 
angle which the incident radiation makes with the face of 
the briquet. 
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angle @ changes so little because of lattice con- 
traction, that changes in the absorption factor, 
A, may be neglected. Only the 2131 and 1015 
lines were used in this part of the work. The ratio 
of the intensity at liquid-air temperature to that 
at room temperature was determined for each 
line, in some cases by measuring the areas under 
the microphotometer traces, sometimes by com- 
paring the heights of the two curves, and occa- 
sionally by both methods. The sums of the ratios 
obtained by the area method and the height 
method were 23.52 and 23.04, respectively, for 
the sixteen cases where measurements were made 
by both methods. Each line was micropho- 
tometered from one to nine times and the average 
value used. The mean values of the ratios for 
these two lines were then substituted in the 
Debye-Waller and Zener formulas (Eqs. (3.0) 
and (4)) and the values of ar and br calculated. 


4. RESULTs 
Table I shows the relative atomic structure 
factors obtained, and Fig. 2 is a plot of these 
values against (sin @)/X. 
Debye has shown that for a given reflection, 


Ip=Ice-™7, (3.0) 


where /, is the intensity of the reflection, J» is the 
intensity which one would get if the crystal were 
at the absolute zero of temperature, and Mr, is a 
function of the temperature. Let 7,4 =intensity 


TABLE I. Atomic structure factors. 


| 
Line | sin@/Al 1-19 | 4-20 | 4-39 | A-40 | 1-41 | 
1010 0.207 | 90° 9.10| 8.50| 
0002 | .234/ 0° 21.0 | 22.6 | 19.5 
1011 | .238| 60° 55’ | 46.7 47.0 
1012 | .312| 41° 55’| 6.12] 6.08| 6.62 
1120 | .363| 90° 5.73| 4.84] 5.73] 5.00 
1013 | 408) 31°41’ | 8.68] 8.56] 7.83] 7.51 
2020 | .417 | 90° 
1122 | .432| 57° 19’ 8.01 
2021 | .434| 74° 31’ 3.98 
0004 | .468/ 0° 1.47 
2022 | .478| 60° 55’ 1.45 
1014 | .512| 24° 5’ 1.43 
2023 | .547| 50° 8’ 
2131 | .566| 78° 7’ 4.92| 4.79 | 4.46 
1124 | .592| 66° 40’ 4.51| 4.99 | 4.44 
2132 | .601| 67° 10’ 1.98| 2.24 | 2.25 
1015 | .621| 19° 37’ 4.32 





























| 


| RELATIVE INTENSITIES CORRECTED FOR ABSORPTION (//A) } | 


——— | 
| A-48 | 4-49 | 4-50 | psrece) | fr 


A-46 fo 

43.9 | 15.2 41.73 | 16.5 

44.9 | 23.2 +0.32 | 25.2 

194. | 16.8 +0.16 | 18.6 

7.00 8.74 | 10.6 35.1 | 15.7 40.50 | 18.5 

5.36| 5.73 | 5.53 | 5.90 50.7. | 11.2 40.10 | 14.3 

0.2 | 7.98 | 7.30 | 7.97 59.7 | 12.6 +0.16 | 16.7 
1.44(?) 9.66 

7.18] 6.93 | 6.06! 6.11 73.4 | 10.3 40.18 | 14.3 

4.08} 3.56 | 4.22 | 4.40 54.7 9.2340.11 | 13.1 

1.22| 1.40 | 1.37 | 1.26 11.2 | 11.8 40.14 | 16.5 

1.62 | 1.76 | 1.70 16.4 | 10.6 +0.16 | 16.0 

1.50] 1.73 | 1.62 | 1.28 16.3 | 10440.19 | 15.8 
3.07 | 3.37 | 3.06 oy 

4.51| 4.50 | 4.33 | 4.13 114. 6.75+40.05 | 12.2 

4.29] 4.77 | 4.17 | 4.32 90.7 7.57+0.06 | 14.3 

2.41] 2.15 | 2.11 | 2.37 49.2 7.224+0.06 | 14.0 

4.19| 4.27 | 4.25 | 4.48 96.8 7.18+0.03 | 13.2 























See SE 
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of reflection when crystal is at temperature of 
liquid air, and 7pr=intensity at room tempera- 
ture. Then 


Tpa=Le "44, (3.1) 
Trr=1e?™ 87, (3.2) 
R= 4/Trr=eOM RTM 1a), (4) 
and 
log R=2(Mrr—M za). (5) 
From the equations derived by Zener,’ 
Mr=(ar cos? Y+br sin? y)(sin 6/d)?Q(x), (6.0) 


where y is the angle which the normal to the 
reflecting plane makes with the C axis of the 
crystal, Q(x) is the quantization factor and is 
nearly equal to unity, x= ©/T7, © is the charac- 
teristic temperature of the metal, 7 is the 
absolute temperature, and adr and br are func- 
tions of the temperature only. 

For a cubic crystal of atoms of one kind, 
Debye and Waller have shown that 


sin 0\? 
(=) | | 


Since Mr is a linear function of 7 for cubic 
crystals, we shall assume that it is also a linear 
function of 7 for anisotropic crystals. Letting 
ar=a-T, and by=)b-T, we get 


Mr=Q(x)T(sin 6/d)*(a cos? y+ sin? y). 


6h?7 
Mr=—— 
mk or 


~ 
— 


(6.1) 


Actually © is a function of both y and 7, but 
this has such a small effect on the value of Q 
that © is assumed to be constant. Moreover, 
© is not known for this alloy, but it can be taken 
equal to that of pure zinc to a sufficiently good 





approximation. Combining Eqs. (5) and (6.1), 
we get 

log R ; ; 
——— = 2(QrrTrr—QraT 1a) (a cos? p 
(sin 0/d)? 

+5 sin? y) =2A(QT)(a cos?, y+6 sin? y). (8.0) 
For the 2131 line, 

Terr log R, 

ar cos? ¥,+7 sin? y= , (8.1) 


2A(QT)(sin 6:/d)? 


°C. Zener, Phys. Rev. 49, 122 (1936). 
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Fic. 3. Relative structure factor for an atom at T=0. 


and for the 1015 line, 
Terr log Ro 


. (8.2) 
2A(QT)(sin 62/2)? 





ar cos? Yo+br sin? yo= 


The results of fifteen films and seventy-four 
microphotometer traces are used in obtaining 
solutions for ar and br. We get 


ar=1.51+40.118A? 


and br =1.84+0.098A°. 


DISCUSSION 


The individual values of R; and Re vary over 
a considerable range; so that the values obtained 
for ar and br can have no high accuracy. Two 
facts are to be noted, however. In the first place, 
the probable errors are sufficiently small that the 
limiting values of ay and br do not overlap. 
Secondly, if one uses (a7—br) as a measure of 
the anisotropy as Jauncey! prefers to do, he gets 


(ar—br) = —0.3340.155A2, 


so that the limiting values of the anisotropy are 
always negative. Therefore, while the results 
cannot be said to prove absolutely that dr is 
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Fic. 4. Atomic structure factor of Cd 80 percent, 
Ag 20 percent (atomic percent) vs. (sin @)/X. 


greater than a7, it seems extremely probable that 
this is the case, and our hypothesis about the 
vibrations being greater in a direction normal to 
the c axis of the crystal than along it is ap- 
parently confirmed. 

If we now use a7 and br to calculate the rela- 
tive structure factors reduced to an atom at 
absolute zero temperature (no correction for 
dispersion) we get the results shown in Fig. 3. 
It will be noticed that the results are not 
essentially different from those for fr, and the 
question immediately comes up, ‘“What causes 
the f values to vary so violently as a function of 
the orientation angle, y?” 

The thing which suggests itself first is distor- 
tion caused by grinding the crystals. Elastic 
distortion we can rule out immediately because 
the effect of elastic distortion is to broaden the 
lines, and the lines are no broader than those 
obtained from annealed filings of the alloy or 
from pure sublimated zinc. Barrett,!® and 
Hengstenberg and Mark" have also found that 
cold working a metal has very little effect on 


the line width. Plastic deformation, in which 


© Barrett, Phys. Rev. 47, 333 (1935). 
eee and Mark, Zeits. f. Physik 61, 435 
(1930). 


slipping takes place along certain planes in the 
crystal, could give rise to highly strained surface 
layers at points along the slip planes or at 
fractures, and possibly even to a layer of 
amorphous material between the slip planes. 
Amorphous material along the slip planes would 
introduce an additional absorption factor which 
would be a function of both @ and y. It is also 
possible that some part of the scattering of the f 
values is due to electron asymmetry in the atom. 
Such an idea is consistent with the fact that the 
spread of the points decreases with increasing 
(sin @)/A, which is just what is to be expected if 
the electron distribution in the atom is asym- 
metric. If the whole effect is due to electron 
distribution asymmetry, the effect must neces- 
sarily reach down into the inner electron shells, 
as an asymmetric distribution of the outer two 
electrons of zinc is insufficient to give the ob- 
served change between the f values for y=90° 
and y=0°. 

Figure 4 shows some unpublished fr values 
obtained by Dr. E. O. Wollan® for a Cd-Ag 
alloy containing approximately twenty atomic 
percent silver. This alloy also has an axial ratio 
of 1.56. The situation is almost exactly parallel 
to that of the Zn-Cu alloy, in that we replace 
part of the atoms of an element whose crystals 
have a high axial ratio (1.89) by an element 
whose atomic number is one less. The alloys 
were both prepared in the same manner, and 
the powders treated in the same way. It will be 


noticed that here again the spread of the points 


increases as (sin 6)/X decreases, although the 
percentage change is somewhat less. Whatever 
the explanation of this large change is, it is 
evident that it will suffice for both alloys. 
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The characteristic temperature for MgO has been directly measured by the x-ray method 
and a value of @=743°K for the range between room temperature and the temperature of 
liquid air obtained. A value of © was also obtained with the aid of theoretical fo values and 
the result so obtained was found to be in qualitative accord with the above value. In this 
connection it was noticed that the theoretical fo values for MgO are considerably in error. 
The characteristic temperature determined by the usual method from specific heat data 
shows a considerable dependence on temperature but the mean of these values over the range 
under consideration is in good accord with the x-ray value. 





INTRODUCTION 


HE theory of the effect of temperature on 
the intensity of reflection of x-rays from 
crystals which was developed by Debye and later 
corrected by Waller involves a quantity which is 
referred to as the characteristic temperature of 
the crystal. This quantity was introduced in con- 
nection with the theory of specific heats, and in 
the Debye form of this theory the characteristic 
temperature is defined by the relation 0 =h»v,,/k 
where v,, is the limiting frequency of the atomic 
vibrations considered propagated as _ elastic 
waves. Some uncertainty still exists, however, as 
to the conditions under which a value of the 
characteristic temperature calculated from spe- 
cific heat data can be used in the Debye-Waller 
theory for the determination of the effect of 
temperature on the intensity of x-ray reflection. 
Measurements which bear on this point have 
been made by James and Firth! on NaCl, by 
James and Brindley ? on KCl, by James, Brindley 
and Wood 3 on Al and by Shonka ‘ on NaF. The 
theory of the temperature effect is strictly valid 
only for simple cubic crystals containing atoms 
of one kind and one sees that the above examples 
represent varying degrees of departure from these 
conditions. The experiments to be described in 
this paper deal with measurements of the effect of 
temperature on the reflection of x-rays from MgO 
1R. W. James and E. Firth, Proc. Roy. Soc. Al17, 62 
(1927). 


2R. W. James and G. W. Brindley, Proc. Roy. Soc. 
A121, 155 (1928). 

3R. W. James, G. W. Brindley and R. G. Wood, Proc. 
Roy. Soc. A125, 401 (1929). 

‘J. J. Shonka, Phys. Rev. 43, 947 (1933). 


and a comparison of the characteristic tempera- 
ture so obtained with that which one gets from 
specific heat data. 


THEORY 


Although the Debye-Waller theory is strictly 
valid only for crystals consisting of atoms of one 
kind, we shall assume for the present that it can 
be applied to the case of MgO and later we shall 
discuss the limitations of this assumption. Ac- 
cording to this theory the atomic scattering fac- 
tor at a temperature T is related to the scattering 
factor fo for the atom at rest by the equation 


fr=foexp (—K sin? @/?*) (1) 
K=6h?TQ(x)/mkO* (2) 


1 7 tdE x 
Q(x) -|-f — +] (. 
xJy ef —-1 4 


is sometimes referred to as the quantization 
factor. Q(x) can be obtained from tables ° for all 
desired values of x= 0©/T. It is a quantity which 
departs appreciably from unity only for low tem- 
peratures (i.e., 7< ©). The exponent K sin? @/»? 
is commonly denoted by the symbol /, and the 
expression e~” is known as the Debye-Waller 
temperature factor. 

From Eq. (1) it is seen that K, and hence 0, 
may be evaluated from a knowledge of corre- 
sponding values of fr and fo. If theoretical values 


5R. W. James and G. W. Brindley, Phil. Mag. 12, 81 
(1931); p. 107. Note correction p. 729. 


where 


and 
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CHARACTERISTIC 


of fo have been calculated from the wave mechan- 
ics, experimental values of fr at one temperature 
will suffice to determine K. 

Let us first consider the more directly experi- 
mental determination of M from measurements 
of fr at two different temperatures. Since the 
intensity of reflection from powdered crystals is 
proportional to the square of the atomic scatter- 
ing factor we can write Eq. (1) in the form 


frye /fre? = py/pg=erMe-MD = g2AM (4) 


where AM is written for (M.—M,) and p; and pe 
represent the integrated intensities at the two 
temperatures. Taking logarithms of both sides 
we obtain 


log (p:1/p2) =2AM =2AK sin? 6/)? (5) 
or (\?/sin? 6) log (p:/p2) =2AK, (6) 
where 
2AK =2(K2—K;) 
= (12h? /mk ©?) | T2Q(x2) —T1Q0(x1)}. (5.1) 


Experimental values of the left-hand side of Eq. 
(6) afford values of 2AK from which © may be 
calculated. Alternatively, if in Eq. (5) we con- 
sider log (p:/p2) to be plotted vs. sin? 6/X? we may 
regard the quantity 2AK as the slope of the re- 
sulting straight line. 

From this same equation it is seen that the 
quantity which must be determined experi- 
mentally is the ratio p;/p2 of the intensities at 
two temperatures as a function of the scattering 
angle. In order to obtain an absolute value of this 
ratio it is necessary that the measurements at 
each temperature be made under identical condi- 
tions and that the total energy falling on the 
crystals in each case be the same. In the photo- 
graphic method of recording intensities this might 
be accomplished by taking identical exposures in 
the two cases. This method was used successfully 











PLANE | pi’/p2’ loge (e1’/p2’) sin? 6/d? 
200 | 41.024 0.024 | 0.057 
220 | 1.050 0.050 0.114 
222 | 1.030 0.030 0.171 
420 1.109 0.104 0.284 
422 1.070 0.068 0.341 
440 1.167 0.156 0.453 
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Fic. 1. Weighted mean of p,;’/p2’ as a function of sin® @/X?*. 


by Shonka ‘ in his measurements on Nak. The 
following consideration will however show that a 
knowledge of the absolute ratio p;/p2 is unneces- 
sary. If the exposure of a film at temperature 7, 
is some unknown constant c times as great as the 
exposure of another film at temperature 7» 
then the correct ratio of the intensities for equal 
exposures is p;’/cpe’ where p;’ and pe’ are the 
observed integrated intensities for the lines con- 
sidered. If we now substitute this correct ratio 
into Eq. (5) it can be written in the form 


log (p1'/p2") =2AK sin? 6/X+log c (7) 


and we see that the unknown constant c enters 
as an intercept and the slope 2AK of the line has 
not been affected. 


EXPERIMENTAL PROCEDURE 


The apparatus used in this research was the 
same as that used previously by one of us ® in 
connection with another problem. It consisted of 
a cylindrical camera the top of which formed a 
metal Dewar flask. The sample of powdered MgO 
in the form of a flat briquet was held against the 
inner copper wall of the Dewar so that the face 
of the briquet was always in the center of the 
camera. Copper Ka radiation was used and pic- 
tures were taken with the sample held at room 
temperature and at the temperature of liquid air. 
Seven pairs of films at the two temperatures were 
chosen for analysis and for each pair of films 
values of p;'/p2’ were obtained as a function of 
sin? @/)? for all lines of even order (h+k+/=even 
number). The odd order lines (h+k+/=odd 


*E. O. Wollan and G. G. Harvey, Phys. Rev. 51, 1054 
(1937). 
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Fic. 2. Values of Mevyen vs. sin? 0/d?. 


number) were not used for reasons to be discussed 
below. The data from each pair of films were then 
combined in such a way as effectively to reduce 
all the measurements to a common value of the 
unknown constant c. 

The results are shown in Fig. 1, where the 
weighted mean of the seven values of p;’/p2’ have 
been plotted. The slope of the line, which is a 
measure of the temperature effect, is seen to be 
very small. p,’ (for liquid-air temperature) is 
only sixteen percent greater than pe’ (for room 
temperature) for the largest value of sin @/) at 
which a line was obtained. The slope was deter- 
mined by the method of least squares to be 
0.299. 


THE CHARACTERISTIC TEMPERATURE 


From direct x-ray measurement 

On the basis of the simple theory we can now 
calculate the characteristic temperature of the 
crystal by substituting the slope which we have 
obtained experimentally into Eq. (5.1). Since we 
are dealing with a crystal containing two kinds of 
atoms we use an average value of the mass of 
these atoms. By this method we obtain the value 
© =743°K for the characteristic temperature. 

Since, however, MgO does not consist of 
atoms of one kind as the simple theory requires, 
let us consider what modifications we must make 
in the treatment. Waller has shown that in the 
case of such a crystal each kind of atom should 
be assigned its own temperature factor. For this 
reason, we must consider separately the even 
order spectra (reflections from planes containing 
both kinds of atoms in equal numbers) and the 
odd order spectra (reflections from alternate 
planes of each kind of atom). For even order re- 
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flections we thus have that 

Feven = Fi + Fo=fie~™'+foe-™ (8) 
and for odd order reflections 

Foaa = Fi — Fo=fie~™' —foe—™?, (9) 


where the small f’s refer to the atoms at rest and 
the subscripts 1 and 2 refer to Mg and O, respec- 
tively. We may then define an effective tempera- 
ture factor for the purposes of the simple theory 
by writing for the even reflections 


fiem™ +foe—™2 = (fi +foje—Meven (10) 
and for the odd reflections 
fe-M\— fue Mt=(fr—fa)e-Moas, (11) 


Since the M’s are small quantities we can re- 
place e~” by 1—M, obtaining 


Meven = (fi: MitfeMe)/(fitfe), (12) 
Moaa= (f1 Mi —feMe)/(fi—fe) (13) 


from (10) and (11), respectively. 

We see from this that Meyen is an average of 
M, and Mz weighted according to the relative 
scattering powers of atoms 1 and 2, while M,.aa 
is nothing like an average. Hence Moyen may 
legitimately be used for the unsubscripted M in 
the simple theory and e~™“eves considered as an 
TABLE II. Room temperature values of f for magnesium oxide 


(even orders only). 























i 
PLANE sin 0/X | WOLLAN | FROMAN | R. AND W. | MEAN 
200 0.238 | 13.40 13.25 13.32 
220 0.337 10.15 10.15 10.15 10.15 
222 0.413 8.20 8.19 8.29 8.23 
400 0.476 6.94 6.41 6.99 6.78 
420 0.533 5.92 5.83 6.19 5.98 
422 0.583 5.35 4.99 5.57 5.30 
440 0.673 4.00 3.68 3.84 
442+ 600 0.714 3.78 3.66 3.72 
620 0.752 3.22 3.07 3.14 
622 0.790 3.00 3.03 3.01 
444 0.825 2.90 3.22 3.06 
640 0.858 2.60 2.85 Be 
642 0.891 2.52 2.77 2.64 
800 0.952 2.21 2.71 2.46 
820+ 644 0.982 2.05 2.31 2.18 
822+ 660 1.010 1.86 2.26 2.06 
662 1.049 2.11 2.11 
840 1.063 2.03 2.03 
842 1.088 1.80 1.80 
664 1.117 1.80 1.80 
| 
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average temperature factor for the two kinds of 
atoms. .V/.uaa is definitely not suitable and hence 
we have excluded data from odd order reflections 
in calculating V7. It should be, however, possible 
from separate values of Moyen and M,aa to obtain 
M, and Me separately by the use of equations 
(12) and (13). We have not done this because the 
odd order reflections were too weak for sufficient 
accuracy. 

Our method of determining Moye, for MgO is 
based on the assumption that this quantity is 
of the form (const.) Xsin? 6/d? as is the WM of 
the simple theory. To test this point we have put 
room temperature values for the f’s in Eq. (10) 
since this introduces a negligible error and these 
f's are available in the literature (see Table IT). 
For M, and Mz we have taken Brindley’s esti- 
mates’ of 0.33 sin? 6/X* and 0.20 sin? 6/d*, re- 
spectively. Plotting Moyen vs. sin? @/* gives the 
very good straight line through the origin shown 
in Fig. 2, establishing that Meyer is of the required 
functional form. 


From theoretical f, values 

As was mentioned above the characteristic 
temperature can also be determined with the aid 
of Eq. (1) from measurements of fr at one 
temperature if theoretical values of fy are avail- 
able. Various values of the structure factor for 
even orders of reflection from MgO at room 
temperature are listed in Table II. The first 
column gives the absolute measurements of 
Wollan* who used an ionization method, the 
second gives the data of Froman® obtained by a 
photographic method and the third column gives 
the data which we have obtained in connection 
with this investigation. The values in the last 
two columns, being relative, have been fitted to 
those in the first column at the 220 reflection. 
The average of these data is plotted in the lower 
curve of Fig. 3. 

Values of the structure factor for the atoms at 
rest can be obtained from the work of James and 
Brindley.» These values were obtained by a 
method of interpolation in connection with solu- 
tions by Hartree’s method of self-consistent 
fields for neighboring atoms in the periodic table. 

*G. W. Brindley, Phil. Mag. 9, 1081 (1930). 


*E. O. Wollan, Phys. Rev. 35, 1019 (1930). 
* DP. K. Froman, Phys. Rev. 36, 1330 (1930). 
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Fic. 3. Upper curve, theoretical fy values. Lower curve, 
experimental fp values. 


The upper curve in Fig. 3 represents the even 
order structure factor for the atoms at rest as 
obtained from the James and Brindley paper. 
One is immediately struck by the fact that the 
two curves cross at about sin 6/A=0.3. This 
means that either the theoretical or the experi- 
mental results, or both, are in error. An inspec- 
tion of the results of Table II shows such good 
agreement between the various observers that it 
would seem unlikely that a large error could 
exist in the fr curve. On the other hand there 
are reasons to believe that the fy values as deter- 
mined for atoms not having closed shells may 
be considerably in error. It is known that meas- 
urements on the scattering of x-rays from oxygen 
and nitrogen in the gaseous form are also not in 
accord with the calculated fy values." 

In view of these facts this method of obtaining 
the temperature factor would seem to be rather 
unreliable. It seems likely, however, that the 
errors in fo should tend to become smaller as 
sin 6/X increases because the outer electrons 
contribute very little to the curve at large angles. 
In accord with this is the curve of Fig. 4 in which 
K, defined by the relation 


fr=fo exp (—K sin? 6/)*), 


is plotted against sin? 6/\*. The trend of the data 
make it seem reasonable to suppose that at large 
angles the curve will approximately asymptote 
to the solid line, which represents the charac- 
teristic temperature ©=743°K which we ob- 


10 A. H. Compton and S. K. Allison, X-Rays in Theory 
and Experiment (Van Nostrand, 1936), pp. 161-4. 
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peratures. 


From specific heat data 


The usual method of obtaining the charac- 
teristic temperature from specific heat data is to 
find that value of © for a given temperature 
which makes the observed value of C, equal 
that given by the Debye theory. Using experi- 
mental data for C, obtained from tables" we 
have calculated © for various temperatures. 
The results are shown in Fig. 5, the crosses and 
the circles representing data taken from two 
different sources. The agreement between the 
data is seen to be far from satisfactory. One 
notices also a considerable variation of © with 
temperature. In our determination of © from the 
x-ray data we have assumed that it does not 
vary appreciably with temperature over the 
range which we have used, and it can be shown 
that the x-ray results are inconsistent with a 
large variation of © with 7. It is, however, 
interesting to note that the mean of the specific 
heat values over the range of temperature used 
gives a value of ©=744°K which is in sur- 
prisingly good agreement with the value ob- 
tained from the x-ray measurements. Too much 
significance should not be attached to this fact, 
however, since both the x-ray and the specific 
heat values carry a rather large probable error. 

It should be mentioned that Zener and 


1 International Critical Tables V, pp. 91, 98. 
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Fic. 5. Values of @ from the experimental data for C, are 
calculated for various temperatures. 


Bilinsky” have given a theoretical justification 
for a certain amount of variation of © with T 
due to the expansion of the lattice. In the case 
of MgO this turns out to be of the order of only 
3°K which is negligible to our degree of accuracy. 


CONCLUSION 


A value of the characteristic 
equal to 743°K obtained by the x-ray method for 
MgO agrees well with the mean of the specific 
heat values over the range from room tempera- 
ture to the temperature of liquid air. It seems to 
vary inappreciably with the temperature. 

The difference between the amount of the 
temperature variation of the specific-heat 0 
given by the Zener-Bilinsky method and that 
given by the usual method implies that MgO 
fits the simple Debye model badly, since the 
former method arises from the application of 
straightforward thermodynamics to the Debye 
model of a solid. 

The wave mechanical fy values available at 
present for Mg and O are not satisfactory, 
especially at small angles. The values become 
more reliable at the larger angles. 


temperature 


2 C, Zener and S. Bilinsky, Phys. Rev. 50, 101 (1936). 
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Intensity Relations in the Low Voltage Helium Arc 


Joun J. HEILEMANN 
Randal Morgan Laboratory of Physics, University of Pennsylvania, Philadelphia, Pennsylvania 
(Received April 14, 1938) 


The linear amplifier photo-cell unit from a Razek- Mulder 
spect rophotometer is used with a high speed monochromator 
to measure the relative intensity of helium spectrum lines 
under various conditions of excitation; the discharge is 
developed in a copper ball 10.5 cm in diameter. A linear 
relation is found between arc current and line intensity for 
a series of voltages between the excitation potential and 
100 volts, and from this family of curves a family of 
excitation functions is drawn for a range of currents up to 
thirty milliamperes. The voltage intensity relation for eight 
helium lines (4388, 4438, 4471, 4713, 4922, 5016, 5876, 


6678) was determined directly. The results check those of 
Razek up to 50 volts and those of Elenbaas over the entire 
range, although secondary maxima are indicated in the 
region between the ionization potential and 45 volts. 
Additional experiments show that the cathode surface 
plays a part in determining the character of the discharge, 
and that the rise in Razek’s curves was caused by the 
discharge changing to the self-maintaining form; the 
hysteresis observed by him is eliminated by the use of a 
two-grid electron source. 





HE function which describes the relation 
between the intensity of a spectral line and 
the accelerating voltage of the exciting electrons 
is called the “optical excitation function”’ of the 
line. This function in helium has been the 
subject of considerable investigation.'~* If the re- 
sults of these investigations are to be interpreted 
in the light of the usual microscopic theories it 
is necessary to insure that (1) the electrons are 
homogeneous as to velocity, (2) each electron 
produces a single photon emission, (3) the 
number of electrons projected into the gas can 
be obtained from the measured current, and (4) 
the energy of the electrons is truly indicated by 
the voltage measuring device. These conditions 
will not be met if (1) the mean free path of the 
electrons is too short, (2) if large space charge 
effects are present, (3) if the atoms are excited 
by more than one collision, or by the combined 
absorption of energy from electrons and photons, 
(4) if the structures in the tube alter the char- 
acter of the discharge, (5) if there is self-absorp- 
tion of the radiation by the gas, or (6) if second- 
ary electrons are present. A linear relation 
between intensity and both current and pressure 
is usually accepted as evidence that the con- 
ditions enumerated above have been met. 
! Elenbaas, Zeits. f. Physik 59, 289 (1929). 
* Hughes and Lowe, Proc. Roy. Soc. A104, 480 (1923). 
’ Bazzoni and Lay, Phys. Rev. 23, 327 (1924). 
‘Hanle, Zeits. f. Physik 56, 94 (1929). 
* Michels, Phys. Rev. 36, 1362 (1929). 
* Thieme, Zeits. f. Physik 78, 412 (1932). 


7 Lees, Proc. Roy. Soc. 137, 173 (1932). 
8 Razek, Phys. Rev. 37, 1252 (1931). 


The photographic method of intensity meas- 
urement has been used by all previous investi- 
gators except Razek.* The greatest accuracy 
claimed for a photographic determination of 
intensity is 5 percent. However, Hanle* reports 
that variations on the plates used in his work 
were between 10 and 20 percent, consequently 
irregularities in the functions he investigated 
which are of an order of magnitude less than 10 
or 20 percent cannot be considered significant. 
Furthermore, the inherent slowness of the photo- 
graphic method detracts from its usefulness in 
excitation function determinations where the 
“fine structure’ of the phenomenon is being 
studied, since it is difficult to keep constant the 


” 


discharge tube currents and voltages for a 
sufficient length of time. Most of the curves in 
the literature for which data are given indicate a 
lack of smoothness, but the limited accuracy of 
the photographic photometry does not justify 
the drawing of anything but a smooth curve 
through the points. Furthermore, the number of 
voltages used is usually insufficient to make the 
curves which are drawn really indicative of the 
detailed character of the function. 

The spectrum tubes used by other workers*: * 
were usually of small dimensions and were 
furnished with nonequipotential cathodes, the 
inequality ranging from 0.3 volts to 5 volts. 
The envelopes were of glass and the electrodes 
and shields were large and placed unsymmetri- 
cally with respect to the discharge. It is well 

® Lay and Cornog, J. Opt. Soc. Am. 24, 149 (1934). 
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Fic. 1. Diagram of discharge tube for low voltage 
helium excitation. 


known in vacuum tube practice that the dis- 
position of the electrodes plays an important 
part in determining the tube characteristics, 
owing to the accumulation of surface charges and 
similar effects, consequently it might be expected 
that the behavior of a small discharge in a 
relatively large system of electrodes would be 
strongly influenced by the geometry of the system. 

Razek used a photoelectric spectrophotom- 
eter!’ |! the response of which was shown to be 
linear with intensity ; this instrument, with which 
the inherent slowness of the photographic 
method may be overcome, was shown by him to 
be applicable to the problem in hand. The linear 
amplifier from Razek’s spectrophotometer and a 
spectrum tube similar to his were used in the 
work here reported. From the description of this 
tube it will be seen that the discharge takes place 
in a ““hohlraum”’ relatively free from electrodes. 
The present work was undertaken with the idea 
of using this kind of tube to determine in greater 
detail the character of the excitation function, 
as well as to demonstrate that the sources of 
error listed above were not present at the rela- 
tively high pressures used in this type of equip- 
ment. The pressure of the helium was 0.27 mm 
of mercury. 

10 Mulder and Razek, J. O. S. A. and R. S. I. 18, 466 
(1929). 


1 Razek and Mulder, J. O. S. A. and R. S. I. 19, 390 
(1929). 


DESCRIPTION OF THE DISCHARGE TUBE 


The discharge takes place in a spun copper 
ball 10.4 cm in diameter, the inside surface of 
which is silvered and polished. A slit, over 
which is placed a tube and a glass plate, allows 
observation of the discharge. The electron source 
(Fig. 1) was mounted on the press from a com- 
mercial vacuum tube by spot welding to the 
leads a heater-type oxide coated cathode (diam- 
eter 1.7 mm) such as is used in the type '27 tube, 
a grid from the same type tube (diameter 4 mm) 
and a second grid from a type '24 tube (diameter 
7 mm). The grid nearer the cathode is the 
“control-grid’”’ and the other is the ‘‘accelerator- 
grid.”’ The wiring is shown in Fig. 2. 

The electrons are emitted by the cathode, and 
the number which is admitted into the space 
between G, and Gz is determined by the potential 
difference between G,; and C; the sense of this 
potential difference may be changed by means of 
the reversing switch. Thus the number of elec- 
trons available for acceleration by the field 
between G,; and Ge may be regulated without 
changing the temperature of the cathode; the 
cathode temperature is maintained by the usual 
insulated filament heated by a storage battery. 
A radio B battery supplies the potential differ- 
ence between G; and C, which potential controls 
the arc current; a bank of storage batteries is 
used to supply the accelerating potential between 
G,; and Gz. The discharge takes place in the 
equipotential space between the ball and G2, and 
the current in the arc is measured by the milli- 








= 


Fic. 2. Wiring diagram of discharge tube for low voltage 
helium excitation. 
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ammeter AG». For each observation in the course 
of the experiment a reading of each meter was 
taken so as to follow the behavior of the discharge 
tube as closely as possible. It was thought that 
a study of the relations between the various 
currents and voltages might yield information as 
to the function of the parts of the electron gun; 
as a matter of fact data were obtained (to be 
mentioned later) which apparently shed light 
upon the action of the control grid G,. 


DESCRIPTION OF EXPERIMENTAL SET-UP 


The discharge tube was permanently con- 
nected to the vacuum system through a stop- 
cock; between this stopcock and the ball was a 
charcoal trap which was immersed in liquid air. 
The initial pumping was continued for several 
hours, during which time the cathode material 
was activated. The presence of waxed joints 
prevented degassing, and the system was con- 
sidered clean when the electron current between 
the elements shcwed no luminosity and when the 
characteristics of the tube remained constant. 
Helium was admitted by the conventional 
arrangement of bulb and stopcocks, having been 
allowed to stand in two liquid-air traps for a 
half-hour before admission to the ball. Examina- 
tion of the discharge by means of a direct vision 
spectroscope showed no trace of impurity such as 
mercury or hydrogen. 

The large aperture (f: 1.4) double prism 
monochromator was designed by Dr. Joseph 
Razek and was constructed in this laboratory 
with funds provided by the Faculty Research 
Committee of the University of Pennsylvania. 
The entrance slit of this monochromator was 
placed as close as possible to the slit in the dis- 
charge tube, leaving just enough room for the 
sliding shutter 5S; (Fig. 3). The spectrum is 
formed on the exit slit 7, which may be moved 
laterally to admit any desired line to the photo- 
cell which, with its amplifier, is attached to the 
monochromator. The photoelectric current is 
amplified, and the amplified current is indicated 
by a galvanometer. A 6-8 volt headlight lamp 
provided a standard of intensity; its brightness 
was kept constant by the Richardson’s bridge 
arrangement indicated in Fig. 3 and described in 
detail by Razek.’ A piece of plate glass P 
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Fic. 3. Optical arrangement for comparison of helium line 
intensities with standard source. 


reflected a portion of the light into the system so 
that the lamp was virtually in the same position 
with respect to the collimator lens as the entrance 
slit Z. Opening the shutter S_ resulted in the 
formation of an jmage at J of the filament of 
the same color as the line whose intensity was 
being measured. 

The experimental procedure was simply to 
open shutter S; and to note the galvanometer 
current in the photo-cell amplifier, then having 
closed S;, to open S: and to make a similar 
reading for the standard intensity. The relative 
intensity plotted in the results is the ratio be- 
tween these two galvanometer currents. Such 
intensity determinations were made (1) for con- 
stant voltage and varying arc current and (2) for 
constant arc turrent and varying accelerating 
voltage. The constant voltage determinations 
were made to test the assumption that the arc 
current was a reliable measure of the number of 
collisions occurring in the field-free space and to 
test the reliability of the photometry. The 
constant current determinations gave the optical 
excitation functions of the lines under the con- 
ditions of the experiment. 

Previous workers have controlled the arc 
current by varying the temperature of the 
cathode, but in this work the distribution of the 
potentials on the grids was used as the con- 
trolling agent. This method was deemed ad- 
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cathode temperature. 


visable since it seemed that the form of the dis- 
charge changed with cathode temperature; 
“hot-spots”’ formed on the cathode as the result 
of operating at low temperatures and high 
accelerating voltages. Improved technique in the 
matter of applying and forming the cathodes 
resulted in such active coatings that above 
certain accelerating voltages the arc became 
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self-sustaining and continued to operate with the 
cathode cold. This tendency to become self- 
sustaining began to show itself in most cases 
when the accelerating voltage was in the neigh- 
borhood of sixty volts. 

Introduction of the control grid made it 
necessary to conduct a series of experiments to 
show that the intensity was not a function of the 
distribution of the potentials on the grids, but 
depended solely upon the current to the acceler- 
ator grid. Although the temperature of the 
cathode was varied it was possible to keep the 
current and voltage between G,; and G»2 constant 
by changing the control voltage between G, 
and C. Having decided upon a constant current 
and voltage, the heater current was varied so as 
to change the initial supply of electrons; the 
intensity of a line was measured as a function of 
cathode temperature. The results are shown in 
Fig. 4. The curves show that as long as the 
cathode was kept at a moderately high tempera- 
ture, the intensity remained constant, regardless 
of the control grid voltage or the initial supply 
of electrons, but an increase in intensity was 
noted when the cathode temperature decreased 
to the point where a self-sustaining arc began as 
evidenced by instability in the arc current. 
In the subsequent work, the cathode was always 
heated to a temperature well above that at 
which instability became evident. 


By making the control grid positive, the 
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Fic. 6. Family of voltage intensity relations derived from 
current intensity curves. 
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Fic. 7. Voltage-intensity relations in the low voltage helium 
discharge: lines 5016 (2'S—3!P) and 4438 (2!P—5'S). 


electrons arrived in the field-free space with the 
energy corresponding to the sum of the potential 
drops between the two grids, and the full 
spectrum always appeared when these two 
voltages added to 24.5, the ionization potential 
of helium. Characteristic curves were plotted in 
which the “‘breaks’’ were found to occur very 
close to the ionization potential, regardless of the 
distribution of the potentials on the grids. These 
experiments were considered as indicating that 
the voltages read on the meters were the voltages 
applied to the electrodes and that the intro- 
duction of the control grid gave rise to no 
significant complication. The drop across the 
milliammeters was read by the voltmeters, but 
this was less than the least count of the meters 
and could be neglected. The absence of contact 
potentials of notable magnitude was indicated 
by the agreement between the accepted value of 
the first ionization potential and the voltage at 
which the full spectrum appeared. The absence 
of voltages due to oscillations was ascertained by 
checking the readings of the voltmeters against a 
vacuum tube voltmeter; the check was within 
the precision of the instruments. 


RESULTS 


The linear relation between current and in- 
tensity is demonstrated by the curves in Fig. 5; 
the open circles represent data taken with in- 
creasing current and the solid circles with 
decreasing current. This linearity shows, as has 
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been mentioned, that the current to the acceler- 
ator grid is a measure of the number of electron 
collisions and that the photometry is reliable. 
By making the measurements at a number of 
voltages, the fact that the linear relation holds 
for all voltages between the ionization potential 
and 100 volts is demonstrated. A family of 
curves such as is shown in part in Fig. 5 may be 
used to plot a similar family of voltage intensity 
curves by simply using the ordinates corre- 
sponding to the constant current value of the 
excitation function. The curves shown in Fig. 6 
were drawn from the data supplied by the curves 
in Fig. 5. 

Although the linearity is demonstrated to well 
within the limits of precision of the measure- 
ments, the greatest deviations occur in the 
vicinity of 55 volts, where peak intensities of 
considerable magnitude have been reported by 
Cornog.” The curves indicate good linearity up 
to current values of thirty milliamperes, beyond 
which a slight tendency to saturation is noticed, 
as might be expected; this would indicate that 
currents as great as 30 ma may be used without 
changing the character of the excitation function. 

The excitation functions (Figs. 7-10) check the 
work of Razek as far as 60 volts; where Razek’s 
curves show an increase in intensity, the present 
results show either a decrease or a flattening. 
The increase in Razek’s functions is probably 
related to the effect noted above, viz., that at low 
cathode temperatures the arc tends to change its 
character ; as a result of this change the intensity 
increases. The collisions which produce the light 
whose intensity is being measured are pre- 
sumably in the positive column of the discharge, 
and should be caused solely by electron-atom 
encounters, but when the arc tends to become 
self-sustaining, other light producing mechanisms 
are involved, as is evidenced by the so-called 
“negative glow’ in the cold discharge. In the 
present arrangement, in which a dense cloud of 
electrons is maintained around the cathode, the 
positive ion bombardment of the cathode is 
minimized, consequently the negative glow 
should not be observed. The unusually long life 
of the cathodes in the two-grid electron source 
as compared to the one-grid type also indicates 
that destructive positive ion bombardment has 

2 Cornog, Phys. Rev. 32, 746 (1928). 
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been considerably reduced. The “hysteresis” 
effect observed by Razek has also been elimi- 
nated ; this is indicated in the curves, where the 
solid circles represent intensities measured during 
a run of decreasing voltages, while the open ones 
represent increasing voltages. It was in the hope 
of eliminating this hysteresis effect that the two- 
grid electron gun was developed, since it would 
then be possible to carry the tube through a 
voltage cycle without changing the cathode 
temperature and consequently without disturb- 
ing the condition of the cathode surface. That 
the cathode surface has something to do with the 
character of the discharge is evidenced by the fact 
that there was no hysteresis as long as the 
cathode was unchanged during a measurement. 

The action of the control-grid may be visual- 
ized by reference to the curves in Fig. 4. In this 
experiment the grid was used to keep the supply 
of electrons in the accelerating field constant; 
it will be noticed that the current to this grid is 
constant and that the grid serves to regulate the 
concentration-gradient between itself and the 
cathode so as to keep constant the electron 
density in the acceleration space. 

The voltage intensity curves are in good agree- 
ment with those of Elenbaas which show 
maxima in the vicinity of sixty volts as well as 
a dip around forty volts. This sixty volt maxi- 
mum was attributed by Lees’ to a change in 
the shape of the beam, and the inability of 
Elenbaas’ apparatus to integrate the intensity. 
The present apparatus is not open to this ob- 
jection since the silvered sphere does accomplish 
the integration, so that the sixty volt maximum 
seems to be real. Some significance may attach to 
the fact that Elenbaas’ apparatus is in general 
similar to that used in the present work; the 
cylindrical anode with a slit for observation 
constitutes a kind of “hohlraum”’ as does the 
sphere in the present work, and the entire 
current to the anode is also measured by Elen- 
baas. Although the results given here are only in 
general agreement with the results of the other 
workers,‘ *. 7 it will be noticed that a tendency 
toward an inflection at forty volts is exhibited by 
the data; the curves are drawn smoothly how- 
ever, since the precision of the photographic 
method is not great enough to warrant attaching 
significance to the indicated irregularities. There 
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is also in the present results a marked indication 
of secondary maxima between the ionization po- 
tential and about 45 volts; this is most noticeable 
in the lines 4471 (28P —4*D) and 5876 (2°P —3*D) 
(Fig. 10). A somewhat similar effect has been 
reported for mercury by Siebertz;" there is a 
possibility, however, that this effect in mercury 
is of a different order of magnitude from the 
effect reported here. 

The sudden changes which occur in the 
neighborhood of the ionization potential make it 
rather difficult to set the currents and voltages 
simultaneously, consequently no attempt was 
made to set the current at the indicated constant 
value, but intensities corresponding to the indi- 
cated value could be calculated on the basis of 
the experimentally proved linear relation be- 
tween current and intensity. Beyond about 35 
volts these corrections amounted to less than the 
precision of the measurements. 


CONCLUSIONS 
From these results it may be concluded that 


the discharge tube used by Razek, modified so as 
13 Siebertz, Physik. Zeits. 31, 141 (1930). 
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to control the supply of electrons without 
changing the cathode temperature, gives excita- 
tion functions which check with those determined 
by other workers who used the photographic 
method. The hypothesis that the “hysteresis” 
effect noticed by Razek is associated with 
changes in the condition of the cathode surface 
is substantiated, since keeping the cathode at a 
constant temperature and protecting it from 
positive ion bombardment eliminated the effect. 

The linear relation between current and in- 
tensity shows that pressures as high as those 
used do not lead to complications in the large, 
high intensity source used in this work. 

In addition to checking the form of the 
excitation functions determined by Elenbaas, the 
present work indicates secondary maxima in the 
region between the ionization potential and 
forty-five volts. 

The writer wishes to express his thanks to 
Professor C. B. Bazzoni, under whose guidance 
the work reported here was done; to Mr. F. R. 
Banks, who assisted in taking the readings, and 
to Dr. Joseph Razek for valuable aid and 
suggestions. 
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The Infra-Red Absorption by C'O,'° at 4.375 u 


Atvin H. NIELSEN* 
Mendenhall Laboratory of Physics, Ohio State University, Columbus, Ohio 


(Received April 22, 1938). 


The fundamental vibration-rotation band v; of C'80,'* has been found with its center at a 
wave-length of 4.3754 or 2284.5 cm™'!. The individual lines in the band have been fitted by the 
formula v3 =2284.5+0.780 N—0.0027 N? where N takes the values —2, —4, —6, ---, etc. for 
the P branch and +1, +3, +5, +---, etc. for the R branch. On applying the appropriate mass 
correction to v; and to the anharmonic constant x3; given by Adel and Dennison, these constants 
take values of 2295.8 cm™! and —11.4 cm™, respectively. Together these constants give the 
position of the center of the band as 2284.4 cm™, which is virtually in perfect agreement with 


observation. 


I. INTRODUCTION 


EASUREMENTS in the infra-red region 
near 4.374 in the atmosphere made by 
earlier investigators in this laboratory have 


* Assistant Professor of Physics, University of Tennessee, 
Knoxville. 


shown the existence of a weak absorption 
maximum with resolvable rotation lines on the 
low frequency side of the very intense funda- 
mental vibration v3; of the CO. molecule. Since 
part of this absorption region overlaps with the 
P branch in v3 of COz and the lines appear to 
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Fic. 1. The fundamental vibration-rotation band pv; of C’QO,.!*, 


have the same frequency separation! it seemed 
plausible that it might be attributed to an 
upper stage transition in the energy of the CO, 
molecule. The existence of this band was there- 
fore considered sufficiently interesting to warrant 
a thorough study of this region. 


II. EXPERIMENTAL 


The region of the infra-red near 4.374 has been 
very carefully mapped with the aid of a prism- 
grating spectrometer equipped with an echelette 
grating of 3600 lines per inch ruled for this 
region by R. W. Wood, and described in previous 
reports from this laboratory. The slits were 
narrowed down to include a spectral interval of 
about 0.4 cm~', and deflections of the galva- 
nometer were read at intervals of 0.35 cm. 
Fig. 1 shows the P branch, the center, and a 
few lines in the R branch of the band. The 
number of lines in the R branch is limited be- 
cause of overlapping with the P branch of the 
CO, fundamental v3. It may be seen that these 
lines are all resolved very well and the intensities 
in the P branch follow the usual regular envelope 
for such bands. The lines are numbered from 
the center outward. Every other line is missing 
because of the zero spin of the oxygen nucleus. 


1 P. Martin and E. F. Barker, Phys. Rev. 41, 291 (1932); 
Donald Cameron and Harald H. Nielsen, Phys. Rev. 53, 
246 (1938). 


III. Discussion 


It was possible to write a formula which would 
predict the frequencies of the lines in this band 
almost perfectly. Table I gives a list of the lines 
with the appropriate ordinal numbers, the ob- 
served positions in cm~', and the calculated 
positions in cm~ as obtained from the formula 


v3 = 2284.5+0.780N —0.0027N?; 


in which NV takes the values of —2, —4, —6, —8, 
—---, etc. for the P branch and +1, +3, +5, 
+9, +---, etc. for the R branch. It may be 
seen from Table I that in only one case in the P 
branch is the disagreement between the observed 
position and calculated position more than 
0.2 cm. It may also be seen from Table I 
that the separation of the lines measured in cm~! 
is the same as that in v3 for CO. which makes it 
almost certain that the band belongs to the CO, 
band system. Although the population of the 
upper energy levels is small, it was first thought 
possible that the band was the result of an 
energy transition from V;=2 to V3;=3. The 
center of this band may be predicted by the 
use of the constants »3;=2362.8 cm and 
X33 = —12.7 cm™ derived by Adel and Dennison? 
and lies at 2299.2 cm. Since the observed 
center of the band is at 2284.5 cm™ the differ- 


2A. Adel, and D. M. Dennison, Phys. Rev. 43, 716 
(1933) and Phys. Rev. 44, 99 (1933). 
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ence between the observed position and this 
calculated position is 14.7 cm~. 

It was felt that both the measurements and 
the theory were much more exact than this dis- 
crepancy indicates in view of the beautiful 
positional agreement of all the COs bands 
following the theory of Adel and Dennison. 
Because of this belief the identification of the 
band as an upper stage transition was abandoned 
in favor of its assignment as v3 of the isotopic 
molecule C"O,"*. Following this assumption the 
individual terms in the expressions for v3 and x33 
in Adel’s and Dennison’s paper were investigated 
for their dependence on the reduced mass of the 
molecule. On applying the appropriate mass 
correction for the isotopic molecule to the 


TABLE I. Frequency positions of the lines in v3 of C¥O,'*. 











N | OBSERVED »v | CALCULATED v | Av 
—44 | 2245.1cm™| 2244.9cem™| 40.2 cm” 
—42 | 2247.1 2246.9 | +0.2 
—40 | 2249.0 2248.9 | +0.1 
—38 2250.8 2250.9 | —@.! 
— 36 2252.8 2252.9 a * 
—34 2254.8 2254.8 | 0.0 
—32 2256.6 2256.7 | 0.1 
—30 2258.6 2258.6 0.0 
— 28 2260.6 2260.5 +0.1 
— 26 2262.4 2262.4 | 0.0 
—24 2264.2 2264.2 0.0 
—22 2266.0 2266.0 | 0.0 
—20 2267.8 2267.1 | +0.7 
—18 2269.7 2269.5 | +0.2 
—16 2271.5 2271.3 | +0.2 
—14 2273.2 2273.0 +0.2 
—12 2274.9 2274.7 +0.2 
—10 2276.5 2276.4 +0.1 
- 8 2278.2 2278.1 +0.1 
— 6 2279.7 | 2279.7 0.0 
-— 4 2281.3 | 2281.3 0.0 
- 2 2282.7 | 2282.9 —0.2 

0 2284.5 | 2284.5 
+ 1 2285.5 | 2285.2 +0.3 
+ 3 2286.5 | 2286.8 —0.3 
+ § 2288.1 | 2288.3 —0.2 
+ 7 2289.5 2289.8 —0.3 
+9 2290.9 2291.3 —0.4 
+11 2292.4 2292.7 —0.3 
+13 | 2294.1 2294.2 —0.1 
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individual terms these expressions become: 


u\' 3m heu\*? up 
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and 


3u Cp OC pf w} 
, 
x33 = -f — 
2) ; 


2a dw ww 4 





wy? —4ws3"(u/u’)- 


where » and yw’ are the reduced masses with 
C=12, and C=13, respectively. Substitution of 
Adel’s and Dennison’s constants w), we, ws, a, 0, ¢, 

etc. in the above expressions gives v3’ and 
x33’ the values 2295.8 cm=! and —11.4 cm“. 
These together give the position of the band as 
2284.4 cm™ and since the observed position is 
2284.5 cm~! the agreement is perfect. The evi- 
dence is considered conclusive enough to permit 
identification of this band as the frequency vs of 
the molecule CQO,"*. Since carbon is the central 
atom in COs, replacement by a heavier atom 
should not alter the moment of inertia, and this 
is verified by the fact that the observed spacing 
is the same as in‘the corresponding band for 
CXO,'6, 

In conclusion the author expresses his grati- 
tude to Professor David M. Dennison of the 
University of Michigan who very kindly con- 
sidered this problem and suggested the identity 
of the band reported here. The author also 
acknowledges his indebtedness to Professor 
Alpheus W. Smith, and to Professor Harald H. 
Nielsen who suggested the problem and without 
whose cooperation these measurements could not 
have been made. 
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The Energy Loss of Positive Electrons in Passing Through Aluminum* 


B. R. Curtis 
University of Michigan, Ann Arbor, Michigan 
(Received April 23, 1938) 


The energy loss of positive electrons between the ener- 
gies of 0.3 and 1.5 Mev has been determined in a hydrogen- 
filled cloud chamber for three thicknesses of aluminum. 
The positive electrons, which were obtained by bombard- 
ing iron with deuterons, were passed through aluminum 
absorbers of 0.0275, 0.053, and 0.114 cm thickness, and 
the values of the average energy loss were compared with 
the values predicted by the theoretical formula of Bloch. 
The value of the average energy loss for the 0.0275 cm 


INTRODUCTION 


HEN iron is bombarded with deuterons! a 
radioactive isotope attributed to cobalt is 
formed. This emits positive electrons of energies 
up to 1.6 Mev. The problem in this study is to 
determine the amount of energy lost by these 
particles in passing through various thicknesses 
of aluminum. At the same time a direct compari- 
son is made with negative electrons to ascertain 
whether or not there is any difference in the 
amount of energy lost by the two particles. 
Within the past few years a great deal of work 
has been done with negative electrons. The 
experiments show good agreement with theory 
for substances of low atomic number ; whereas for 
high atomic numbers the experimental values are 
much higher than theory would predict.?~* Ab- 
sorption measurements with positive electrons of 
low energy indicate that they lose energy in much 
the same manner as the negative electrons.® In 
the case of the high energy positrons found in 
cosmic radiation, the evidence seems to indicate 
that the amount of energy lost by the positive 
and negative electrons is the same,* although one 
experiment has seemed to indicate a difference.’ 
*A preliminary report of this investigation was given 
at the Chicago meeting of the American Physical Society 


in November, 1937. 

1B. T. Darling, B. R. Curtis and J. M. Cork, Phys. Rev. 
51, 1010 (1937). 

2 J.J. Turinand H. R. Crane, Phys. Rev. 52, 63, 610 (1937). 

3A. J. Ruhligand H. R. Crane, Phys. Rev. 53, 618 (1938). 

‘J. Laslett and D. Hurst, Phys. Rev. 52, 1035 (1937). 

§ J. Thibaud, Phys. Rev. 45, 781 (1934). 

*P. M. S. Blackett and J. G. Wilson, Proc. Roy. Soc. 
A160, 304 (1937). 

7 J. Crussard and L. Leprince-Ringuet, J. de phys. et 
rad. 8, 213 (1937). 


thickness was 0.123 Mev, which is roughly 20 percent 
greater than the theoretical value; whereas for the 0.053 
cm thickness, the loss was 0.274 Mev, an increase in 
difference to 40 percent. A direct comparison of the energy 
lost by the negative electron with that lost by the positive 
electron was made using electrons obtained from radio- 
active phosphorus. The results show that there is no 
essential difference in the amount of energy lost by the 
two particles. 


APPARATUS 


The cloud chamber used in these experiments 
has been described by Crane.* It was 15 cm in 
diameter and 4 cm deep, filled with hydrogen and 
ethyl alcohol vapor at a pressure of 100 cm of 
mercury. The magnetic field was obtained from a 
pair of air core solenoid coils. In the earlier 
experimental work, the current through the coils 
was turned on and off, but in the later experi- 
ments, the current was allowed to flow con- 
tinuously. This removed the possibility of any 
errors due to the contacts of the circuit breaker or 
to misreading the swing of the ammeter needle. A 
collimated beam of parallel light, obtained from a 
carbon arc, illuminated the center 1.5 cm section 
of the chamber. A diagram of the experimental 
arrangement is shown in Fig. 1. By locating the 
source S at the correct position, it was possible to 
obtain a partial magnetic separation of the 
positrons so that particles of only a certain 
definite energy range entered through the thin 
window in the wall of the chamber. This pre- 
vented the appearance in the pictures of too 
many extraneous tracks which did not strike the 
absorber at approximately normal incidence. 

The aluminum absorbers were coated with 
paraffin and lamp black and were placed across 
the center of the chamber parallel to the light 
beam. By measuring the radii of curvature of the 
incident and emergent tracks, the amount of 
energy lost by the particles in traversing the 
absorber could be determined. 

The method used for determining the curva- 


8H. R. Crane, Rev. Sci. Inst. 8, 440 (1937). 
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ture of the tracks was to project the full-sized 
image upon a white card and then to measure the 
radii of the tracks by means of a celluloid or 
cardboard template upon which circles of radii 
varying by 0.5 cm had been drawn. In order to 
be included in the data actually used, a track was 
required to satisfy the following conditions: 

(1) It must come from the window in the side 
of the chamber and must strike the absorber at 
approximately normal incidence (within 10°). Of 
course this criterion does not exclude the chance 
that an electron may be ejected from the absorber 
by a gamma-quantum and come back toward the 
window. 

(2) The track must be of sufficient length to 
make possible an accurate measurement of its 
curvature. Since only half of the 15 cm of the 
chamber is being considered, at least 6 cm of 
track must be measured. 

(3) The track must emerge from the absorber 
within 15° of the normal and must not make an 
angle greater than 10° with the plane of the 
chamber. Any track which satisfies condition (2) 
will satisfy this condition. 

(4) A track must not be visibly scattered. 
When working with particles of energy as low as 





thin window 











magnetic 
field colls 


=~ 
source —-(_] | Pb ] chamber light beom 


- - 





























Fic. 1. Experimental arrangement of cloud chamber and 
electron source. 
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Fic. 2. Energy losses plotted against the number of 
positive electrons, for 0.0275 cm aluminum and initial 
energies of 0.48 to 0.86 Mev. 


0.5 to 1.0 Mey, it is necessary to fill the chamber 
with hydrogen in order to keep down scattering 
to a minimum. The probability for scattering is 
inversely proportional to the square of the kinetic 
energy of the particle and directly proportional 
to the square of the atomic number of the 
scattering material. Thus if hydrogen is used, the 
scattering for these energies is very small. 


RESULTS 
Positive electrons 

The principal source of positive electrons was 
iron bombarded in the cyclotron with deuterons. 
These electrons have an upper limit in energy of 
1.6 Mev and a half-life of 18 hours. 

For the 0.0275 cm aluminum, approximately 
200 positrons were obtained which satisfied the 
necessary criteria. In Fig. 2 the momentum loss 
curve is shown for positrons ranging in initial 
energy from 0.48 to 0.86 Mev. The curve contains 
a total of 123 positrons. For convenience, the 
scale of energy loss is plotted below the Ap 
(change in radius of curvature) scale. To obtain 
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the change in momentum, the latter may be 
multiplied by the value of the magnetic intensity, 
H. The following values are recorded in the 
figure: number of positrons observed, thickness 
of aluminum used, intensity of the magnetic field, 
H, interval of initial energies, average value of 
the initial energies, E, and the average loss of 
energy in Mev, (AE),. 

It will be noted that the curve is roughly 
symmetrical and fairly sharp. A large number of 
tracks appear to have passed through the ab- 
sorber with no change in momentum, and a few 
appear to have acquired actual gains in mo- 
mentum. This can be attributed to several 
factors: (1) Coincident effects. (a) Occasionally 
a gamma-quantum may eject, from the wall of 
the cloud chamber, an electron which hits the 
absorber at a point opposite that at which a 
positron has entered. If this electron has an 
energy equal to, or greater than, the initial 
energy of the positron, it will seem that the 
positron has actually gained energy. (b) Another 
positron which is at first out of the light beam 
may be scattered into view on the far side of the 
absorber in such a way as to appear to be the 
emergence of a first particle. (2) Small angle 
scattering. When working with particles of such 
low energies as those used in this experiment, 
scattering is certain to occur. For every large 
angle scattering which can be easily observed, 
there is a considerable amount of small angle 
scattering which cannot be noticed. (3) Errors in 
measurement. 

Of these three factors, the last two are by far 
the most important. Only one case was observed 
in the 800 tracks obtained in which a positron 
appeared to gain so much energy that it was 
definitely a coincident effect with an electron 
emitted from the wall. Since it was required that 
the tracks appear to within 0.5 cm of the 
absorber on both sides, then from the geometry 
of the apparatus, it would be impossible for any 
track to have sufficient path length in the light 
beam to pass either over or under the absorber. 
It is still possible, of course, that a particle might 
be scattered in such a way as to appear to be the 
continuation of a track which has been stopped 
by the absorber. Small angle scattering is a factor 
which is always a great difficulty in this type of 


CURTIS 


experimental work. Its effect is to broaden the 
general curve. 

Undoubtedly the errors in measurement are 
the chief cause of the large number of tracks 
which appear to have lost no energy. Since the 
accuracy of measurement is not high enough to 
distinguish any loss of energy less than 0.05 Mev, 
a large number of tracks will be obtained for thin 
absorbers which seem to have lost no energy. 

The number of tracks which have definitely 
appeared to gain energy is not very large, 
constituting less than 5 percent of the total 
number. These cases are probably due to a 
combination of the two factors, small angle 
scattering and errors in measurement, although a 
few of the cases may be due to coincident effects. 

The most probable loss of energy may be found 
by inspection of the curve to be approximately 
0.09 Mev. The average energy loss may be found 
by two different procedures: (1) by determining 
the energy loss for each particle separately and 
then averaging; (2) by determining the average 
Ap and, from the value of A//j and the initial 
energy, finding the corresponding (AZ),. 

In order to test these two methods of con- 
sidering the data, 200 tracks were chosen having 
energies between 0.50 and 1.00 Mev. The first 
method gives an average energy loss of 0.128 
Mev. It is undoubtedly fortuitous, but the value 
determined by the second method was again 
0.128 Mev. 

All of the data were treated by both methods. 
For energies above 0.5 Mev, the two values were 

TABLE I]. Energy loss. (—) indicates a negative electron. 
An asterisk indicates that tracks have already been included 
in the data and that the initial energy of all the particles was 
the same. (?) This low value ts attributed to the fact that the 
initial energies of the particles are so low that large energy 


losses are discriminated against in measurement. 

















AVERAGE Loss 
(MEv) 
NUMBER INITIAL Most ‘cinta 
THICKNESS OF ENERGY PROBABLE 
(cM) TRACKS (MEv) Loss Exp | Tenen 
0.0275 | 123 |0.48to0.86| 0.09 |0.123 | 0.098 
77 ~—-|0.80to1.0| 0.10 |0.123 | 0.105 
96(—)|0.48 to 0.86, 0.08 | 0.127 | 0.098 
0.053 209 0.48 to 0.86) 0.19 | 0.266 | 0.205 
101 0.30 to 0.55} 0.21 | 0.230?! 0.214 
52* 10.60 10.270 | 0.204 
73* 0.74 10.278 | 0.204 
161 (0.67 to 0.93} 0.21 | 0.265 | 0.205 
114 0.90 to 1.40} 0.22 [0.285 | 0.206 
0.114 30— «(0.8 to 1.40) | 0.480 | 0.434 
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Fic. 3. Energy losses plotted against the number of 
positive electrons, for 0.053 cm aluminum and initial 
energies of 0.90 to 1.40 Mev. 


usually the same; when any difference occurred 
the value obtained from the loss of energy of the 
individual particles was considered to be the final 
result. Method (1) was found to be advantageous, 
for there is no necessity in this case for de- 
termining the average initial energy of the 
particles. 

The apparent gains in momentum are attrib- 
uted to instrumental error. They should there- 
fore be regarded as the accuracy of the measure- 
ment, and a corresponding error allowed for the 
values of the large energy losses. If the data are 
regarded in this way, then the average energy 
loss for the 0.0275 cm thickness of aluminum is 
0.123 Mev. (See Fig. 2.) 

The 0.053 cm aluminum was used in obtaining 
the greater part of the data because it proved to 
be sufficiently thick to prevent cases of apparent 
gain in energy without being thick enough to 
stop an excessive number of particles. 

The particles have been grouped into intervals 
of initial energy, which are kept as small as 
possible and yet have a sufficient number of 


tracks so that statistical variations are not too 
important. Some of the groups overlap slightly, 
but they have been kept as nearly as possible 
within 0.2 Mev intervals from 0.3 to 1.5 Mev. 
These are tabulated in Table I. 

The general shape of the curve (see Fig. 3) 
obtained from this thickness of aluminum illus- 
trates two points: first, that the curve is not 
quite so symmetrical as those for lesser thick- 
nesses but extends farther to the right; and 
second, that the width of the curve is in general 
greater. Theoretically, the half-width of the 
curve is dependent upon the thickness of the 
material traversed and the incident energy of the 
particle considered. Since this thickness is double 
that first used, the width of the curve will be 
more pronounced, due to the increased scattering. 

The unsymmetrical character of the curve may 
be attributed to several processes which cannot 
be distinguished from one another : (1) scattering, 
which tends to increase the straggling in the 
energy loss values since the particle actually 
traverses a greater thickness of material; (2) 
electron collisions, in which a large fraction of the 
initial energy has been lost; (3) radiative col- 
lisions. The last two of these have very little 
effect, for the probability of radiative collisions is 
very small at the energies considered in this 
work, and electron collisions in which a particle 
suffers a head-on impact are relatively rare, 
making the contribution from this process 
small. 

In Fig. 3 the distribution for 114 positrons is 
shown for particles whose initial energy ranges 
from 0.90 to 1.40 Mev. The most probable energy 
loss occurs at an energy of approximately 0.20 
Mev. The average energy loss is equal to 0.28 
Mev. This illustrates the nature of the curves 
obtained when sufficient numbers of particles are 
used so that statistical variations are not large. 
The other curves for this thickness are simi- 
lar in nature, the results of which appear in 
Table I. 

The amount of data obtained for the 0.114 cm 
aluminum is too meager to make any very 
accurate estimate of the energy loss, but from the 
results available, the most probable energy loss 
has been found to be approximately 0.35 Mev, 
and the average energy loss, 0.48 Mev. 
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Fic. 4. Energy losses plotted against the number of 
negative electrons, for 0.0275 cm aluminum and initial 
energies of 0.48 to 0.86 Mev. 


Negative electrons 


To compare the energy loss of positrons 
directly with the energy loss of electrons, radio- 
active phosphorus was used as a source of 
particles. The upper limit of the energy obtain- 
able is 1.5 Mev as compared with the limit of 
1.6 Mev for the cobalt (55) positrons, so that a 
direct comparison is possible without changing 
the experimental arrangement. There was a 
sufficient number of tracks obtained to show that 
there is no essential difference between the energy 
lost by the positive electrons and that lost by the 
negative electrons. The results procured from 96 
particles are shown in Fig. 4. The value of the 
energy loss is 0.121 Mev, in good agreement with 
the energy loss of the positrons. It can be 
definitely stated that, within the accuracy of 
these experiments, the energy lost by the 
positive electron is the same as that lost by the 
negative electron, at least for the energies of the 
particles considered here. 


CURTIS 


DISCUSSION 


The theoretical values of the average energy 
loss have been calculated from Bloch’s formula.® 
The energy used in this formula has been taken as 


center of the absorber. These values are com- 
pared with the experimental ones in Table I. 

Theoretically, there is no difference between 
the phenomena of energy loss of the positive and 
the negative electron, except for the fact that the 
positive electron may be annihilated somewhere 
along its path and thus decrease its effective 
range. It can be seen from the values given in 
Table I that the average energy loss determined 
by the experiment is about 20 percent higher 
than the theoretical value for the 0.0275 cm 
aluminum and that this difference increases to 
nearly 40 percent when the thickness of the 
absorber is doubled (0.053 cm aluminum). This 
fact is interesting, for it would seem to indicate 
that the discrepancy is increasing with the 
thickness, and that some process must be taking 
place which causes this. One explanation might 
be multiple scattering within the absorber, which 
would tend to make the actual path traversed by 
the particle greater than the actual thickness of 
the absorber. Whether this is the case or not 
cannot be determined until further data are 
obtained. 

It must be remembered that none of the tracks 
which are stopped by the absorber are included in 
the calculation of the average energy loss. These 
are excluded on the assumption that the particles 
have been scattered out of the light beam as weil 
as stopped by the absorber. If these tracks were 
included in the calculation, the value of the 
energy loss would be much greater. 

From the results of this investigation, it may 
be stated that the experimental value of the 
energy loss of the positive electron in passing 
through aluminum seems to be greater than the 
theoretical formula of Bloch would predict. The 
direct comparison with the negative electron 
shows that there is no difference between the 
amounts of energy lost by the two particles. This 
is in agreement with the conclusions of Anderson, 


®Heitler, The Quantum Theory of Radiation (Oxford 
Press, 1936), p. 218. 
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Blackett and J. G. Wilson, and Ruhlig and 
Crane, but is in disagreement with the results of 
Crussard and Leprince-Ringuet for cosmic-ray 
electrons. 

The author wishes to express his sincere 
gratitude to Professor J. M. Cork and Dr. R. L. 
Thornton for the privilege of working under their 
guidance with the University of Michigan 
cyclotron, and for many helpful suggestions 
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which have facilitated the progress of this work. 
He is deeply indebted to Dr. H. R. Crane for the 
original discussion which led to the undertaking 
of this particular problem and for his constant 
encouragement. The financial support which 
made possible the carrying out of this investi- 
gation came through the courtesy of the Horace 
H. Rackham Endowment Fund, which the 
author gratefully acknowledges. 
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Calculations on a New Neutron-Proton Interaction Potential 


R. G. SAcHs AND M. GorPpPpEeRT-MAYER 
Johns Hopkins University, Baltimore, Maryland 
(Received April 16, 1938) 


On the basis of a neutron-proton interaction potential of the form Ce~’ 


a 


r numerical calcu- 


lations for the deuteron and for the scattering of neutrons by protons have been made. 


INTRODUCTION 


HERE has recently been suggested,' on 

theoretical grounds, a new neutron-pro- 
ton interaction potential of the form J(r) 
= —(C/r) exp (—r/a), where a=h/Moc, Mo 
being the mass of the heavy electron. Of the two 
unknown constants, a determines the range, 
A=C-a the strength of the interaction. If this 
potential is to give the correct energy’ for the 
ground state of the deuteron (Ey)=2.17 Mev) a 
relation between A and @ must be fulfilled. For 
other two-constant potentials, namely the spher- 
ical well, inverse fifth power, exponential, and 
Gaussian type interaction, the dependence of the 
magnitude of the forces on the range of the forces 
and the neutron-proton scattering cross sections 
have previously been calculated.’ In this paper 
the same calculations are made for the new 
potential. 


1H. Yukawa, Proc. Phys-Math Soc. Japan 17, 48 (1935); 
19, 1084 (1937). N. Kemmer, Nature 141, 116 (1938). 
H. J. Bhabha, Nature 141, 117 (1938). 

*H. A. Bethe, Phys. Rev. 53, 313 (1938). 

3H. S. W. Massey and R. A. Buckingham, Proc. Roy. 
Soc. 163, 281 (1937). Morse, Fisk and Schiff, Phys. Rev. 
50, 748 (1936); 51, 706 (1937). 


1. THe GROUND STATE OF THE DEUTERON 


The differential equation for the eigenfunction, 
Yo, of the ground state of the deuteron is‘ 


h? a7uo ors 
— —+(—E))uj= —A Uo; Ug=TYo, 
M dr’ ria 


where Ey is taken positive. 
The transformation r=ax yields 


u’+al(e*/x— eo) u=0, (1) 


where a=Aad’M/h?, e =Ey/A; (2) 


u(x) =uo(r), and primes indicate differentiation 
with respect to x. The problem is now to de- 
termine, for given values of a, the constant €» 
which yields a wave function mu satisfying the 
boundary condition. 

In the range between x =0 and x=} (the value 
} is chosen for convenience) the integration of 
(1) is made by using the first six terms of a 
power series development of u(x). The boundary 
condition that u be zero (yo be finite) at x=0 is 


*H. A. Bethe and R. F. Bacher, Rev. Mod. Phys. 8, 82 
(1936). 
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TABLE I. Values of €0, a and A for various values of a. 
Values in parentheses are extrapolated. 

a 1.683 2.00 2.30 2.50 2.70 (2.78) 

€0 0 0.010 0.032 0.050 0.073 (0.108) 

a (X 10") 0.62 1.18 1.54 1.93 (2.39) 

A (Mev) 217 67.8 43.4 29.7 (20) 

then automatically satisfied. For the range 


beyond x=} Eq. (1) may be put in the form 
b'+p?+a(e-*/x—€o) =0, p=u'/u. (3) 


This equation was integrated numerically from 
large distances in to x=} by the method of 
Runge-Kutta.> The value of p for large x must 
approach — (deo)! in order that the wave function 
decrease exponentially at large distances. This 
was used as the boundary condition at x=7.5, 
where e~*/x is 10~* and can be neglected. Values 
of €9 were tried for a given value of a, until a 
value was found that would make the numerical 
solution join smoothly with the power series 
solution at x= 4. The values of ¢€9 for given values 
of a are listed in Table I. By using the known 
energy”? of the ground state of the deuteron, 
E y=2.17 Mev, the corresponding values of A 
and a have been calculated and are also listed 
in Table I. The values in parentheses are 
interpolated or extrapolated. 

The equation for the singlet state is: 
u’’+a(e-*/x)u=0. This is obtained from (1) by 
setting ¢)=0 since the binding energy for the 
singlet state is approximately zero. Eq. (3) is 
replaced by: 


where 


p'+p?+a(e*/x) =0. 


Different values of a were tried until one was 
found that made the two solutions join smoothly 
at x=}3. This vaiue of a gives a fixed relation 
between a and A for the singlet state. Its value is 
a= 1.683. 


2. SCATTERING 


In order to calculate the scattering for neutrons 
on protons, the wave equation is written for 
positive energy: 
u’’+a(e*/x) +6e)u=0; a=Aa?M/h?, 

e=E/A, (4) 


°C. Runge and H. Konig, Numerisches Rechnen (J. 
Springer, 1924) p. 286. 


GOEPPERT-MAYER 


TABLE II. Values of 59 calculated for the triplet and 
singlet states. 





En aX 108 ; R 
(Mev) (cm) Of 98 

me 0.62 2.382 1.418 
1.93 2.304 1.243 

(2.39) (2.277) 

(2.24) (2.286) 

5.4 2.39 1.877 

2.24 1.849 

13.7 2.39 1.520 

2.24 1.525 
rs 0.62 1.688 1.248 
1.93 1.487 0.862 

(2.39) (1.417) 

(2.24) (1.439) 





E being the energy of the neutrons in a coordinate 
system fixed with respect to the center of gravity. 
(If E, be the energy of the neutrons in a system 
fixed relative to the protons,t £,=2E.) The 
solution®: * of Eq. (4) for large values of x is 
asymptotically wu =c sin [ (ae)'x+6 ];and the total 
scattering cross section’ is c=47(h?/ME) sin? 6. 
The phases, 6, were calculated for different 
values of E,, by using the same type of series 
solution as before, out to x=3, and integrating 
numerically from there out to x=6.5 by the 
Runge-Kutta’ method for second-order differ- 
ential equations. The values of a and A were 
obtained from Table I for the triplet state; and, 
for the singlet state, A was calculated for those 
same values of a from a=1.683. Table II gives 
the results. The subscripts ‘‘t’’ and “‘s’’ refer to 
triplet and singlet scattering, respectively. 
Table III gives a comparison of the results 
with the results of Massey and Buckingham’ for 
the exponential and spherical well potentials. By 


TABLE III. Comparison of the present results for the phases 
6 with the results of Massey and Buckingham for the expo- 
nential and spherical well potentials. 





J(r) @X10%cm | En=2.2 5.4 13.7 21.5 Mev 
J=C r<a (3.107)| (2.277) (1.914) (1.461) (1.226) 
J=0 r<a 
(C/r)e"** 2.39 (2.277) 1.877 1.520 (1.417) 
Centrie 2.4 2.277. 1.930 =: 1.532 








6 Mott and Massey, The Theory of Atomic Collisions 
(Oxford, 1933) Ch. II. 


7 See reference 5, p. 311. 
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interpolation, the values of @ for the e~*/x and 
spherical well potentials were adjusted to give the 
same phase and therefore scattering cross section 
as the exponential interaction for 2.2 Mev 
neutrons. The phases for higher energies were 
then calculated for these values of a. Only the 
triplet phases were considered since the effect of 
the singlet state can be neglected* above £,,=2 
Mev. According to the table the values of the 
phases for the spherical well and exponential 
potentials diverge slightly but run almost parallel 
for different values of E,. However, the phases 
for the e~*/x interaction first go below and then 
cross over to go above the phases of the other 
two potentials. 

The observed® cross sections for the scattering 
of neutrons by protons are considerably smaller 
than those calculated with the usually accepted 
range of about 2 10-" cm. Ranges as small as 
0.5x10-" cm would give agreement, but the 


8]. Chadwick, Proc. Roy. Soc. A142, 1 (1933). W. H. 
Zinn, S. Seely and V. W. Cohen, Bull. Am. Phys. Soc. 13, 
14 (1938). R. Ladenburg and M. H. Kanner, Phys. Rev. 
51, 1022 (1937). E. T. Booth and C. Hurst, Proc. Roy. Soc. 
Al61, 248 (1937). 
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experiments are at present too uncertain to 
draw conclusicns. 


CONCLUSION 


It appears that there is little possibility of 
distinguishing between this and other potentials 
even by means of scattering experiments over a 
large range of neutron energies. With the cross 
sections for the exponential and e~*/x potentials 
adjusted to be approximately equal for 21.5 Mev 
neutrons, the total cross sections differ by only 6 
percent for 2.2 Mev neutrons. 

The ranges of force considered correspond to 
heavy electron masses ranging from about 
700m .(a=0.5 X10-" cm) to 150m,.(a=2.410-" 
cm). Since it is generally accepted‘ that the range 
of the interaction in any case cannot be much 
larger than 2X10~-" cm, it appears that a lower 
limit to the mass of the heavy electron would be 
about 100 electron masses, on the assumption 
that this potential is the correct one. 

The authors wish to thank Professor E. Teller 
for the suggestion of the problem and the interest 
taken in it. 
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The Scattering of Fast Electrons in Gases * 


J. B. Horner Kupert 
Palmer Physical Laboratory, Princeton University, Princeton, New Jersey 


(Received January 27, 1938) 


The angular distribution of elastically scattered electrons of energies between 49.5 and 


87.7 Kev has been measured in argon, neon, and helium over the range 0.3° to 6°. Comparison of 
the results with the theoretical curves calculated from the Hartree functions with the ‘“‘Born 


” 


approximation 


shows good agreement in the case of argon, and in neon except at the smallest 


angles, but there is no agreement at all in the case of helium. The discrepancies are very similar 
to those found at the smallest angles in studies of scattering of electrons of a few hundred volts 
energy, and it is felt that the explanation in terms of polarization of the atomic field by the 


passing electron may be accepted here also. 


INCE the pioneer work of Lenard! carried 
out more than forty years ago, various 
problems in connection with the passage of 


* The work here reported forms part of a dissertation 
presented to the faculty of Princeton University in candi- 
dacy for the degree of Doctor of Philosophy. The complete 
dissertation is on file in the Princeton University library. 
A preliminary report was given at the Washington meeting 


electrons through matter have attracted the 
attention of physicists down to the present time. 
The scattering of electrons (elastic and inelastic) 


of the American Physical Society, 1937; Kuper, Phys. Rev. 
51, 1024A (1937). 

+t Now at Washington Biophysical Institute, c/o Indus- 
trial Hygiene Division, National Institute of Health, 
Washington, D. C. 

1 Lenard, Wied. Ann. 51, 225 (1894); 56, 255 (1895). 
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Fic. 1. Sketch of scattering apparatus. 


in gases has been studied very thoroughly up to 
energies of the order of 2000 volts or so, but at 
higher energies (up to 100 Kev) interest has 
been concentrated instead on diffraction effects 
and on scattering in metal foils. While experi- 
mental difficulties (mostly attributable to the 
very small cross sections) prevented the com- 
pletion of the full program originally planned it 
was found possible to study the elastic scattering 
in argon, neon, and helium over a small angular 
range. Comparison of the experimental results 
with the theoretical curves calculated with the 
aid of the ‘Born approximation’ from the 
Hartree functions reveals agreements and dis- 
crepancies that are strikingly similar to those 
found at very small angles with electrons of a 
few hundred volts energy. Although various 
other possible explanations of the discrepancies 
must receive consideration, it seems likely that 
the effects are due to polarization of the atomic 
field by the passing electron? as in the low 
voltage experiments. 


? Mott and Massey, The Theory of Atomic Collisions, 
(Oxford, 1933). Chap. IX; Chap. X, Sec. 10. 
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KUPER 


The apparatus was similar in general design 
to that used by S. N. Van Voorhis* but built on 
a larger scale and modified where necessary to 
permit the application of higher voltages. In 
this type of apparatus an electron ‘‘gun’’ fur- 
nishing a well-collimated beam of electrons of 
the desired energy is mounted so as to pivot 
about an axis in the center of a region in which 
the scattering is to take place. The gas under 
investigation is admitted to this region and 
differential pumping is employed. A pair of slits 
defines the scattered beam selected, which then 
passes through an electrostatic energy analyzer* 
before reaching the collector. 

Electrons from the thin oxide-coated platinum 
ribbon filament F (see diagram, Fig. 1) are 
concentrated on the first slit, S;, by the ‘‘re- 
flector’’ electrode C, which is held about 25 
volts negative with respect to the filament. 
The latter was about 160 volts negative with 
respect to S;. The accelerating slit system, 5S; to 
S3, was designed to act as an electron lens, 
producing an approximately parallel beam. All 
three slits were 0.5 mm wide, the first two 
being 7 mm long and the third 5 mm. A ring 
shaped insulator of quartz Q served to maintain 
a separation of about 3 mm for the first gap, 
and S; was 25 mm below So. S3; and S, are 
mounted on the ends of a brass tube, 292 mm 
long, which is insulated from the rest of the 
apparatus. A compromise had to be made in 
the design of slits S; and S; in order to obtain 
sufficient pumping resistance and still have a 
reasonable beam which to work. The 
dimensions finally chosen were: length, 5 mm; 
width, 0.1 mm; and depth, 5 mm, with the 
distance between them 343 mm. A pressure 
ratio of about 1000 could be maintained between 
the scattering volume and the analyzer. The 
entire electron gun could be rotated about an 
axis perpendicular to the plane of the drawing 
through the center of V. Ss, located 71 mm 
below S;, was of the knife-edge type, 12.3 mm 
long and 0.2 mm wide. The deflecting plates for 
the energy analyzer A, and A: are arcs of circles 
of radii 20.4 and 20.9 cm, respectively, sub- 
tending an angle such that the electrons are 


with 


3S. N. Van Voorhis, Phys. Rev. 46, 480 (1934). 
* Hughes and Rojansky, Phys. Rev. 34, 284 (1929). 
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swung through 127° as required for refocusing.‘ 
The exit slit S; was 0.5 mm wide and 12 mm 
long. Beyond it was another fixed slit 1 mm 
wide, covered by an aluminium foil 0,01 mm 
thick which shielded the collector from stray 
fields. The collector was mounted on a reentrant 


‘ glass insulator G3; which was carefully shielded 
‘to prevent building up of surface charges. 


The electrical connections are indicated sche- 
matically in the diagram, Fig. 2. The high 
voltage was obtained from a Keleket trans- 
former-rectifier outfit provided with an extra 
filter consisting of two, 50 kv, 0.25 mf Pyranol 
condensers charged through 0.1 megohm re- 
sistors. The divider made of one hundred 0.1 
megohm, 50 watt, vitrified wire-wound radio 
resistors furnished the intermediate potential for 
the electron lens and potentials for the analyzer. 
Between 9 and 10 percent of the accelerating 
voltage was applied between S; and S». The 
system of condensers across the portion n—>p of 
the divider and from the midpoint of the high 
voltage apparatus to ground was necessary to 
bypass a small alternating current flowing from 
the main transformer through stray capacity to 
ground. Residual fluctuations in analyzer voltage 
limited the resolving power to 1 in 580, in 
contrast to a resolving power of 1 in 700 expected 
from the geometry. 

The electrometer tube used was a Western 
Electric D-96475, operated in the modified 
Barth circuit as described by Penick.’ Unusual 
care had to be taken with the shielding on 
account of the surges from the high voltage, and 
radiofrequency filters in the galvanometer and 
battery leads were indispensable. The sensitivity, 
which remained remarkably constant over a 
period of months, was 37,600 mm/volt. 

Leads to the pumps were connected at the 
various points marked P, a separate pump being 
used for the analyzer. The gas under investiga- 
tion was admitted to the apparatus at J after 
passing through an adjustable capillary leak and 
a liquid-air trap. The purity was checked 
spectroscopically, and on completion of the 
experiments the helium sample was analyzed 
with the mass spectrograph by Dr. L. G. Smith, 


® Penick, Rev. Sci. Inst. 6, 115 (1935). 
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who reported that no heavy impurity was 
present to as much as 1 part in 7000. 

Under the conditions of focus and slit align- 
ment finally chosen more than 90 percent of the 
electrons were concentrated into a beam less 
than 0.5° 
“wing’’ present on one side of the main beam.°® 


wide, but there was a rather large 


Occasional check tests were made and as no 
difference in the net scattered currents on the 
two sides of the beam could be found all readings 
were taken on the clean side of the beam. 
Although the currents obtained at scattering 
angles as small as 0.6° were of the order of 10~* 
times the direct beam or less the possibility of 
multiple scattering was carefully checked and 
excluded. Readings were taken with argon at 
various pressures, up to a limit considerably in 
excess of any ever used. Over a range greater 
than 5 to 1 the scattered current was very 
accurately proportional to the pressure. 

Before taking any readings the apparatus was 
always thoroughly outgassed by operation at a 
voltage well above that to be used. The back- 
ground current was measured over the range of 
angles to be studied then gas was admitted to 
the scattering chamber and readings were again 
taken at the same set of angles. In almost all 
cases the readings were taken by the rate of 
drift method (with 1 mm/sec. equal to 7.2 
x<10-'® amp.). After the background had _ been 
subtracted and allowance made where needed 
for differences in pressure or filament emission 
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CAPACITIES IN MICROFARADS 
RESISTANCES IN MEGOHMS 


Fic. 2. Schematic diagram of-electrical connections. 


®It was possible to obtain a beam with the ‘‘wings”’ 
symmetrically arranged, but as the scattered currents to 
be measured were exceedingly small it was thought best 
to use the distorted beam and have a smaller background 
current to allow for. 
































Fic. 3. Cross sections for scatter- 
ing in argon for electron energies 
from 49.5 to 87.7 Kev. The ordi- 
nates are as shown for the 49.5 Kev 
curve. The curves for higher volt- 
ages are successively displaced up- 
ward through 0.5 units. 


Fic. 4. Cross sections for scatter- 
ing in neon, for electron energies 
from 49.5 to 87.7 Kev. The ordi- 
nates for the 49.5 Kev curve are as 
shown. The curves for higher volt- 
ages are successively displaced up- 
ward through 0.7 units. 
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Fic. 5. Cross sections for scatter- 
ing in helium, for electron energies 
from 49.5 to 78 Kev. The ordinates 
for the 49.5 Kev curve are as shown. 
The curves for higher voltages are 
successively displaced upward 
through 1.4 units. 











the results of various runs were averaged. 
Different runs (made usually several days apart) 
always agreed with each other within the ex- 
pected errors of the individual points, about 
3 percent in the most favorable part of the 
angular range. A correction for the change in 
effective scattering volume with angle had to be 
applied. Although this type of apparatus is not 
suitable for determining absolute scattering 
cross sections, it is legitimate to assume that, 
provided the slits above and below the scattering 
chamber are parallel, the effective volume varies 
inversely as the sine of the angle of scattering. 
Accordingly the results were multiplied by a 
factor 1000sin 0, which gave numbers of a 
convenient size. Finally the points were fitted 
to the theoretical curves using an empirically 
determined multiplier. In the case of helium the 
experimental curves departed so widely from 
theory that it was necessary to use the average 
value of the multiplier found for argon and neon 
to locate the curves with respect to each other. 
Determination of this factor is practically equiv- 


alent to a measurement of the scattering volume, 
and confidence in the experimental results was 
increased by the fact that it remained nearly 
constant during the course of the work. Under 
standard conditions (pressure in scattering cham- 
ber 4X10-? mm and current at 0° equivalent to 
a drift of 210° mm/sec.) the average value of 
the multiplier was 21 to give numerical agree- 
ment with the cross sections given by Mott and 
Massey in units of 5.66 10-*° cm?. 

The results plotted on a logarithmic ordinate 
scale are shown in Figs. 3 to 5 where the solid 
lines are the theoretical curves. For clarity the 
curves have been separated from each other 
vertically and the order reversed. Actually, of 
course, the curves for the higher voltages should 
lie slightly below the low voltage ones. The 
curves for argon were successively displaced by 
0.5 units on the logarithmic ordinate scale, for 
neon 0.7, and for helium 1.4 units. The scale is 
correct in each case for the 49.5 Kev curve. 
The theoretical curves were obtained by graph- 
ical interpolation from the tables of cross sections 
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published by Mott and Massey.’ It is interesting 
to nqte that the measurements in argon at lower 
voltages, although extending over less than 6 
degrees, cover a range greater than that tabu- 
lated; the curves of Mott and Massey were 
extrapolated a small amount at each end. 

It is difficult to make a reliable estimate of 
the probable errors in this work, especially as 
they varied so much along a single curve. The 
uncertainties in the pressure measurements could 
hardly exceed 4 or 5 percent, and in any case 
would only affect the constancy of the empirical 
multiplier mentioned above. At the smallest 
angles the possible error of about 0.05° in the 
angle could account for an error of 16 percent, 
but this would of course diminish on going to 
larger angles. Because of unsteadiness in the 
high voltage and other factors affecting the grid 
current of the electrometer tube an uncertainty 
in the natural rate of drift of 0.3 mm/sec. was 
not impossible, and as the curves were often 
followed down to currents that small the points 
at the largest angles cannot be trusted very 
much. 

In general it may be said that as long as the 
theoretical curves are concave upward (as in 
argon or in neon except at the smallest angles) 
reasonably good agreement is found, but as soon 
as the theory predicts a bending over of the 
curve towards the axis (as it must in all cases 
if followed to small enough angles) the agreement 
disappears. The results for helium bear a strong 
resemblance to those for neon at small angles, 


7 Reference 2, p. 124, for argon and neon, p. 120, for 
helium. I have corrected four obvious errors in the second 
column of this table, and also expressed the cross sections 
in the same units. 


but unfortunately it was impossible to follow 
them further to see if the curves would merge. 
Numerous attempts to obtain curves for hydro- 
gen were all unsuccessful, but it was observed 
that the variation of cross section with angle 
was even more rapid than in helium, so the 
agreement with theory would have been poorer. 
The discrepancy found in this work seems to be 
exactly the same as that observed in helium at 
small angles with medium energy electrons.* It 
must be conceded that because of the limited 
resolving power of the energy analyzer some 
inelastically scattered electrons could have been 
recorded, and that they would tend to give rise 
to an excess of scattering at very small angles, 
but it is hard to understand why they would not 
be at least as abundant in argon as in the other 
gases, if not more so. It seems best to accept 
tentatively the qualitative explanation in terms 
of the polarization of the atomic field as appli- 
cable here also.® 

In conclusion the author wishes to express his 
great appreciation to his colleagues and the 
staff of the Palmer Physical Laboratory for 
many helpful discussions and for technical 
assistance. In particular it is a pleasure to 
thank Professor G. P. Harnwell, who originally 
suggested the problem, for his continued interest 
and support. 

5 Reference 2, p. 122. 

®An attempt was made to apply the calculations of 
Massey and Mohr (Proc. Roy. Soc. A146, 880 (1934)) to 
the curves for helium although their results are not given 
in a convenient form for this purpose. They predict a 
departure from the Born approximation curves at angles 
smaller than about 2.2° as observed, but the experimental 
increase of scattering at smaller angles appears to be much 
more rapid than expected; in the neighborhood of 0.6° the 


scattering is about six times the Born approximation value 
instead of twice. 
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The Saturation Requirements for Nuclear Forces 
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The inequalities between the coefficients M, H, B, V, representing the proportions of 
Majorana, Heisenberg, Bartlett and ordinary interactions in the symmetric Hamiltonian, 
arising through considerations of saturation and of the instability of odd-odd nuclei heavier 
than N" are represented graphically (Fig. 1). The allowed values of M, H, B, V correspond toa 
part of a plane bounded by three straight lines. These inequalities are sufficient for saturation at 
infinitely high density of nuclear particles and all conditions derivable from the necessity of 
saturation for the potential energy in the high density condition are derivable from those listed. 
The limits set on the interactions in the *P and 'P states by the inequalities are discussed. The 
results are summarized in Table I, Eqs. (7) and (8) and in the adjoining section. The sufficiency 
of the conditions is discussed in Section 2. The physical limitations are gone into in Section 3. 














T is customary at present to use in nuclear 
theory a symmetric Hamiltonian consisting 
of interactions between pairs of nuclear particles 
represented by potentials of the following types: 
Majorana exchange (exchange of space coordi- 
nates), Heisenberg exchange (exchange of space 
and spin coordinates), Bartlett exchange (ex- 
change of spin coordinates), ordinary non- 
exchange potential. It is, of course, not certain 
that the actual forces are of this type and that 
the saturation of the nuclear forces owes its 
origin to a preponderance of exchange forces 
rather than to saturation phenomena of other 
types.' It appears to be, nevertheless, of interest 
to systematize the available information? re- 
garding the possible proportions of the different 
types of exchange forces as obtainable from 
considerations of stability of heavy nuclei (satu- 
ration) and the instability of nuclei heavier than 
N* with an odd number of neutrons N and an 
odd number of protons Z. A brief elementary 
discussion is given below of the inequalities that 
must exist between the four kinds of potentials 
together with a diagrammatic representation of 
the values that are possible according to present 
evidence. A proof is then given by means of 
group character formulas that the elementary 
method contains all of the information derivable 
from requirements of saturation. 


1G. Gamow and E. Teller, Phys. Rev. 51, 289 (1937). 

2. Breit and E. Feenberg, Phys. Rev. 50, 850 (1936); 
D. Inglis and L. A. Young, Phys. Rev. 51, 525 (1937); 
E. Wigner, Phys. Rev. 51, 947 (1937); N. Kemmer, Nature 
140, 192 (1937); E. Feenberg, Phys. Rev. 52, 667 (1937); 
cf. also related discussion by H. Volz, Zeits. f. Physik 105, 
537 (1937). 





1. NECESSARY CONDITIONS 


The potential energy between a pair of nuclear 
particles will be represented by 


[MP™+HP#4BP84+V-1]J(r); J(0)<0, (1) 


where M, H, B, V are ordinary dimensionless 
numbers and P™”, P#, P® are, respectively, the 
Majorana, Heisenberg and Bartlett exchange 
operators. The function J(r) is taken to be 
negative in the vicinity of r=0. The coefficients 
M, H, B, V will be supposed to be connected by 


M+H+B+V=1. (2) 


The potential energy in the stable *S state of 
the deuteron is J(r). In the 4S state scattering 
experiments indicate a smaller degree of attrac- 
tion represented by the potential g/J(r) with 
g~1/2. Thus 


M-H-B+V=¢<1. (3) 


If a definite value of g, say g=3/5 is assumed, 
the above relations restrict the four quantities 
M, H, B, V so that only two of them are inde- 
pendent. The quantities M and H will be used 
as independent variables and a diagram using 
M and H as coordinates will be employed. The 
interactions in the *P and 'P states will be 
written as (?P)J(r), (\P)J(r). Positive values of 
(®P) and ('P) correspond to attractions in these 
states. These quantities are given by 


@P)=-—-M-—H+B+V=-2M-2H+1, 
(P)=—-M+H—-B+V=-—2M+2H+4. 
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Fic. 1. Graphical representation of inequalities. Origin of M, H coordinate 
system is at 0. Origin of B, V coordinate system is at JJ. Origin of (®P), ('P) 
coordinate system is at P. To any point in the plane there correspond the values 
of M, H, B, V, in these coordinate systems. The inequalities c,=0, H+2B—M=0, 
H—M=0 exclude all points above the lines c,=0, H+2B—M=0, H—M=0. 
The lines c,=0 are marked by c, in the figure. The region allowed by all the in- 
equalities considered here lies below the lines c3, c: and M=H. The shaded area 
I, II, III, IV is that part of the allowed region for which M, H, B, V are positive. 
The scale of (*P), (\P) differs from that of M, H, B, V and is marked in the 
coordinate system (°P), (‘P). Diagram is drawn for g=3/5. 


One can express also M, H, B, V in terms of 
(®P) and ('P) as independent variables: 


M=i711+q—(@P)—('P)], (3.3) 
H=%[1—q—(@P)+('P)], (3.4) 
B=({(1-—q+(P)—('P)], (3.5) 
V=3l1+¢+@P)+('P)]. (3.6) 


Denoting the projection of total spin along an 
arbitrary axis by of/2 and distinguishing be- 
tween protons and neutrons by suffixes 7 and v 
one has the following inequalities 


oq =4V—M—2H+2B=—-5M-—4H+3+q=0; 
o,=0,=0; Z=N (4.1) 


c2=2V—M-—2H+B=—3M—3H+}(3+q)=0; 
Z=0 (4.2) 


g¢,=¢,=0; 


¢c3=2V—-M—H+2B=—-3M—-3H+2 
=¢2.+}(1—g)=0; 
o,2=N; Gg=Z; 


Z=N (4.3) 


¢s= V-M-H+B=-—-2M-2H+1 
= +(—1+2c;) = 0; 
o,=N: Z=0 (4.4) 


¢;5=2V—-M—H+B= —3M—2H+}3(3+ q)=0; 
o,=N; o,=-—Z; Z=N (4.5) 


co=8V—-3M—4H+5B 
= —11M—9H+}(13+3q) =0; 
Z=N. (4.6) 


o,=N; o,=0; 


The first four inequalities are not new® and some 
additional conditions are known as well.’ It will 
be shown below that with the value of g indicated 
by present experimental data all the inequalities 
are consequences of (4.1), (4.3) and (5.2). If one 
assumes the inequalities (5), conditions (4.2), 
(4.4), (4.5) and (4.6) do not give any additional 
information. On the other hand, if the in- 
equalities (5) are not used then (4.5) gives new 
information inasmuch as it is sufficient for the 
validity of (4.1) and (4.3) in the region of 
small M and large /H. 
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The above six inequalities may be obtained by 
using trial wave functions contained in a volume 
that is sufficiently small to make J(r) equal to 
J(0). The trial wave functions give an expecta- 
tion value for the energy which is higher than 
that corresponding to an accurate solution of 
the wave equation. The numbers of neutrons 
and protons, N and Z, are contained in this 
expression as N? and Z? as well as in lower 
powers of N and Z. The sum of all terms in N? 
and Z? must, therefore, be positive or zero. 
Otherwise the energy in the most stable condition 
will go to — © as —N? or —Z*. The conditions 
(4.1) to (4.6) are obtained in this manner and 
by neglecting the Coulomb energy which is 
justifiable? in view of its relatively small value. 
In addition the nonexistence of stable odd-odd 
nuclei heavier than N' gives? 


H—-M+2B=—-M—H+1-q=0 (parallel 


spin of vy and 7 in unstable nucleus) (5.1) 
H—M=0 (antiparallel spin of vy and x 
in unstable nucleus) (5.2) 


The variables H, M are used in Fig. 1 as Car- 
tesian coordinates of a point in a plane. Any 
one of the inequalities c=0 has the geometrical 
significance of dividing the 7, M plane into an 
allowed part below the line c=0 and a forbidden 
part above that line. Only the part of the plane 
below the six lines given by the inequalities 
(4.1)---(4.6) and the two lines corresponding to 
(5.1) and (5.2) can represent a satisfactory set 
of M, H, B, V. The six lines corresponding to 
C,x=0 are marked as c, in the figure. The two 
lines corresponding to Eqs. (5.1) and (5.2) are 
labeled by these equations. It is seen that for 
qg=3/5 the permissible region of the plane is 
bounded by the three lines 


H=M. (6) 


The axis of H is also the negative of the axis of 
B with B=0 which comes at H=3(1—g). Simi- 
larly the values of V increase upwards along the 
M axis with V=0 coming at M=}3(1+gq). In 
the same figure a rotation through 7/4 gives 
lines of constant (*P) and (‘P). The values of 
(®P) increase as one moves up and to the left, 
those of (‘P) as one moves up and to the right. 
The scales for (?P) and (‘P) are shown on these 


c3=0, c,=0, 


AND E. 


WIGNER 


axes. Unity is seen to be represented by a smaller 
length in these coordinates than in the —J// 
system. As has been noticed by Feenberg,? the 
inequality c3;=0 determines an upper limit of 
—1/3 for (®P) which is of interest for the 
scattering of protons by protons. According to 
this the repulsion between protons in the *P state 
should be greater than that corresponding to a 
potential —J(r)/3. The condition is apparent in 
the diagram. The line c;=0 is one of the bounding 
lines of the allowed region and for it (@P) = —1/3. 
An upper limit is set for ('P) by the condition 
H=M. For this line ('P)=q and only smaller 
values of ('P) are allowed. If it is supposed that, 
in addition, the interaction is such*® that c,;=0 
then ('‘P) must be greater than that corre- 
sponding to the intersection of c;=0 with c;=0, 
which is —3g. The latter condition [('P) = —3q] 
was first given by Feenberg. If one similarly 
restricts the allowed region for (®P) by the 
requirement c,=0, then /J7<M imposes the 
additional condition (*P) = —%4(4q¢+3). Thus for 
6,:=0 
—1/32=(P)= — (4q+3)/9, (7) 
q=(*P)= —3¢. (8) 


Still closer limits on (*P) are set by requiring 
that V, M, H, B should all be positive. The 
corresponding region is shaded in Fig. 1 and is 
seen to be a trapezoid. Its vertices are the 
points J, JJ, IJJ, IV and the points at the 
intersection of c,=0 with cs;=0 and 7=M are 


Values of (3P), (P), V, M, H, B for points 
(1), U1), (ID, (IV), (V), (VD. 


TABLE I. 





























| " ‘ 
| INTERSECTING | | | 
Point | LINES (3P) ('P) V| Mid#H B 
(I) —_ —2%4—1| 3¢—6 | 3g—1| g+3 | 0 I-¢ 
H=0 5 5 | 0] 5 7 
(II) H=0 —q —1 0 I+q| 0 1—q 
V=0 2 | 2 
(IIT) V=0 —1 —q 0 1+q | l-@ 0 
B=0 2 2 
(IV) =0 —2—q | 3—6q| 2—q | 1+3¢| 1-9 0 
a=0 | a y ) ois i 3 
(V) a=0 1 | 3, | 1—3q | 14+3¢ | 1—3¢ 1+3q 
= ef pre | b- | es 
e3=0 3 | rie | 3 3 6 
| 
(VI) | a=0 —4q—3)| | 7g+3} q+3 | g+3 3—1lq 
| H=M | a | | is | 9 | 9 18 


3D. Inglis, Phys. Rev. 51, 531 (19.7); E. Feenberg, 


Phys. Rev. 51, 777 (1937). 
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designated by V and VJ respectively. In Table I 
are given the values of (°P), ('P), V, M, H, B 
for these six points. The values of V, M, H, B 
are very much restricted if one excludes negative 
values of these quantities. The simultaneous 
requirements of V, M, H, B=O and c,=0 
determine (*P) particularly well and place it 
between the values for the points J and JV. For 
q=3/5 these limits give —11/25>(®P) > —13/25 
and for g=1/2 they are —2/5>(8P)>-—1/2. 
If one removes the requirement c;=0 but keeps 
V, ---, B=O, then the points J and JJ set the 
limits on (@P) which give for g=3/5, —11/25 
> (@P) > —1. Similarly removing the requirement 
c,=0 but keeping V, ---, B=O the value of 
(‘P) is seen to lie between the limits set by 
points JJ and JV. Thus —1<('P) <}(3—69) 
which for g=3/5 becomes —1<('P) <—3/25. 
The requirement of positiveness of M, H7, B, V 
is seen to restrict (#P) within a narrower region 


than ('P). 


2. SUFFICIENT CONDITIONS 


The inequalities (4.1)---(4.6) have been ob- 
tained by considering states with special spin 
orientations that have been listed. It is not 
obvious that they correspond to the most 
stringent conditions on V, M, H, B because by 
a suitable change of the wave function it may 
be possible to lower the expectation value of 
the energy. It will be shown in this section that 
the requirement of the absence of terms in the 
expectation value of the potential energy, that 
go to — as —N?* and —Z’, in the approxi- 
mation of long range forces, gives conditions on 
V, M, H, B that follow from the inequalities 
(4.1), (4.3), (5.2). This means that these in- 
equalities are sufficient for saturation of the 
potential energy at high particle densities where 
by saturation is meant the impossibility of 
finding a trial wave function which gives, at 
high particle densities, a negative expectation 
value for the potential energy varying quad- 
ratically with the number of particles. For 
certain functions J(r) it can be shown then by 
an argument analogous to one previously used‘ 
that the same holds for every wave function 
independently of its extension in space. Since 


4E. Wigner, Proc. Nat. Acad. 22, 662 (1936). 
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the expectation value of the kinetic energy is 
positive these conditions are also sufficient for 
the saturation of the total energy. 

For wave functions with a_ small 
extension the J(r) in (1) may be replaced by 
J(0). The problem of finding a trial wave 
function with the lowest expectation value for 
(1) is then reduced to that of finding the eigen- 


spatial 


function of 
—MP™“—HP#—BP®—-YJ\-1 


which corresponds to the lowest eigenvalue. 
Since the four operators in this sum commute, 
there is a complete set of eigenfunctions common 
to all four of them. The lowest eigenvalue is 
obtained for that eigenfunction for which the 
sum of the four eigenvalues is lowest. 

It has been shown by Dirac® that the eigen- 
functions of P® correspond to definite values of 
the total spin S. Similarly, according to Feenberg 
and Phillips,® the eigenfunctions for P” have 
definite isotopic spin?’ JT and those for P™ are 
known’ to have definite “multiplicity.” The 
part of the potential energy which is proportional 
to the square of the total number of particles 
becomes 


E=—A*J(0){-3V +E +3(P+P?+P”) M 
+(3+7T?)H—({+S*)B}. (9) 


Here’? A=N+Z=atomic number, S=S/A, T 
=T/A, S=particle spin in units h, T=isotopic 
spin, P=greatest S possible for the multiplet, 
P’ =greatest T possible for this S, P’’ =greatest 
Y/A possible for these S, 7. 

Since 1/>0 (see diagram), it will suffice to 
consider the following cases: 


(1) M>0, H>0, B<0; 
(2) M>0, H>0, B>0: 
(3) M>0, H<0, B>O. 


The case M>0, H<0, B<0 need not be con- 
sidered because B+J7=}(1—q)>0, which is 
impossible for 17<0, B<0O. For each of the 


'Cf. P. A. M. Dirac, The Principles of Quantum Me- 
chanics, second edition (Oxford, 1935) §61. 

6 E. Feenberg and M. Phillips, Phys. Rev. 51, 597 (1937), 

7E. Wigner, Phys. Rev. 51, 947 (1937). 

SE. Wigner, Phys. Rev. 51, 106 (1937). The notation 
used in the present paper differs from that of this reference 
in referring to S, 7, Y of Eqs. (11) of the latter as PA, 
Paek ae 
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above cases the expression for the energy will be 
minimized by choosing suitably a7 PP. . 
It will then be required that E, as given by 
Eq. (9), be >0. For any other choice of s. f. 
P, P’, P” it follows that E>0. In this way one 
obtains, therefore, sufficient conditions for E>0 
and it is necessary for nuclear stability that all 
such conditions be fulfilled. Since the actual 
energy may be lower than (9) these conditions 
may not be sufficient for stability but it appears 
likely that they are,‘ since this has been shown 
to be the case for a special class of functions J(r). 

For case (1) the minimum £ is obtained for 
P=P’=P"=T=S=0 which leads on substitu- 
tion into Eq. (9) to c,<0 and is therefore con- 
tained in the preceding inequalities. For case (2) 
one may consider: (a) M—2B>0and (b) M—2B 
<0. For (a) the combination 


1P?M —5°B = P(}M—B)+(P?—8*)B 


is a minimum when P=S=0. The minimum of 
E corresponds, therefore, to P=S=P’=P"”=T 
=0 which gives c; <0 as in case (1). For (b) the 
combination }P?é—S?B is a minimum when P 
has its maximum and P?—§? has its minimum 
value. Since P is the largest possible S, the 
minimum P?—§?=0 and the maximum P=1/2 
which corresponds to a parallel orientation of 
all the particle spins. The best values are thus 
P=S=1/2. For these values the minimum E 
corresponds to P’ = P” = T =0. Substituting into 
Eq. (9) one finds then c;<0 in agreement with 
the fact that in the inequality (4.3) all particle 
spins were taken to be parallel. 

For case (3) it is convenient to rearrange the 


expression for E: 


E=—A*J(0)(§M+3H—§(q+3)J+E£., (9.1) 
E,=—A*J(0)[$}M(P?+P2+P") | 
+H(T?+8*) —3(1—g)S?]. (9.2) 


The variable terms are contained in £,. Here 
H<0, —3(1-—q) <0. Hence, for fixed P, P’, P”, 
a minimum is obtained for T=P’, S=P’”. 
For these values Eq. (9.2) becomes 


E,=—A*J(0)|[H+4M—}(1-9) 1" 


+(H+3M)P2+4MP’}. (9.3) 


In this expression the coefficient of P”’ is positive 
since M>0. The coefficient of P? is necessarily 


AND 
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less than that of P”. Therefore, there are three 
possibilities: (a) H+3M—3(1-—q)>0, H+3M 
>0; (b) H+3M—3(1-q) <0, H+3M<0; (c) 
H+4M—3(1-—q)<0, H+3M<0O. For (a) the 
minimum £ is obtained if P=P’=P’”=0. In 
this case, therefore, S=7=0 as well. The 
minimum E£ is the same as in case (1) and one 
obtains c; <0. For (b) the minimum E£ is obtained 
by using maximum P=} and minimum P’ =P” 
=0 as in (2b). The condition c;<0 is, accord- 
ingly, obtained again. 

For (3b) the coefficients of P? and P” are both 
negative. In order to obtain a low E£ it is now 
advantageous to use large values of P and P’ 
and at the same time P”’ must be made as small 
as possible. It is not possible, however, to vary 
P, P’, P” independently of each other as the 
minimum of £E is approached. In fact the 
minimum for Eq. (9.3) with independent varia- 
tions of P, P’, P’”’ corresponds to P=0, P’=1/2, 
P"”"=0 which is an impossible set of these 
quantities since P=P’=P”. The way in which 
P, P’, P” may vary can be inferred*® from their 
expression in terms of Ag=A;=A,=A,=0. 


1 
P=—(A4—A3—A2—A)); 
2A 
1 
P’ =——(Ag—A3+A2—A}) 
2A 
1 


P" =—(Aq—A3—Ao+A)); 
2A 


A =A4yt+A3+Aot+A. 


If here A; >0O one can decrease A; and increase 
A, by an equal amount. This leaves A, P” 
unchanged and increases both P and P’ which 
gives a lower energy. Thus A,;=0 will include 
the most favorable cases. With this restriction 


P"” =P+P’—3 and 


E,= —A?J(0)} (H+M—}(1—q) ]P? 


M 
+(M+H)P"+M(PP'+1)-(P+P)}. (9.4) 


Further 


Ay/A =P+P’>A;3/A=}-—P’>A:2/A 


=}—P>0. (9.5) 
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The above inequalities (9.5) confine the possible 
P, P’ to a triangular region in the P, P’ plane 
bounded by the lines P=}, P=P’, 2P’+P=}. 
In this triangle the values of P, P’ can be varied 
arbitrarily. If one uses the facts that in the 
present case /7<0, +2H<0 while according 
to conditions already derived M+J/>0 one 
finds by means of Eq. (9.4), after a somewhat 
lengthy but quite elementary discussion, that 
the minimum £ corresponds to P=1/2, P’=0, 
P”=0 which gives c;<0 as in (3b). This ex- 
hausts the possibilities. 


3. FINAL REMARKS 


It should be observed that the conditions (4) 
have been derived on the assumption of the 
absence of other nuclear forces than those given 
in Eq. (1). In particular forces of the ‘“‘many 
body” type could bring about saturation in 
spite of the inequalities (4) not being fulfilled. 
Even if under ordinary circumstances these forces 
should be negligible they may possibly be of 
special importance in the state of high density 
considered for saturation. A similar remark 
applies? to the influence of a velocity dependence 
of the potential. 

It should also be pointed out that the in- 
equalities (4) are obtained by a consideration of 
an idealized state of high density which would 
be the actual state, if the inequalities (4) were 
not satisfied. On the other hand, inequalities (5) 
follow from the consideration of actual nuclei 
and involve several approximations. One can 
derive them on the basis of the high density 
model also but this certainly is not the condition 
in actual nuclei. All other considerations on this 
question are based, apart from some, rather 
rough approximations, on the assumption that 
H and B are small compared’? with M and 


®J. A. Wheeler [Phys. Rev. 52, 1083 (1937), especially 
p. 1106], has investigated the saturation properties of 
many particle forces. 

10 The following qualitative arguments suggest that H 
and B separately are small compared with M and hence 
neither of them has a large negative value. (1) The 
empirical mass differences N¥%—C¥8+C'"—N"“=0.0026, 
Cu — Bu+ Be!®— B!°=0.0027, He*— H*+ He*— Li* = 0.0039 
give approximate values of —2(H+B).J(r) in N™, B®, Li’. 
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should not be used, therefore, to establish the 
inequality J7<./. It is rather the apparent 
agreement of these considerations with experi- 
ment that gives one confidence in the inequality 
IT<M. The underlying experimental rule" holds 
without exception for’ A>14. In fact, this 
seems to indicate that // is materially smaller 
than ./, as otherwise accidental factors would 
be likely to make at least some of the heavier 
odd-odd nuclei stable. 

On the other hand, information regarding the 
approximate symmetry of the nuclear Hamil- 
tonian arises out of experimental material with 
ordinary rather than high nuclear densities. It 
is more closely related to the inequalities (5) 
than to (4) and one may be skeptical about the 
latter from this point of view. 

Very little is known about the space de- 
pendence of the function J(r) and still less about 
J(0). It is possible that the coefficients M, //, 
B, V should be multiplied by different functions 
J(r) and there is a chance that the four J(r) 
are roughly the same for distances of the order 
2X10-" cm but that the J/(0) are widely differ- 
ent. If such should be the case the saturation 
inequalities (4) will have very little meaning. 

One of us (E.W.), would like to express his 
indebtedness to the Wisconsin Alumni Research 
Foundation for its support of this work. 


The last value is not very reliable experimentally and can 
be expected to be great on account of the loose packing of 
He®. The first two are seen to be practically the energy 
difference 'S—*S in H*®. If H and B were separately large 
such close agreement would not be expected for then the 
wave functions would be changed appreciably from the 
symmetry required by the preponderance of M. (2) Con- 
siderations of the kind as the one leading to Table II in 
footnote 7 and the apparent existence of S* and Ca*® [A. O. 
Nier, Phys. Rev. 53, 282 (1938)] may be considered as 
qualitative arguments supporting the above view. From 
this point of view (5.1) is a more restrictive assumption 
than (5.2) inasmuch as in the region of positive H and B 
a greater part of the plane is eliminated by (5.1). Small 
negative values of H or B cannot be excluded by the above 
arguments. Nevertheless, the set of values M=14/12, 
H=-—7/12, B=10/12, V=—5/12 proposed by Volz, 
reference 2, contradicts the inequality c; <0. 

1 W. D. Harkins, Phys. Rev. 19, 135 (1922). J. Am. 
Chem. Soc. 39, 856 (1917); 45, 1426 (1923); J. Frank. Inst. 
194, 165, 329, 521, 645, 783 (1922); 195, 67, 553 (1923); 
Phil. Mag. 42, 305 (1921). 

122 For A=14, the conditions are different and odd-odd 
nuclei can be stable even if H is quite small. 
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A guard ring method for the measurement of the thermal 
conductivities of liquid metals in which the furnace auto- 
matically attains equilibrium is described. The design of 
the furnace permits frequent and easy replacement of 
specimens; a layer of paraffin over the liquid metal pre- 
vents oxidation by the air. The thermal conductivities of 
mercury, of sodium, and of sodium amalgams ranging in 
concentration from 70 to 94 atoms percent sodium are 
measured over these temperature ranges: mercury, 40°C- 


220°C; sodium 85°C-210°C; sodium amalgams, 100°C- 
150°C. Sodium behaves as a normal metal, but mercury 
and the sodium amalgams are abnormal. Some evidence 
for the existence of the intermetallic compounds Na;Hg, 
and Na;Hg in the liquid amalgams at 100°C is presented 
by breaks in the thermal conductivit y-concentration curve. 
An eutectic at 86.5 atoms percent sodium is evident from 
the same curve. Calculations of Lorentz’ number for 
sodium and for mercury are given. 





ERCURY, sodium, and sodium amalgams 

are of particular interest when in the liquid 
state. The behavior of the electrical conductivity 
of mercury with temperature causes mercury to 
be classed with the abnormal metals; the same 
thing may be said of its thermal conductivity, 
though there is too much disagreement in the 
existing measurements to be certain of this fact. 
Sodium and mercury are known to form many 
intermetallic compounds in the solid state which 
may yet exist in the liquid state. The thermal 
conductivity of liquid sodium has never been 
investigated over any range of temperature. 
Hence, this investigation was carried out. 


APPARATUS ARRANGEMENT 


A method based on the use of a tall column of 
liquid for measuring the thermal conductivities 
of the liquid metals was developed. The major 
details of the furnace, including the significant 
dimensions, are shown by Fig. 1. The liquid 
metal was contained in a cylinder of sheet 
asbestos which had been formed about a bar of 
the proper diameter, and which had been coated 
thinly on both sides with Insalute cement. This 
cylinder rested on a bar of iron of the same 
diameter in the furnace. Two thermocouple tubes 
of Pyrex glass about 1.2 mm in diameter were set 
in the cylinder at right angles to the axis, with 
the aid of the cement. The tubes, which were 
approximately two centimeters apart, were set in 


*Now at Naval Research Laboratory, Washington, 
~~ <. 


from the ends of the liquid metal column by 
about one centimeter in order to insure a uniform 
flow of heat throughout the section of liquid 
metal whose temperature gradient was measured. 

The heater which supplied the heat that was 
conducted down the column was wound uni- 
formly of 38-gauge chromel wire on a sheet of 
mica. This was sealed in by the Insalute cement 
between two circular plates of stainless steel that 
were bolted together. After completion, the 
heater, which was only 5 mm thick, was allowed 
to float on the liquid metal. A layer of paraffin 
three or four millimeters in thickness above the 
heater protected the metal from oxidation or 
evaporation. 

The details of the guard-ring arrangement are 
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Fic. 1. The furnace. 
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of interest. About the asbestos cylinder were four, 
concentric, sheet metal cylinders which had the 
upper ends closed. Each one was held in position 
through bolts at the lower ends of the cylinders 
to a baseboard of pressed asbestos. The cylinder 
C was the guard ring proper, for it supported a 
uniform winding of nichrome wire around the 
wall and across the upper end which was used to 
maintain the average temperature of the cylinder 
C equal to that of the liquid metal column. A 
layer of sheet mica extended over each of the 
cylinders A and B; on each layer was wound 
nearly 25 feet of 30-gauge pure platinum wire 
which was uniformly spaced on the wall and 
across the top of the cylinder. Leads to the two 
resistance thermometers thus fashioned passed 
out of the top of the furnace through glass tubes. 

Of course, the two sheet metal cylinders A and 
B reduced the lateral radiation from the liquid 
metal. But their main function was to support 
the two resistance thermometers which were used 
automatically to prevent heat transfer to or from 
the liquid metal. This was done simply by 
keeping the ratio of their resistances equal to the 
value which existed at room temperature; i.e., 
to the value at which the heat transfer was zero 
with all of the furnace at the same temperature. 
As long as the ratio had this value the two 
cylinders A and B were at equal temperatures, 
which meant that there was no heat transfer 
from the column of liquid metal. 

It is true that the heat transfer by radiation or 
conduction between two surfaces whose average 
temperatures are equal will be zero only if 
Newton’s law for the transfer of heat holds; i.e., 
IZ=k(7T,—T:2). But at the operating tempera- 
tures of the furnace, which varied between 40°C 
and 220°C approximately, Newton's law would 
hold quite well for the small temperature differ- 
ences which could have existed. 

To facilitate changing one liquid metal for 
another, the construction of the furnace per- 
mitted the upper part which held the four 
cylinders to be raised away from the asbestos 
cylinder. The lower part of the furnace consisted 
essentially of the iron bar which held the 
asbestos cylinder and which served to adjust the 
mean temperature of the liquid metal to any 
desired value. 

A voltage regulator was employed to regulate 


the voltage at the heater to within 0.3 percent 
while the furnace was in use. This was necessary 
because of the time required for the furnace to 
attain equilibrium. 

The automatic apparatus for bringing the 
furnace to equilibrium functioned as follows : The 
short platinum leads of the two resistance 
thermometers were connected at a_ terminal 
board on the upper furnace to manganin resist- 
ances which made up the other two arms of a 
Wheatstone bridge. Connected across the bridge 
was a sensitive galvanometer. The mirror of this 
galvanometer reflected the light rays from a lamp 
filament to a lens which formed a bright image of 
the filament on a box housing a photoelectric cell. 
When the bridge circuit was balanced, the image 
of the lamp filament was just at the edge of an 
opening in the box before the photoelectric cell. 
But, when the bridge was unbalanced due to the 
passage of heat outward from the liquid metal, 
the galvanometer mirror reflected light into the 
photoelectric cell. After amplification, the re- 
sultant photoelectric current caused a relay to 
close, thus shunting part of the resistance in the 
power line to the furnace winding on the cylinder 
C. This increased the current to the furnace 
winding by about five or ten percent. As soon as 
the bridge circuit approached balance the 
galvanometer mirror started to rotate back to its 
first position, and by the time a very slight 
unbalance had been reached in the reverse 
direction caused by the passage of heat inward to 
the liquid metal, the galvanometer mirror had 
turned far enough so that light no longer passed 
through the opening to the photoelectric cell. 
Then the amplifier and relays operated in a 
reverse manner and removed the shunt from the 
previously mentioned resistance. A calculation 
indicated that a temperature difference of about 
one-hundredth of a degree centigrade between 
cylinder A and cylinder B was sufficient to cause 
the relays to operate. 

A standard resistance of one ohm was included 
in one of the current leads to the heater for the 
measurement of the current. Leads to a high 
resistance voltage divider were connected directly 
to the heater terminals within the furnace 
eliminating the necessity for a correction to the 
voltage measurements caused by the resistance 
of the current leads. 
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The chromel-copel thermocouples which were 
used were calibrated at the boiling points! of 
water, napthalene and benzophenone. Six couples 
that were calibrated could be represented by the 
following equation to within 0.1 percent: 


E=4.463 X 10-*t+5.017 X 10~*#? — 3.556 X 10-4 


millivolt, where the temperature is in degrees 
centigrade. 

For the purpose of getting the distance between 
the thermocouples as accurately as possible ‘a 
traveling microscope was bolted to a frame that 
was in turn bolted to the lower part of the 
furnace. The microscope was directed into the 
asbestos cylinder in line with the axis and just to 
one side of the two parallel thermocouple tubes. 
A very small mirror that had been attached a 
short distance from the objective end of the 
microscope turned the light rays from a thermo- 
couple tube so that the microscope could bring 
them in focus. The microscope, without the 
mirror, was used to get the cross-sectional area of 
the asbestos cylinder simply by measuring the 
rise in the mercury level in the cylinder which 
resulted when a weighed quantity of pure 
mercury was poured in. 

It is estimated that the maximum error in the 
absolute value of the thermal conductivity 
measurements which resulted from the neglected 
thermal expansion of the asbestos cylinder 
amounted to 0.15 percent. 


THERMAL CONDUCTIVITY OF MERCURY 


Many measurements of the thermal con- 
ductivity of mercury were made for the purpose 
of checking the operation of the furnace as well 
as to determine the true manner of variation of 
TABLE I. Effect of temperature gradient on the thermal conduc- 


tivity of mercury (cross section of liquid metal = 18.70 cm?). 




















| CONDUCTIVITY 
| (Watt/cm®C) 
MEAN Temp. | TEMP. GRADIENT | POWER MEAS- CALCU- 
(°C) (°C/cM) |((WattTs) URED LATED 
74.7 4.90 | 9.82| 0.0999 | 0.0999 
78.9 9.52 17.93 | 0.1007 | 0.1008 
72.6 4.96 9.20 | 0.0991 | 0.0991 
66.9 9.82 17.87 | 0.0973 | 0.0978 
| 





1 Nat. Bur. Stand. Technologic Paper No. 170. 
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the thermal conductivity with temperature. The 
mercury used here and in the sodium amalgams 
came from the Mallinckrodt Chemical Company 
and was certified pure. 

It is to be expected that if there were more 
than a negligible heat transfer laterally from the 
liquid metal it could be detected in a change in 
thermal conductivity which would result from a 
change in the temperature gradient down the 
liquid metal at a given mean _ temperature. 
Several of the measurements in which the 
temperature gradient varied by about a factor of 
two are given in Table I. The calculated value 
was arrived at by using the observed tempera- 
ture coefficient from the data shown in Fig. 2. 
From the close agreement exhibited in this table 
and the other measurements it was decided that 
there was no heat transfer through the guard 
ring. 

By taking the thermal conductivity of asbestos 
paper as 0.0006 cal./em cube —°C, a calculation 
gave the ratio of the heat conducted by the 
asbestos cylinder to that conducted by the 
mercury at 0°C as 0.2 percent. This was the 
maximum value; hence it was neglected in all 
calculations. 

The results of 51 
conductivity fit the equation 


observations of thermal 


K,=0.0819(1+0.00314¢ —0.000002957*) 


in watt/em —°C, with an average deviation, 
regardless of sign, of 0.6 percent from the curve. 
These results are in agreement with the limited 
results of Nettleton? and of Berget,* but disagree 
with the results of Gehlhoff and Neumeier,* 
except as to the fact that mercury has a positive 
temperature coefficient of thermal conductivity. 


THERMAL CONDUCTIVITY OF SODIUM AMALGAMS 


The sodium amalgams which lie between 70 
and 94 atomic percent sodium were investigated. 
The sodium, which had been stored in benzene, 
was first freed from the scale that enclosed it by 
melting it under paraffin; then it was amal- 
gamated with mercury under paraffin. A special 


2H. R. Nettleton, Proc. Phys. Soc. 26, 28 (1913); 27, 
129 (1915). 

3 Berget, 1900; Glazebrook, 
Physics, Vol. 1 (1922). 

4 Gehlhoff and Neumeier, Verh. d. 
(1919); 15, 876, 1069 (1913). 


of Applied 


Dictionary 


phys. Ges. 21, 201 
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THERMAL 


glass pipette was used to transport the various 
amalgams to the furnace so that the amalgams 
were never in contact with the air. 

The ordinary method of analysis by titration 
of samples of the amalgams that had been added 
to beakers of distilled water was resorted to, 
using standardized hydrochloric acid with methyl 
orange as an indicator. 

A condensed summary of the data over the 
amalgams is given in Table II. 

A plot of the results for the first 12 alloys has 
been made in Fig. 3 to show that the points for a 
given alloy fall on a straight line over the 
100°C to 150°C. 
From these lines the thermal conductivities, both 
at 100°C 150°C The 
results are given by Table III. 

The behavior of the thermal conductivity of 
the sodium amalgams with concentration changes 
at 100°C is indeed interesting, as Fig. 4 shows. 
For sodium amalgams which contain less than 80 
atomic percent sodium by up to a considerable 
amount, the thermal conductivity lies below that 
of either component metal alone. Also, there are 


temperature interval from 


and at were obtained. 


certain breaks in the curve which, if real, must 
throw some light upon the chemical structure of 
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Fic. 2. The thermal conductivity of mercury. 


the amalgams. Two of these breaks, existing at 
71.4 and 75.0 atomic percent sodium, respec- 
tively, may be caused by the presence of the 
intermetallic compounds Na;Hge and Na;Hg. 


TABLE II. Thermal conductivity of sodium amalgams (thermocouple separation/cross section of liquid metal=0.1103 cm™~'). 


ee DIFFER- 


THERMAL AToMK 
TEMP. ENCE POWER , CONDUCTIVITY | PERCENT 
ALLoy| (°C) | (°C) (Watts)| (Watt/cm —°C) SODIUM 
1 147.6 20.55 19.02 0.1021 70.0 
127.6 19.72 17.29 0.0967 
97.4 16.85 13.42 0.0879 
2 148.6 20.21 18.84 0.1027 71.4 
125.8 19.21 16.89 0.0972 
100.2 17.08 14.09 0.0910 
127.1 19.61 17.28 0.0972 
150.5 20.23 18.98 0.1034 
3 151.2 20.36 19.25 0.1042 73.1 
127.2 19.03 16.85 0.0977 
102.7 17.04 14.07 0.0911 
4 149.4 19.44 18.56 0.1052 73.8 
127.7 19.36 17.39 0.0990 
104.0 17.17 14.47 0.0929 
5 154.1 20.35 19.82 0.1074 75.1 
102.6 16.96 14.26 0.0928 
126.7 18.74 16.98 0.0998 
6 153.7 19.98 19.72 0.1088 76.1 
127.2 18.74 17.16 0.1010 
101.6 16.50 14.07 0.0940 
7 153.2 19.76 19.92 0.1112 77.6 
131.8 17.60 16.92 0.1060 
105.5 17.06 15.08 0.0975 
106.0 17.03 15.11 0.0978 
8 152.0 19.24 19.65 0.1126 79.1 
127.3 17.84 17.11 0.1058 
104.3 16.35 14.76 0.0995 




















| MEAN | Dirrer- | THERMAL ATOMK 
| Temp. ENCE | Power | CONDUCTIVITY PERCENT 
| ALLoY eC) | (°C) (Watrts)| (Watt/cm —°C) SopIUM 
9 153.2 | 19.06 19.88 | 0.1150 80.0 
| 126.3 | 17.82 17.35 | 0.1073 
101.2 15.82 14.34 0.0999 
10 156.9 18.91 20.34 0.1186 80.4 
125.9 17.43 17.39 0.1101 
103.7 15.79 14.75 0.1030 | 
| | 
il 151.7 17.82 19.71 | 0.1219 82.1 
125.8 16.41 17.17 0.1154 
102.2 15.01 14.65 0.1078 
102.3 15.02 14.64 0.1077 
136.1 18.82 20.17 0.1181 
| 
12 155.8 17.09 20.50 0.1323 83.9 
128.5 16.16 18.09 0.1235 
103.2 14.43 15.15 0.1158 | 
13. | 150.1 15.25 19.88 | 0.1437 | 86.4 
124.0 14.05 17.18 0.1348 | 
14 151.6 14.73 20.52 0.1537 | 86.9 
131.8 14.51 19.39 0.1474 
15 160.5 15.04 23.20 0.1701 88.4 
125.6 13.11 18.6¢ 0.1569 
107.2 12.41 16.77 0.1490 
16 | 160.3 14.18 24.04 | 0.1870 90.6 
131.4 13.23 21.00 0.1750 
17 155.5 11.32 25.54 0.2486 93.7 
126.6 10.95 23.01 0.2317 
100.2 13.92 26.93 0.2134 
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Fic. 3. The thermal conductivities of the sodium amalgams. 
Ihe numbered curves are for the alloys listed in Table II. 


Johnson and Fernelius® record the presence of 
these compounds in the solid sodium amalgams 
within the concentration range of Fig. 4. Several 
experimenters® have sought to detect the exist- 
ence of these compounds, and others, in the 
liquid state but to no avail. The upper break in 
the curve occurs at 86.5 atomic percent sodium, 
which concentration agrees well with that of the 
eutectic which Vanstone’ found at 85.5 atomic 
percent sodium in the thermal diagram. 

Perhaps the most reasonable interpretation of 
the drops which may be found in the thermal 
conductivity concentration curve at 100°C at 
71.4 and at 75.0 atomic percent sodium when the 
sodium content is increased, results by assuming 
that they are caused by decreases in the free 
electron concentrations when the compounds are 
formed. At the eutectic there is an increase in the 
number of free electrons. 

With regard to the possibility of the liquid 
sodium amalgams being colloidal in nature as 
suggested by Paranjpe and Joshi,° it is evident 
that these data present nothing in favor of sucha 
view. The breaks in Fig. 4 occur either at 
compounds or at an eutectic, and at no other 
points. 


THERMAL CONDUCTIVITY OF SODIUM 


The same procedure that was employed hereto- 
fore was used in the measurements on sodium. It 
was necessary, of course, to increase the buoyancy 


5 W. C. Johnson and W. C. Fernelius, J. Chem. Ed. 7, 
981-99 (1930). 

6 Vanstone, J. Chem. Soc. London 105, 2617 (1914); 
Bornemann and Miiller, Métallurgie 7, 396 (1909). 

7 Vanstone, Trans. Faraday Soc. 6, 42 (1911). 

8G. R. Paranjpe and R. M. Joshi, J. Phys. Chem. 36, 
2474 (1932). 
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Fic. 4. The thermal conductivities of sodium amalgams. 


of the heater on the molten sodium by means of a 
small support. Because of the high thermal 
conductivity of sodium the temperature gradient 
was quite small and varied from 1.15°C/cm to 
3.45°C/cem at the extremes. Fig. 5 gives the 
results which were obtained. 

The thermal conductivity of sodium in the 
range from 100°C to 210°C may be represented 


by the equation 


K,=0.860[ 1 —0.338 & 10-*(¢— 100) 
—2.6X10-*(t— 100)? ] 


TABLE III. Thermal conductivities of sodium amalgams at 
100°C and 150°C. 














THERMAL CONDUCTIVITY 
(Watt/cmM —°C) 
ATOMIC PERCENT E 
ALLOY SODIUM 100°C 150°¢ 
1 70.0 0.0887 0.1029 
2 71.4 0.0907 0.1031 
3 73.1 0.0903 0.1039 
4 73.8 0.0916 0.1052 
5 75.1 0.0921 0.1063 
6 76.1 0.0934 0.1076 
7 77.6 0.0960 0.1105 
8 79.1 0.0983 0.1121 
9 80.0 0.0996 0.1141 
10 80.4 0.1022 0.1169 
11 82.1 0.1073 0.1219 
12 83.9 0.1146 0.1303 
13 86.4 0.1266 0.1437 
14 86.9 0.1362 0.1532 
15 88.4 0.1463 0.1662 
16 90.6 0.1619 0.1827 
17 93.7 0.2140 0.2456 
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where K, is the thermal conductivity at (°C 
expressed in watt/em —°C. The average devia- 
tion, regardless of sign, of the ten experimental 
points from this equation is slightly less than 
0.4 percent. 

Below the freezing point of sodium three 
measurements, which are included in Table IV, 
were made. As a check may be cited the value of 
thermal conductivity obtained by Bidwell’ for 
sodium at 75°C of 1.13 watt/cm —°C. Extrapola- 
tion of the curve of Fig. 5 down to 75°C yields a 
value of 1.14 watt/cem —°C. 

If one assumes that the experimentally de- 
termined curves for the solid and the liquid state 
can be extended exactly to the freezing point, 
which shall be taken at 97.5°C, one gets 0.861 
watt/em —°C for the thermal conductivity of 
liquid sodium at the freezing point, and 1.120 
watt/cem —°C for solid sodium. The ratio of the 
thermal conductivity of the solid to the liquid 
state is thus 1.33. Porter and Simon!® obtained 
1.31 for the value of the same ratio. 


DISCUSSION 
Table V presents the results of calculations of 
the Wiedemann and Frantz ratio and of Lorentz’ 


TABLE IV. Thermal conductivity of solid sodium 
(cross section of sodium =18.70 cm?), 





MEAN | THERMAL 

Temp. Temp. GRADIENT POWER CONDUCTIVITY 

(°C) (°C/cm) (WaTTs) (WatTtT/cmM —°C) 
rr | — sail 

84.7 + 15 24.47 1.133 

88.2 1.27 26.78 1.128 

91.5 1.37 | 28.92 1.127 


TABLE V. Wiedemann and Frantz ratio and 
Lorentz’ number. 


METAL MERCURY SopIUuM 





Temperature, °C 0 120 200 100 200 
Thermal ¢ ‘onductivity, K 
(watt/cm —*C) 


| eaees 0.1093 | 0.1237 0.860 
Electrical Conductivity, 


0.808 


(ohm! cm~!) 10630 9484 8726 102,000 73,600 

Wiedemann and Frantz ra- 
tio, (K oX106 = watt- 
ohm/°C) 

Lorentz’ number, (K/oT 
X 105 watt-ohm/(°C)*) 


7.70 11.53 14.17 8.43 11.0 








2.82 2.93 2.99 | 2.26 2.32 





°C. C. Bidwell, Phys. Rev. 28, 584 (1926). 
1” Porter and Simon, Proc. Phys. Soc. 27, 307 (1915). 
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Fic. 5. The thermal conductivity of sodium. 


number using the values of electrical conductivity 
obtained by Williams"! for mercury, and listed in 
the International Critical Tables for sodium. 

Mercury is an exception to most metals, 
whether liquid or solid, in that its thermal 
conductivity increases with temperature. One 
can interpret this fact to mean that the effective 
number of free electrons increases with tempera 
ture. Mott" suggests that this would not be ar 
unreasonable behavior for a divalent metal; 
however, many divalent metals exist which are 
normal in behavior. If one fails to take this 
interpretation, one must assume that the lattice 
conduction in mercury is large and that it 
increases with temperature. Such a view is 
inconsistent with the theory developed by Bloch, 
and others. 

The thermal conductivity of liquid sodium 
decreases in a normal manner as the temperature 
increases. But the liquid sodium amalgams are 
abnormal in their behavior, even when of high 
sodium content. In fact the conductivity in- 
creases with temperature at a more rapid rate as 
the sodium content increases. 

In conclusion, the author wishes to express his 
thanks for the cooperation which was extended 
by the members of the department of physics at 
the University of Kansas, and to Dr.C. V. Kentin 
particular, who offered many useful suggestions. 


1 F, I. Williams, Phil. Mag. 50, 589 (1925). 
2 N. T. Mott, Proc. Roy. Soc. 146, 465 (1934). 
8 Bloch, Zeits. f. Physik 59, 208 (1930). 
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Internal Friction in Solids 
V. General Theory of Macroscopic Eddy Currents 


CLARENCE ZENER 
College of the City of New York 
(Received April 30, 1938) 


In vibrating ferromagnetic metals macroscopic eddy 
currents tend to shield the interior of the metal from 
changes in magnetic induction. The dissipation of energy 
by these eddy currents contributes to the internal friction 
of the metal. This internal friction has previously been 
investigated theoretically only for longitudinal vibrations, 
and only the asymptotic formulae for high and low fre- 


quencies have been found. In this paper this internal 


§1. INTRODUCTION 


ART of the internal friction of ferromagnetic 
metals arises from the magnetic-elastic 
coupling. Under certain conditions this part is 
of a larger order of magnitude than the internal 
friction of nonferromagnetic origin. In spite of 
numerous experimental investigations,' only one 
source of ferromagnetic internal friction has been 
theoretically investigated? for small strains, 
namely that due to the macroscopic eddy 
currents. These eddy currents tend to shield the 
interior of the sample from changes in magnetic 
induction. Further, this effect has only been 
investigated for longitudinal vibrations, and 
only the asymptotic values for high and low 
frequencies have been obtained. In this paper 
the internal friction due to macroscopic eddy 
currents is calculated at all frequencies for the 
longitudinal and transverse vibrations of circular 
rods, and for the transverse vibrations of reeds. 
It is found that the internal friction due to 
the eddy currents may be written in the form 


Q*=[(Es—Ex)/Ez |f(»). (1) 


Here Q™' is (1/27) times the fraction of the 
vibrational energy dissipated per cycle. Eg and 
Ey are Young’s moduli for constant magnetic 
induction and constant magnetic field, respec- 


1E. Giebe and E. Blechschmidt, Ann. d. Physik 11, 905 
(1931); O. v. Auwers, Ann. d. Physik 17, 83 (1933); J. 
Zacharias, Phys. Rev. 44, 116 (1933); W. T. Cooke, Phys. 
Rev. 50, 1158 (1936); Mary D. Waller, Proc. Phys. Soc. 
50, 144 (1938). 

2M. Kersten, Zeits. f. tech. Physik 15, 463 (1934); 
W. F. Brown, Phys. Rev. 50, 1165 (1936). 


friction is calculated for all frequencies for both longi- 
tudinal and transverse vibrations. The theory of internal 
friction due to macroscopic eddy currents is shown to be 
formally identical with the theory of internal friction due 
to macroscopic thermal currents. The methods developed 
by the author for the study of thermoelastic internal fric- 
tion are thus directly applicable to the study of the macro- 
scopic electric eddy currents. 





tively. The first factor may be calculated by 
means of Eq. (14). The second factor is given 
in §3 for the three cases, and is plotted in Fig. 1. 
The asymptotic expressions for f at high and 
low frequencies are shown in Fig. 1 by broken 
lines, and are likewise given explicitly in §3. 
The frequency of maximum internal friction is 
determined by the magnetic diffusion constant 
Dy. This is given by 

Dy =108/(0.4ry0), (2) 


where x is the differential magnetic permeability, 
and the electric conductivity ¢ is in ohm™ cm™. 
The area beneath all three curves is the same, 
namely 37 logis e=0.682. This is a direct conse- 
quence of the general formula 


[ foyrar=x/2, (3) 
Jo 


valid for a specimen of arbitrary cross section in 
which the only stress associated with the vibra- 
tion is a tensile stress along the axis of the 
specimen. 

The method used in the present investigation 
is that developed by the author for studying 
the internal friction due to macroscopic thermal 
currents.* The differential equations for the two 
problems become identical when 7 and S are 
replaced by H and B/47, respectively. The two 
problems would be identical if the adiabatic 
boundary condition for thermal flow were re- 
placed by an isothermal boundary condition. 


3C. Zener, Phys. Rev. 52, 90 (1938). 
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Thus the thermoelastic internal friction due to 
macroscopic thermal currents would, if the 
surface were maintained at constant tempera- 
ture, be given by Eq. (1) and Fig. 1, the sub- 
scripts JJ and B being replaced by 7 and S, 
respectively. 


§2. GENERAL THEORY 


The theory for longitudinal and transverse 
vibrations is particularly simple when the wave- 
length of vibration is large compared with the 
transverse dimensions of the specimen. In the 
first place, all stresses are then negligible com- 
pared with the tensile stress parallel to the 
specimen’s axis. If this axis is denoted as the 
axis, the energy loss per cycle per unit volume 
is then 


Aw =v-' time average of Z.de.,/dt. (4) 


In the second place, the transverse components 
of the oscillating part of the magnetic field and 
magnetic induction may then be neglected com- 
pared with the parallel components. Thus Z, 
may be written as a linear function of only e., 
and AH,, the fluctuation of //, from its mean. 


Z.= (02. 0€ 22) H(2)Czz + (OZz, O17.) e¢2n All:. (5) 


We are here neglecting a term proportional to 
the temperature fluctuation. The effect of this 
thermoelastic term has already been investigated 
by the author. 

By substituting (5) into (4) we obtain 


Aw =v"(0Z./0H-z) 22) 
time average of A//.de../dt. (6) 


We must now find the differential equation 
which A//, obeys. This is obtained of course 
from Maxwell’s equations. Here the displace- 
ment current may be neglected compared with 
the real current, since their ratio is given by 
eg/(1.8X10"%c), where o is the electric con- 
ductivity in ohm cm™. Elimination of the 
electric field from Maxwell's equations then gives 


(d/dt)B = —(108/0.42¢) curl curl H. 


If we take the zs component of this vector 
equation, and use our previous approximation of 


neglecting J/, and J/,, we obtain 


(d/dt)B,=(10%/0.420)(0?/0x°+0?/dy")/7.. 
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Fic. 1. f functions associated with macroscopic eddy cur- 
rents. Curve A, transverse vibration of reed. s=Dy/2ra?, 
d=transverse width of reed. Curve B, longitudinal vibra- 
tion of circular rod. z=D,/2ra?, a=radius of rod. Curve 
C, transverse vibration of circular rod. z=Dy/2ra?, a= 
radius of rod. 


On introducing the relation 
6B,= (OB, O1T .) ¢22)611,+ (0B, O@--) 112) 5€ 
we obtain 


(d/dt)AH, = Dy(0?/dx°+0?/dy*)AlT, 
+ (0/1,/0e.-)p .de../dt. (7) 


Here Dy, is defined by Eq. (2) with 
u=(0B, OFT .) (22). 


The first term on the right side represents the 

change in JJ, arising from diffusion of the 

magnetic induction. The second term represents 

the change arising solely from the change in 

strain. This would be the only term in the 

absence of diffusion, i.e., if the electric con- 

ductivity were infinite. 
The function de,,/dt is affected only very 

slightly by the magnetic-elastic coupling. Hence 

very little error is introduced in Aw by regarding 

de.,/dt in Eqs. (6) and (7) as specified by the 

elastic equations in the absence of this coupling. 

The function A//, is then completely determined 

by Eq. (7) and by the appropriate boundary 

condition. 
We have previously assumed the wave-length 

of the vibration to be long compared with the 

transverse dimensions of the sample, and also 

the transverse components of the oscillating 

part of H7 and B to be negligible compared with 

the parallel component. To this same approxi- 

mation we may neglect A//, outside the speci- 








1012 


men. The appropriate boundary condition is thus 
A/lT,=0 at surface. (8) 


Equations (6), (7), (8) contain the physics of 
our problem. The mathematics involved in the 
solution of Eq. (7) subject to the boundary 
condition (8), in taking the time average in 
Eq. (6), and finally in evaluating the integral 


f wdv 


O'= - 


27 energy of vibration 


has been formally solved in the appendix of 


reference 3. The result is 


Q '=[(E,—En) En l|> fi VEV (.7+*), (9) 
where f;, satisfies the condition 
Diem. (10) 


The constants f; and »;, are interpreted in terms 


of the differential equation 


{Dy (02/dx2+02/dy2) +2rv%| Us(x, y)=0 (11) 
together with the boundary condition 
U,.=0 at surface. (12) 


The », are the eigenwert of this equation. The 


f, are the squares of the coefficients of the 
expansion of the normalized function 
,; 
(13) 


tu /| fontéeas] 


in terms of the normalized eigenfunctions U,. 
Here e.. is taken at an arbitrary z and 1. 

The ratio (Ez,—Ey)/Ex may be calculated in 
terms of the magnetostriction constant and the 
permeability by means of the equation? 


(Ep—En)/En 

=4rkK1(de.- OIT..)2z:-)*/(OB. OLT.) 22). (14) 
This equation may be derived in exactly the 
same way as the equation 


(Es—Er)/Es=Er(0e::/0T) 2(2)"/(0S/0T) 2:2). 


§3. APPLICATION TO SPECIAL CASES 
In the previous section we found that the 
Q- of a rod of arbitrary cross section vibrating 


CLARENC 


E ZENER 


either longitudinally or transversely is given by 
(E,y—Ey)/E, times the function 


f(v) =frrnv/(ve2 +r). (15) 


In this section we describe in detail the evalua- 
tion of this function for one particular case, the 
transverse vibration of a circular rod, and then 
give the results for the longitudinal vibration of 
a circular rod, and for the transverse vibration 
of a reed. 
Circular rod: transverse vibration 

Let a be the radius of the circular rod. Then 
those eigenfunctions of Eq. (11) which satisfy 


the boundary condition (12), and whose asso- 


ciated coefficients in the expansion of 
é€--=constant Xr cos ¢ 
do not vanish, are given by 
cos oJ i(q.r/a). 
Here gq; is the kth root of 
J1(q) = {), 
and determines the eigenwert »;, by the equation 


gi = (27, Dy)?*a. (16) 


We may now write f; explicitly as 


9 


(| rT i(qur adr) 
“0 


f¢=—™- 


| rar | rJ\*(qiur/a)dr 
Ss 0 


With the aid of the integral formulae for Bessel 
functions* we obtain, as in reference 3, §2, 


f= 


Since J,(q.)=0, the recurrence formulae for 
Bessel functions® gives 


J2(qQx) = —Jo(q ). 


Hence fic=8/ 12. (17) 
See E. Jahnke and F. Emde, Funktionentafeln (Teub- 
ner, 1933), pp. 213-214. 
5See G. N. Watson, Theory of Bessel Functions (Cam- 
bridge, 1922), p. 17. 
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As a check upon Eq. (10), we note that® 
>a? =}. 


On substituting Eqs. (16) and (17) into Eq. 


(15), we obtain 


f(y) =8s>0 (g2+2*)", (18) 


where 


The asymptotic expressions for small and large sz 
may readily be obtained. From® 


> 9 1 = 1/192 
we obtain 


when 72<1l. 


f(v) =2/24 


By observing that when 2 is large, successive k’s 
differ by 7, we obtain, for large z’s, 


yo 


f(v) = (82 n) | (q*+27)—'dz. 
0 


Hence f(v) =2's when =s>1. 


When z is of the order of magnitude of unity, 
f(v) may be evaluated by means of the rapidly 
converging series (18). In plotting this function 
in Fig. 1, the values of g, were taken from refer- 
ence 5, p. 748. 
Circular rod: longitudinal vibration 

The constant f; is given by 

Si =4 9°; 


® Reference 5, p. 502. 
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where q is the &th root ot 
J (g) =0. 


The constant », is related to q by keg (16 
The function f is given in terms of s=(2ra*®, Dy 


by 


|e Bs (q) '4+<°)~! exactly, 


l 


by 3/8 when s<2zr, 


by (2/s)* when s>2z 


These two asymptotic expressions have previ 


ously been obtained.” 


Reed: transverse vibrations 


Let d be the width of the reed in the plane of 
vibration. Then 


V_= RDy 2xd*. 
The function f is given in terms of s= (27rd*/Dy)v 
by 


62r—s > (k4+2°7)—' exactly, 
k=1 


when s<4r, 


by (3/rv2)s-? when 2>4r. 
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LETTERS TO THE EDITOR 


Prompt publication of brief reports of important discoveries in physics may be secured by 
addressing them to this department. Closing dates for this department are, for the first issue of the 
month, the eighteenth of the preceding month, for the second issue, the third of the month. Because of 
the late closing dates for the section no proof can be shown to authors. The Board of Editors does 
not hold itself responsible for the opinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


The Scattering of Neutrons by Deuterons 


The angular distribution of 328 tracks of recoil deuterons 
has been studied. These tracks have been selected from 
one thousand tracks measured on 3000 photographs, which 
were taken with the aid of an expansion chamber. The 
chamber was filled with 54.3 percent deuterium, 37 percent 
argon, 6 percent nitrogen, and 2.7 percent heavy water 
vapor; this mixture was bombarded with neutrons of 2.6 
Mev energy from the deuteron-deuteron reaction. The 
criterion for selecting the above tracks was established by 
using: (1) the calculated stopping power of the gas mix- 
ture, (2) the pressure variation within which tracks can 
form, (3) the energy-range relationship for deuterons,! and 
(4) the energy variation in the neutrons caused by (a) loss 
of energy of the incident deuteron beam in the thick 
heavy-ice target and (b) the angular aperture of the de- 
tector. To these selected tracks have been applied the 
necessary azimuthal, area, and random track corrections, 
in a manner similar to that used in the study of the scatter- 
ing of neutrons by protons.? 

The results are given in Table I and Fig. 1. 

TABLE I. Observations on the distribution and angle of the recoil deu- 
terons bombarded with neutrons. Ad is the scattering angle interval. N is 
the number of tracks observed in each Ad. C is the total correction factor. 
Ne ts the number of tracks, corrected. N-/(Aw*Ad) is the number of tracks 
per unit solid angle per unit Ad. % 
tracks per unit solid angle per unit Ad. 





Ad N Cc Ne N-/ (Aw: Ad) % / (Aw Ad) 

0-9 17 3.2 55 7050 25 
10-19 31 3.3 102 4600 16 
20-29 38 3.5 132 3310 12 
30-39 47 2.4 114 2280 8 
40-49 53 3.0 158 2525 9 
50-59 69 3.6 247 3470 12 
60-69 54 5.2 281 3560 13 
70-79 18 6.7 121 1447 5 


| 
| 
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Distribution and angle of deuterons recoiling from neutron 
collisions. 
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(Aw Ad) is the relative number of 


Although these results are based on a small number of 
tracks (328), the data indicate that there is present a 
marked deviation from spherically symmetrical scattering 
in the center of mass system. More data are needed before 
the precise shape and size of the peak in the curve (cf. 
Fig. 1) at scattering angle ¢=60° can be determined. 

The marked departure from S scattering in this experi- 
ment can possibly be explained in three ways. Massey and 
Mohr,’ assuming a rectangular potential hole of radius 
4X 10-" cm and one of 6 X10~-" cm, obtained curves which 
are not far different in shape from our curve, except at 
small angles. Their calculations include 1 Mev and 6.25 
Mev neutrons and take into account P scattering. It seems 
likely that the choice of a different potential hole with 
further regard for the structure of H** might allow agree- 
ment with this experiment. A second possibility,® to account 
for the marked difference from S scattering and the scat- 
tering predicted by Massey and Mohr at the small angles 
(0°-20°), would be to include the effect of a polarization 
of the deuteron by the incident neutron. A third possi- 
bility is to include higher orders of partial waves in the 
description of the scattered neutron wave. 

From the stopping power of the gas in the chamber, 
the range-energy relation, and the maximum track length 
extrapolated in the usual way in the scattering angle 
intervals 0°-10°, 0°-20°, 0°-30°, the range of a recoil 
deuteron at 0° scattering’angle has been calculated. This 
has been used to determine the disintegration energy VU 
for the reaction 


1, D?+ ,D?—,He?+ on’ +Q 


and leads to a value of Q=3.40+0.10 Mev. This value, 
because it is based on deuteron recoil atoms, is an inde- 
pendent check on the value recently given by Bonner,® and 
agrees very well. 

A more complete description of this experiment will be 


ready for publication shortly. 
P. G. KRUGER 
W. E. SHoupp 
R. E. Watson 
F. W. STALLMANN 
Physics Department, 
University of Illinois, 


Urbana, Illinois, 
May 17, 1938, 


1H. Bethe, Phys. Rev. 53, 313 (1938). 

2 P. Gerald Kruger, W. E. Shoupp, and F. W. Stallmann, Phys. Rev. 
52, 678 (1937). 

3H. W. S. Massey and C. B. O. 
(1935). 

4L. I. Schiff, “Scattering of Neutrons by Deuterons,’’ Phys. Rev. 52, 
149 (1937). 

5’ This idea was suggested by Dr. J. H. Bartlett, Jr. 

6 T. W. Bonner, Phys. Rev. 53, 711 (1938). 


Mohr, Proc. Roy. Soc. 148A, 206 
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LETTERS TO 


The Molecule of Sulphur and Its Raman Spectrum 


Warren and Burwell! have recently shown that rhombic 
sulphur consists of Ss molecules; the atoms in each molecule 
are located at the corners of a puckered octagon. We will 
now apply the methods based on group theory as developed 
by Wigner,’ Tisza,’ and others for investigating the normal 
modes of oscillation of such a molecule. There are sixteen 
symmetry operations which transform the model into itself 
and these operations form the point group Ss,. There are 
seven irreducible representations for this group and the 
character Table I obtained from the general tables given 
by Tisza is given. 

The significance of h;, x;’, n; and n,’ is the same as that 
in a recent paper by Wilson.‘ The activity or otherwise of 
these oscillations is also given in the last two columns of 
the table. P indicates a strong and well-polarized Raman 
line. D indicates a completely depolarized Raman line and 
f stands for a forbidden line. The table shows that we 
should expect two strong and well polarized lines and five 
depolarized lines in the Raman spectrum of sulphur. 

Venkateswaran® has recently obtained the complete 
Raman spectrum in which eight lines at 88, 114, 152, 185, 
216, 243, 434 and 468 have been recorded. He has given 
good reasons to show that the lowest frequency is to be 
attributed to the lattice and not to the molecule. The 
molecular frequencies are therefore only seven and two 
amongst these, namely 216 and 468, have been found to be 
strong and well polarized in complete accord with our 
These two represent the total symmetric 
vibrations A;. All the other five, namely 114, 152, 185, 243 
and 434 should be completely depolarized according to the 
and 434) have been 


conclusions. 


theory. The stronger of these (152 
found to be so but data in respect of others are not given 
by Venkateswaran presumably because they are weak. The 
Raman spectrum data so far as are available may thus be 
taken to be in complete accordance with the model pro- 
posed by Warren and Burwell. Interpretation of the infra- 
red spectrum and a more detailed analysis of the vibra- 
tional modes and frequencies will be given elsewhere. 


S. BHAGAVANTAM 
Department of Physics, 
Andhra University, Waltair, India, 
April 13, 1938. 


1 Warren and Burwell, J. Chem. Phys. 3, 6 (1935). 

2? Wigner, Géttinger Nachrichten, 133 (1930). 

3 Tisza, Zeits. f. Physik. 82, 48 (1933). 

‘Wilson, Phys. Rev. 45, 706 (1934). 

6 Venkateswaran, Proc. Ind. Acad. Sci. 4, 345 (1937). 
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Radioactive Iodine Isotopes 


We have bombarded tellurium with 8 Mev deuterons 
and have found two new radioactive iodine isotopes, with 
half-lives 
two mechanisms: 
(2) by direct transmutation from stable tellurium. 


of 13 hours and 8 days. The latter is created by 
(1) by the decay of radio-tellurium and 
Process 
(1) is demonstrated by the fact that successive extractions 
of iodine, from the same solution of bombarded tellurium, 
show a growth of the 8-day period; this must therefore be 
associated with either [9 or ['%!, 
a double decay: Te!®: °°(d, p)Te!9. 3; Te! 
—Xe!??. 131. Process (2), 
basis of relative intensities, is the usual type of reaction: 
Te. 13°(d, nm). 31, Conclusive proof of this identification 
and interpretation has been furnished by the extraction of 


as the second activity of 
, dl» ] 129, 131 


which is known to occur on the 


the 8-day iodine from tellurium which has been bombarded 
with neutrons. We cannot as yet state the period of the 
radio-tellurium from which this iodine grows. 

Absorption measurements on the negative electrons 
emitted by the 8-day iodine indicate a maximum energy of 
0.9 Mev; a gamma-ray is also present. 

We have irradiated iodine with the fast neutrons from 

lithium plus deuterons source and confirm the 13-day 
period reported for the same bombardment by Tape and 
Cork,' who surmised it to be due to I'**, We have chem- 
ically identified this activity as an iodine isotope, so that 
the assignment to I'*° appears to be definite. (The antimony 
fraction of this same bombardment was inactive, while the 
tellurium precipitate exhibited a 10-hour half-life which 
can be ascribed definitely to Te’, This tellurium period 
has been reported previously,' following neutron and deu- 
teron bombardment of tellurium, but without definite 
isotopic identification.) 

The yield of the 8-day iodine from Te(d,)I is very 
much larger than that of the strongest 13-day iodine that 
we have been able to produce by the reaction I'*"(m, 2m) 1"*°, 

This research has been aided by grants from the Re- 
search Corporation, The Chemical Foundation and the 
Josiah Macy Jr. Foundation. 

J. J. Livincoop 
G. T, SEABORG 
Radiation Laboratory, Physics Dept. (J. J. L.) 
Cc hemistry De partment (G, S.), 
University of California, 


Berkeley, California, 
June 1, 1938. 


1G. F. Tape and J. M. Cork, Phys. Rev. 53, 676 (1938). 





TABLE I. 

E 28! 2c? 23 Cc 4e, 4C2 nj R&T ny’ RAMAN _ INFRA-RED 
Ai 1 1 1 1 1 1 1 2 2 P j 
Ao 1 1 1 1 1 —1 —1 1 R, 0 
Bi 1 —1 1 —1 1 1 —1 2 7; 1 f M,=0 
Bz 1 1 1 -1 1 -1 1 1 1 f f 
Ei 2 v2 0 —yv2 —2 0 0 3 Ty 2 J M,=0 
E: 2 0 —2 0 2 0 0 3 Ry 2 dD f 
E3 2 —7V2 0 v2 —2 0 0 3 3 D f 
hj 1 2 2 2 1 4 4 
x;' 24 0 0 0 0 2 0 
hjx;" 24 0 0 0 0 x 0 
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Radioactivity of Be’ 


As reported at the Washington meeting of the American 
Physical Society! we have obtained evidence for the 
formation of Be’ in the reaction Li®'+H?—~Be’+n. At 
that time we had found that there was no radioactive 
emission of charged particles. We have since discovered a 
gamma-ray radioactivity with a 43-day half-life in lithium 
targets which have been bombarded by deuterons. Further 
experiments indicate that the gamma-ray activity origi- 
nates in K electron capture by Be’ which leaves the 
residual Li? nucleus in the 450-kv excited level. 

The activity was first observed in a LiF target which 
had been bombarded a month previously with 35 micro- 
ampere hours of 1000-kv deuterons. Similar activity was 
subsequently produced by bombarding a lithium metal 
target with 20 microampere hours of 1000-kv deuterons. 
The activity was followed for a month and showed a 
half-life of 4326 days. The absorption-coefficient in lead 
was 0.169 g cm™!. No charged particles could be detected 
either by thin-walled counters, by ionization chambers, 
or by placing the sample in a cloud chamber, although the 
gamma-ray counting rate indicated the emission of more 
than 100 quanta per second. 

The surface of the target (in which the activity was 
concentrated) was scraped off and subjected to a chemical 
separation, carried out by Dr. M. H. Van Horn of George 
Washington University. The activity was concentrated in 
the final precipitation of 100 micrograms of beryllium 
carbonate which had been added as a carrier. The separa- 
tion would not, however, remove iron or aluminum. As 
iron also shows a 40-day period, we bombarded an iron 
target and looked for a similar gamma-ray. The iron 
impurity necessary to produce the observed gamma-rays 
from the lithium targets is greater than 100 percent. 

Although we were not able to detect the alpha-particle 
group produced in the reaction! B!°+H'!—Be’+ Het, we 
have since observed a long period gamma-ray activity 
from a B,C target which had been bombarded with 11 
,microampere hours of 950-kv protons. The intensity of 
the radiation was 25 percent less than that produced in 
the lithium metal target. Its absorption-coefficient was 
also 0.169 g cm™!. We therefore ascribe both this activity 
and that produced in the lithium targets to Be’. 

This unusual type of radioactivity could originate in 
two nuclear reactions. Either a positron is emitted with 
so little energy that it does not ionize appreciably and is 
detected only by the annihilation radiation, or a K electron 
is captured and the residual nucleus is left in an excited 
state. In the latter case the average number of quanta 
per disintegration is less than one because there exists the 
alternative process of capture in which the nucleus is left 
in the ground state. Assuming as usual that a neutrino is 
emitted, these reactions may be written: 


Be’—Li'+e*+7 
et+e-—hyt+hy 
Be’+e,-—Li’+7 
Be?+e,-—Li™*+7 
Li7*—>Li?+hyp 


hvy=0.5 Mev 
neutrino energy = E Mev 
neutrino energy = E—0.45 Mev 

hv=0.45 Mev 


THE EDITOR 

The absorption-coefficient in lead of the gamma-rays 
from the lithium targets and of the annihilation radiation 
of the positrons from N™ were measured under idential 
conditions and were observed to be 0.169+0.008 and 
0.130+0.002 g cm“, respectively. This indicates that the 
energy of the radiation is 425+25 kv, which agrees well 
with the energy expected from K electron capture but 
differs from the energy of annihilation radiation. 

A comparison of the numbers of quanta observed per 
neutron emitted from the carbon and lithium targets gives 
additional information about the process. Li® and Li? are 
(per atom) equally effective in producing neutrons at 
800 kv and the yield curves for the two reactions are 
known.” The energy distribution of the neutrons from Li® 
shows that approximately one atom of Be’ is formed for 
every two neutrons emitted. In carbon one neutron is 
emitted for each N® atom formed. The comparison of the 
number of neutrons emitted by lithium and carbon has 
been made by Amaldi, Hafstad, and Tuve.* Thus the 
relative numbers of Be? and N® atoms can be calculated. 
Assuming two quanta per N"® we find that, on the average, 
one quantum is emitted for each 10 Be’ atoms produced. 
Positron-emission should lead to a gamma-ray intensity 
20 times greater than that observed. On the other hand 
the data are consistent with the hypothesis of K electron 
capture, because in this case transitions to the ground 
state should be more probable than transitions to excited 
levels. 

As a final check we have looked for double coincidence 
counts. No appreciable change in the background rate of 
0.8 per minute was observed although the simultaneous 
emission of two quanta at 180°, which would be expected 
if positron annihilation occurred, should increase the rate 
by at least a factor of three. 

R. B. RoBerts 
N. P. HEYDENBURG 
G. L. LoOcHER 


Department of Terrestrial Magnetism (R. B. R. and N. P. H.), 
Carnegie Institution of Washington, 
Washington, D. C.., 
Bartol Research Foundation of the Franklin Institute (G. L. L.}, 
Swarthmore, Pennsylvania, 
May 27, 1938 


1R. B. Roberts and N. P. Heydenburg, Abstract No. 78, 1938 
Washington Meeting American Physical Society. 

2 L. H. Rumbaugh, R. B. Roberts, and L. R. Hafstad, to be published 
shortly in Phys. Rev. 

+E. Amaldi, L. R. Hafstad, and M. A. Tuve, Phys. Rev. 51, 896-912 
(1937). 





The Transmutations of the Cosmic-Ray Electrons and the 
Nuclear Forces 


The “‘transmutations of mass” (i.e., the changes of the 
rest mass) of the electrons probably play an important role 
in the passage of cosmic rays through the atmosphere. Thus 
according to Blackett! the cosmic-ray particles at sea level 
are heavy when energetic, but change their mass suddenly 
when slowed down below a critical energy of about 2-3 - 108 
ev and become ordinary electrons. According to Bowen, 
Millikan and Neher? the penetrating heavies are mainly 
secondaries produced in the atmosphere by incoming 
electrons and photons. However this conclusion, does not 


necessarily follow from the data of the authors mentioned; 
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these data are consistent also with the assumption of a 
primary character of the penetrating heavies. 

No adequate mechanism has been suggested up to now 
for these transmutations of mass. The suggestions of 
Jauncey’ and Bhabha‘ can be shown to be untenable. The 
creation of heavy pairs by photons does not account for the 
disappearance of the slowed down heavies. Youkawa's 
hypothesis of a spinless heavy electron seems to lead to 
difficulties® and will not be considered by us. 

Considerations connected with the principle of detailed 
balance and also with the probable dependence of the 
cross section for transmutation on the atomic weight of the 
substance traversed, seem to indicate, that the heavies are 
stable in absence of atomic nuclei and that the transmuta- 
tion of mass can be induced only by a nonelectromagnetic 
short range interaction of an electron with protons and 
neutrons. 

Let us assume, that this interaction is of the 6-type and 
corresponds to terms of the form 


vey (1) 


in the Hamiltonian, where ¥; and y, are the wave functions 
of electrons of mass kmy and Imp, respectively, and V the 
wave function of a proton or a neutron. 

Consider an electron of mass m, and energy E hitting 
a nucleus and scattered by it with a modified mass my. If 
the recoil of the nucleus is neglected 

E=c(m,2c?+ p;?)4=c(me2c?+ p,?)}, (2) 

where p; and fp» are the momenta of the incident and the 
scattered electrons.® 

The cross section o for the transmutation m.—>m, differs 
from that for the transmutation m:—>mz, by the factor 


pr? F—m,rct 


—_ » P ’ 
p? E?—m2ct 





so that it is sufficient to consider only the case for which 
(for instance) m,<mz. In this case o vanishes for E <myc? 
and then increases with increasing £ somewhat faster than 
E*, But when the wave-length of the electron becomes 
comparable with the nuclear radius (EF ~A7~!-7-107 ev, 
where A is the atomic weight) the interference of the partial 
waves scattered by different particles of the nucleus makes 


“a 


« tend with increasing FE to a constant “stationary” value 
o,t and also ensures the preponderance of the forward 
scattering. o.¢ is proportional to A‘, But if E exceeds 
about 2-10* ev the recoil of the nuclear particles can no 
longer be neglected, since the recoil energy exceeds their 
binding energy (if the angle of scattering is not very small). 
In consequence ¢ will again begin to increase in proportion 
to E*, Hence, as long as o has at E= 10" ev any reasonably 
moderate value, its value at E<5-10* ev must be quite 
negligible. 

However, in the domain of very high energies the validity 
of the above calculations becomes questionable. It is 
known, that interactions of the type considered lead to 
divergent results even in such approximations, which can 
be satisfactorily treated by the present quantum electro- 
dynamics. Thus the interaction under consideration implies 
the possibility for a neutron or a proton to create a pair, 
the electron and the positron of the pair having eventually 
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different masses. Hence there must exist certain forces 
between neutrons and protons which correspond to the 
exchange of pairs between these particles. But if one tries 
to evaluate in the first approximation the energy #(r) 
of this interaction one gets a divergent result, even if the 
separation r of the interacting particles is finite.’ 

In absence of any better theory let us introduce in the 
interaction operator (1) some “guillotine factor,” e.g. the 
factor e~7"et PY); py and p; are the momenta of the electrons 
of mass kmty and Imo, measured in that reference system, in 
which the nucleus is at rest. This factor makes the energy 
J(r) finite: at large distances 3(r) varies as r~*, but tends 
to a finite value 3(0) as r tends to zero, 

But the guillotine factor will affect also the cross section 
for the scattering of electrons with transmutation of mass. 
No such factor is needed in the case of the Coulomb scat- 
tering since for the latter only distant encounters are 
essential, while the guillotine factor accounts in some very 
crude manner for the finite size of the interacting particles, 
On account of this factor the probability of transmutations 
of energetic electrons will become very small, but, if 
7=400 myc, it will rapidly increase as the electrons are 
slowed down to about 3-108 ev and at E~2-108 ev o will 


“a ’ 


reach the “stationary’’ value defined above. This is in 
accordance with Blackett’s measurements. 

With this value of r one gets for the energy of inter- 
action of nuclear particles a reasonable estimate (0) ~ 
—3N?(o,,2ir- 1075) X 10° ev where o,;*'7 is the 


value of o in air expressed in cm*, and N the number of 


“ 


stationary” 


allowed values of the mass km of an electron. The sign 
and the spin dependence (if any) of 3(0) are correct. 

The possibility of explaining the nuclear forces by the 
exchange of pairs between neutrons and protons has al- 
ready been suggested by Gamow and Teller’ and by Went- 
zel,? who, however, considered only ordinary electrons 
(m=my) and were unable to establish any connection 
between the proposed explanation of nuclear forces and 
other physical phenomena. As pointed out by these authors, 
the forces due to the exchange of pairs do not depend on the 
electric charge of the nuclear particles and are consistent 
with Wick’s explanation of the magnetic moment of the 
neutron. The only serious difficulty is that they apparently 
do not possess saturation properties. 

A more detailed communication will appear in Comptes 
Rendus Acad. Sct. U. S. S. R. 

Ic. TAMM 

Physical Institute of the Academy of Sciences of U.S. S. R., 


Moscow, U.S. 
May 3, 1938. 


' Blackett, Proc. Roy. Soc. 164, 1 (1938) 

2 Bowen, Millikan and Neher, Phys. Rev. 53, 217 (1938). 

3 Jauncey, Phys. Rev. 52, 1256 (1937) 

4 Bhabha, Proc. Roy. Soc. 164, 257 (1938). 

5’ Serber, Phys. Rev. 53, 211 (1938). 

6 In lectures held in Moscow and Leningrad in February I suggested 
that some of the cloud chamber tracks which apparently violate the 
conservation of energy may be explained on the basis of Eq. (2) 
Afterwards I learned that Zwicky also suggested a similar explanation 
cf., Phys. Rev. 53, 315 and 611 (1938). 

? Although the integration may be arranged in such a way as to give 
a finite result at finite r’s, still if one performs the integrations in a 
somewhat different succession the integrals do diverge even at finite r's 
The reason at the bottom of this is that the interaction corresponds to 
the exchange of pairs of particles, and not of separate particles. 

* Gamow and Teller, Phys. Rev. 51, 289 (1937). 

* Wentzel, Helv. Phys. Acta 10, 107 (1937). 
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Excitation of Nuclei by Bombardment with Charged 
Particles 


Charged particles are able to excite nuclei without pene- 
trating into the nucleus by means of the action of their 
electric field upon the nucleus when they pass nearby. 
The cross section for such a process can be calculated and 
is surprisingly large. An approximate expression for the 
cross section for the excitation of a transition from a level 
A with the orbital momentum /=0 to another level B 
with /=L=2 is given by 
16? L 
™ (2L+1) (2L—1)(2L—2) 

1 = ( hv’ 2L “4(Oap)?* 

LE aot AX2L-3 
Here E and 2X are the energy and wave-length of the 
incident particle of charge ze; E’, 27X’, v’ are its energy, 
wave-length and velocity after the nuclear excitation. 
Ze is the charge of the nucleus, and Q4z is given by: 


(Qaz)®=LS ba*ri“dpdr, 
r 





OL 








(1) 


where 7; are the radial coordinates of the nuclear protons 
and ¢a and 92 are the eigenfunctions, respectively, of the 
nuclear states A and of one of the 2L+1 degenerate 
states B. 

The expression (1) is valid only so long as E is less than 
the potential barrier: E<Ey=2Ze?/p (p nuclear radius). 
If E is greater than that limit, oz, is to be multiplied by 


(2L—1) (/Ey\24? Eo 
aye) a 
(2L—2)\E (2L—1)E 
The exact formula for oz contains rather complicated ex- 
pressions with £, E’ and L and will be presented together 
with the derivation of (1) in a paper by M. S. Plesset and 
the present author which will appear shortly in the Physical 
Review. 

The cross section 2 for the excitation of one quadrupole 
transition with an excitation energy of 1 Mev in a nucleus 
Z=50 by means of a proton of 6 Mev is according to (1): 
o2=0.5-10-%* if one assumes Qag=5-10-%. The cross 
sections for L >2 are much smaller. The dipole transitions 
which are not included in (1) are also very weak because 
of the smallness of nuclear dipole moments. The numerical 
values obtained by (1) are extremely sensitive to the value 
of Qaz. The order of magnitude of Qaz is deduced from 
y-ray evidence. 

It is possible to estimate the total excitation cross sec- 
tion; i.e., the sum of all cross sections for any possible 
excitation in a heavy nucleus. If the density of levels with 
1=2 is given by w(£), the total cross section for excitation 
of a nucleus by means of a charged particle of an energy E 
is approximately (apart from the factor (2) if E> Ep) 


1 E” Qe 
total .. 61/25/28 —_. -—- 7 F) «a(R — Fe’ 
total = 51/2_ - + -T(E)-w(E-E’). 


Here E’=(5/2)7(E) and 2xX’ is the corresponding wave- 
length. The temperature 7(£) is given by 1/T=d lg w/dE 
and (E—E’) is the average excitation energy of the nucleus 
after collision. Q is to be considered as the average value 
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of the quadrupole moments corresponding to all possible 
transitions. One gets for the total excitation cross section 
of a nucleus with Z=50 with a-particles of 14 Mev: 
o,*°tal —3.10-25 and a mean excitation energy of 10 Mev. 
An expression is used here for the level density given by 
Bethe! and Weisskopf,? and Q is taken equal to 0.4-107%, 
This value of Q is computed from the radiation breadth 
of neutron capture levels. These results have to be con- 
sidered as a rough estimate of the order of magnitude 
because of the uncertainty of the value of Q and w(£). 
The process described here may be of some interest for 
the investigations of nuclear excitation levels. Further- 
more, after bombardment of nuclei with charged particles 
of high energy (E>8 Mev) one should expect that the 
excited nucleus emits a particle so that one gets processes 
of the kind (p—p, ), or (p—2p). This effect may be of 
some importance in bombardments of heavy nuclei; proc- 
esses such as (p—p,m), (d—d,n), (a—a,m) are not to 
be expected after penetration of the projectile into the 
nucleus since the re-emission of a charged particle by the 
compound nucleus is much less probable than the emission 
of neutrons.?:? These reactions could, however, be pro- 
duced in the way described here. 
University of Rochester, Rochester, N. Y. V. F. WEISSKOPF 
June 2, 1938 
{1. A. Bethe, Phys. Rev. 50, 332 (1936). 


1H. 
2V.F. Weisskopf, Phys. Rev. 52, 295 (1937). 
3E. J. Konopinski and H. A. Bethe, Phys. Rev. in press. 





Heavy Electrons and §- Decay 


Yukawa! has shown that nuclear forces of the correct 
magnitude and range can be accounted for by assuming 
the existence of new charged particles of mass about 137 
times that of an electron and of zero spin. These heavy 
electrons have recently been observed? in expansion 
chamber photographs of cosmic rays and have consequently 
brought Yukawa’s theory to general notice. If the heavy 
electron were formed from an electron and a mass-less 
neutrino it should be extremely unstable and should 
dissociate with the liberation of energy equal to 136 mc’. 
It is difficult to see how it could ever exist long enough 
to be observed in cosmic rays. We can get over this 
difficulty by assuming the heavy electron to be formed 
from an electron and a “heavy neutrino’ of mass 136 
and spin }. It might now be possible to do away entirely 
with the conception of the mass-less neutrino, for it is 
evident that the emission of one of these heavy electrons 
obeying Bose statistics in 6-disintegration would satisfy 
the statistical requirements of the nucleus. If this heavy 
electron then disintegrated into an electron and a heavy 
neutrino while still close enough to the nucleus to interact 
with it the continuous energy distribution of the §-rays 
would be qualitatively accounted for. 

Bhabha’ has pointed out that if the disintegration of 
the heavy electron is spontaneous it may be regarded as a 
“clock,” and hence it follows merely from relativity 
considerations that the time of disintegration is longer the 
greater the velocity of the particle. A heavy electron of 
mass 137 formed as suggested above from an electron and 
a heavy neutrino of mass 136 would be just on the point 
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of instability since it would have no mass defect. It is of 
interest, therefore, to see if its much greater velocity in 
cosmic rays would increase its time of dissociation suff- 
ciently to enable it to be observed. 

If t% is the time of disintegration of the particle at rest 
and ¢ the time when it is moving with velocity v={c 
we have 

t=to/(1—B*)'. 
From the relativistic energy equation 
eV =300 mce[(1—s*)-$—1], 


it follows that the relative times of disintegration of two 
particles having energies V; and V2 e-volts is given by 
, V2+300 me? /e 
,= foo - ay ’ 
V,+300 mc?/e 








where m is the mass of the particle, here equal to 137 
times the mass of an electron, and 300 mc?/e=7 X10’. 
The relative distances traveled before disintegration is 
then given by d2=dyvele/mfh. 

Taking the mean energy of the end-points of 8-decay 
spectra to be 10° e-volts I find that the ratio of the distance 
traveled by a cosmic heavy electron of 10" e-volts energy 
to that traveled by a #-decay heavy electron before 
spontaneously disintegrating is 10°. To account for the 
continuous energy distribution of B-rays the heavy electron 
would have to dissociate while still close enough to the 
nucleus to interact with it, so that this figure of 10° is 
hardly great enough to allow the heavy electron to be 
observed in cosmic rays, unless the interaction with the 
nucleus caused the slower §-ray heavy electron to be 
dissociated in a time much shorter than its time of spon- 
taneous disintegration. Crane‘ has reported that in the 
course of some cloud chamber experiments a few 8-particles 
appeared to behave in an anomalous way for which no 
satisfactory interpretation could be found. Zwicky® also 
has drawn attention to some peculiar cloud chamber 
tracks which, in his view, may be due to the electron 
suddenly changing its rest mass. This is just what would 
be observed if a #-decay heavy electron disintegrated 
spontaneously outside the atom, the probability of dis- 
sociation of the particle by interaction while it is still 
close to the nucleus is not quite unity. 

Corben® has suggested that the heavy electron may be 
formed by the combination of an electron with one of 
Eddington’s neutral particle (} of a scalar particle) which 
has a mass 135.9. He considers this occurs with the emission 
of a neutrino. In Eddington’s theory,’ however, the neutral 
particles have no objective existence: they form a back- 
ground consisting of the unspecified particles of the 
universe. The highest energy level of this background 
represents the ground level of the “object system” to 
which the vector wave functions of quantum theory apply: 
the scalars apply to the background which forms the 
reference frame. The neutral particles of the background, 
which are all below the ground level of the object system, 
are therefore in negative energy states. It is a vacancy or 
hole in the background that manifests itself as a particle 
in the object system. , 

If the neutron—proton transition in the nucleus interacts 
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with the background as I have suggested,* and causes an 


electron to be created, the hole left in the background 
will represent the creation of a neutral particle with spin }, 
which could then combine with the electron to form a 
heavy electron which obeys Bose statistics, as is suggested 
above. The hole must represent an uncharged particle of 


spin } for charge and spin to be conserved among the 


particles of the object system. On this view the heavy 
electron is created by the neutron—proton transition in 
the atom. 


The University, F. L. Arnot 
St. Andrews, Scotland, 
May 14, 1938. 


1 Yukawa, Proc. Phys.-Math. Soc. Jap. 17, 48 (1935) 

2 Neddermeyer and Anderson, Phys. Rev. 51, 884 (1937) and others. 
3 Bhabha, Nature 141, 117 (1938). 

‘Crane, Phys. Rev. 53, 317 (1938). 

6 Zwicky, Phys. Rev. 53, 611 (1938). 

* Corben, Nature 141, 747 (1938). 

7 Eddington, Relativity Theory of Protons and Electrons. 

8 Arnot, Nature 139, 1065 (1937). 





Some Results of the Search for Super- Novae 


Some time ago Baade and I discussed the existence of a 
new class of novae which surpass the common novae in 
luminosity by a factor of one thousand.':? We proposed 
to call these new stars super-novae, a designation which is 
now in general use. It should be emphasized, however, 
that super-novae in some respects differ fundamentally 
from ordinary novae. 

In our original communications':? we made a_ first 
tentative attempt to estimate how often super-novae 
appear in an average nebula. From a considerable, though 
very heterogeneous set of records from various sources 
extending over the past fifty years we concluded that the 
frequency of occurrence of super-novae in an average 
nebula is of the order of one per several centuries. We also 
suggested that the study of super-novae promises to 
throw new light on the problem of the generation of energy 
in stars and perhaps on the origin of the cosmic rays. In 
conjunction with considerations concerning the tremendous 
liberation of energy in super-novae we suggested the 
formation of collapsed neutron stars: ? as the most powerful 
source of energy. 

These suggestions clearly reveal the great potentialities 
of a study of super-novae. We decided that no effort 
should be spared to track down and to study in detail as 
many as possible of these rare objects. Several seasons of 
unsuccessful work with unsatisfactory telescopic equip- 
ment preceded a period in which the greatest advance was 
made with the aid of the 18-inch f : 2 Schmidt telescope 
which had been built in the meantime with the authoriza- 
tion of the observatory council of the California Institute 
of Technology. 

In the period from September 5, 1936 until January 31, 
1938, about six hundred excellent photographs covering 
seventy square degrees each were obtained with the 
Schmidt telescope. The resulting search is equivalent to 
the continuous control during one year of about 1800 
nebulae whose apparent brightness is greater than m=15., 
Since three super-novae were discovered in the period 
mentioned, the resulting frequency is one super-nova per 
average nebula of our collection in a period of six hundred 
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years, a result which confirms the correctness of our 
original estimate. 

The two brightest of these super-novae, which appeared 
in N.G.C. 1003 and I.C. 4182, have now yielded a vast 
amount of information which will be presented in a series 
of papers to the Astrophysical Journal. Summarizing some 
of these results we may state: 

1. For the first time excellent photographic light curves 
have been obtained for super-novae, including the maxima. 
A large number of photographs taken with yellow filters 
is also available from which red magnitudes may be derived. 

2. Intensive spectroscopic studies by M. L. Humason, 
R. Minkowski and others have shown that super-nova 
spectra are totally different from the spectra of any other 
known stellar objects. 

3. The most luminous super-novae may be brighter than 
the nebulae in which they appear by as many as five 
magnitudes, as the case of the super-nova in I.C. 4182 
has shown. This brightest of all super-novae which is 
consequently the brightest stellar object ever observed, has, 
according to Baade, an absolute photographic magnitude 
close to M= — 16.6, corresponding to about four hundred 
million times the luminosity of the sun. The total radiation 
of energy from a bright super-nova between the wave- 
lengths A=6800A and A=3600A over a period of one 
year after its appearance amounts to about 10* ergs. 

4. The frequency of occurrence of super-novae has been 
determined as already mentioned. 

The most important conclusion which we may draw 
from these new observational results is, that the existence 
of two classes of temporary stars, super-novae and common 
novae, has been established beyond doubt. The new data 
have also furnished for the first time the necessary basis 
for a consistent statistical evaluation of the data available 
for previously observed super-novae. This evaluation has 
been carried out in detail by Dr. Baade and has resulted 
in significant results regarding the average absolute 
magnitudes of super-novae as well as the dispersion in 
magnitudes. 

California Institute of Technology, F. ZWICKY 

Pasadena, California, 

May 18, 1938. 


1'W. Baade and F. Zwicky, Phys. Rev. 45, 138 (1934) and 46, 67 
(1934); Proc. Nat. Acad. Sci. 20, 254 (1934) and 20, 259 (1934). See 
also F. J. M. Stratton, Handbuch der Astrophysik, Vol. 7 (1936), p. 671. 

2F. Zwicky, The Scientific Monthly 40, 461 (1935); Proc. Nat. 
Acad. 22, 457 (1936) and 22, 557 (1936); Publ. Astro. Soc. Pac. 48, 
191 (1936) and 49, 204 (1937). 





Angular Distribution of Slow Neutrons from a 
Paraffin Surface 


In experiments on the properties of slow neutrons and 
the induced radioactivities produced by slow neutrons it has 
been a common practice to use a Rn-Be neutron source 
embedded beneath the surface of a paraffin block and plane 
parallel layers of absorbers and detectors. The analysis of 
such experiments requires a knowledge of the angular 
distribution of neutrons coming from the surface. The 
expected distribution has been derived theoretically by 
Fermi! and found to be f(@)~ cos 6+ +/3 cos?@. The 
distribution has been studied experimentally and found to 
follow the Fermi distribution to within experimental errors. 
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Fic. 1. Experimental arrangement for observing the angular distribution 
of slow neutrons emerging trom a paraffin surtace 
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Fic. 2. Observed angular distribution of slow neutrons from a 
paraffin surface. The observations are shown with their standard 


errors and fit the Fermi distribution, /(@), to within the accuracy ot 
the experiments. 


Figure 1 shows the experimental arrangement. A Rn-Be 
source of neutrons, S, was placed 3 cm beneath the surface 
of the end of a paraffin cylinder of the kind described by 
Amaldi and Fermi.2? The neutrons emerging from the 
surface were detected in a BF; ionization chamber which 
was placed within cylindrical shields of BsC and Cd 
forming a collimator for slow neutrons. The ionization 
pulses due to the slow neutron disintegration of B' were 
observed with a pulse amplifier and counter. The variation 
of slow neutron intensity with angle was measured by 
rotating the BF; ionization chamber about a vertical axis 
in the plane of the paraffin surface. 

The cylinder was covered entirely with Cd (0.88 g 
Cd/cm?) with a hole 7 cm in diameter in the center of the 
end cover. This hole defined the area of the emitting 
surface, and was at the center of, and smaller than, the 
area ‘“‘seen’”’ by the detecting chamber. This insured having 
the same emitting surface effective at all angles. The BF; 
chamber subtended an angle of 6° 30’ from a point at the 
center of the emitting surface. The collimator tube ab- 
sorbed possible stray slow neutrons from the walls and 
floor. 

Counts were taken without and with Cd over the end of 
the collimator tube, the difference bet ween these runs being 
taken as the slow neutron count. Fig. 2 shows J/Jo plotted 
as a function of angle, 6. The observed points are shown 
with the fluctuations based on the standard errors of J and 
Ip. The solid lines give for purposes of comparison the 
functions cos 8, cos? 6, and f(@). The observed points tend to 
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follow the curve for f(@#). Runs taken with the neutron 


source 5 cm below the paraftin surface show the same trend 


Cornell Universit J. G. Horruas 
It! New York 
May 28. 193% M.S. Livincston 
I ber < Seow! 14, 134 1936 si 1] \ Bet 
Rev. Mod. ? 9, 132 (193 
kK. Amald 1 ke. Fer |’ Rev. 50, 901 rr 





Radial Lines in Laue Spot Photograms 


Wadlund! recently published in these columns under the 
above title some Laue photog: phs showing sharp Laue 
spots together with faint radial striae. He found these 
radial lines with ervstals of NaCl and KCI, but not with 
crystals of mica, Iceland spar, zine blend, fluorite, o1 
gypsum. Zacheriasen,® in discussing these lines, considered 
it unrezsonable to ascribe them to “strain,”’ and suggested 
that they were cross-grating spectra from a two-dimensional 
array of atoms. He designated boundary layers of the 
mosaic blocks of the imperfect NaCl and KCI crystals as 
the cross gratings responsible for the lines 

Neither author has mentioned that the cross-grating 
explanation of “asterism” in Laue photograms has already 
received much discussion, Linnik® proposed it as an expla- 
nation of patterns in heated mica, but Bragg! showed that 
the ordinary conception of diffraction from a_ three 
dimensional lattice was adequate to explain all features of 
Linnik’s patterns, if a small deviation of orientation about 
a mean orientation was assumed. He also offered experi 
mental proof of this, reproducing Linnik’s asterism by 
oscillating a normal mica crystal a few degrees in all direc 
tions from its mean position. Berg® has used the cross 
grating conception, but merely as a convenient method of 
analyzing patterns rather than as a fundamental mecha 
nism of diffraction. Bragg and Wirchner®: * pointed out an 
obstacle to the theory as applied to equi-axed crystallites 
by showing that in order for atoms in a row parallel to the 
incident beam to scatter in phase and act as one point of a 
cross grating (for x-rays, not electrons) the crystallite 
containing the row could not be thicker than one or two 
atom layers in this direction, and an equi-axed crystallite 
of this size would be quite incapable of producing sharp 
Cross grating spectra, 

Zachariasen’s argument? avoids the above objection, 
since it applies to diffraction from a surface layer of atoms 
of considerable extent laterally, but of very few atoms 
thickness, a laver in which the periodicity in the direction 
of the beam is presumably abnormal because it is on the 
surface of a mosaic block. But the cross-grating theory 
cannot be tested by Wadlund’s photographs, for these are 
fully explained by ordinary three-dimensional lattice 
diffraction, 

If one interprets Fig. 1 of Wadlund’s article in the con 
ventional way, it may be stated that the crystal consists of 
regions of diverging orientations; and that a majority of 
the regions scatter less than 0.1° from the mean, but a 
few deviate up to 8° in all directions. The horizontal and 
vertical streaks near the center of the photograph are (200 
O.32A upto 


A=0.63A where the effective intensity of the general radia 


reflections which form a spectrum from Ay i; 


tion and the amount of properly oriented material became 
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1 iron crystal. S 


otograph ot ipressed i 
and show a range rienta 


| 
10) and (200) plane 
ti 


small. The 45° streaks are (220) reflections which extend 
from A 0.32A out to a reflecting angle of 8.0° corr 
sponding to A4=0.55A, and are here terminated by the 
limit of scatter of orientations. It will be noted that both 
sets of streaks—as far as can be seen from the re proc 

have at ingential width and a shape cons stent with this 
amount of scatter 

Chere is a striking similarity between the streaks in the 
rocksalt pattern discussed above and the Laue photograph 
which accompanies this note, Fig. 1. The specimen was 
hyvdrogen-purified iron crystal that had been compressed in 
approximately the [100] direction to 1/20 of its original 
thickness, and set up in the camera with the compression 
axis parallel to the beam.* The position, length, intensit 
and sh ipe of every streak in the pattern are explicable Ss 
above, by assuming a preferred orientation of the cryst 
axes such that they filla solid angle around a mean position 
having two cube axes, one horizontal and one verti n 
the figure, and the third parallel to the axis of compression 
normal to the plane of the paper). The maximum deviation 
shown here is about 9°, The change in orientation was fol 
lowed as deformation progressed by means of 
photographs, and gradual changes in both the mean orient 
tion and the range of scatter were observed, The equilil 
rium state is shown in Fig. 1, 

lhe asterism in these cases is fully explaine e three 
dimensional lattice orientations of the tvpe to be expects 
in imperfect and deformed crystals and thus furnishes no 
evidence for a cross-grating eflect 


Metals Research Laboratory CHARLES S. BARRETI 
Carnegie Institute of Technology, 
Pittsburgh, Pennsylvania, 
May 27, 1938 


Wadlund, Phy Rev. 53, 843 (1938 

Zachariasen, Phy Rev. 53, 844 (1948 

Linnik, Nature 123, 604 (1929 
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Magnetic Anisotropy of Ferromagnetic Crystals in Weak 
Fields and Reversible Susceptibility 


The formula given for xy in the previous letter! is a 
particular case of the general formula for the reversible 
susceptibility which may be obtained by calculation of the 
change in the magnetization dJ caused by reversible 
displacement of the boundaries between domains when the 
field is changed by d/I, In iron crystals there are six and in 
nickel crystals eight directions of easy magnetization. In 
weak fields the spins of domains are parallel to these 
directions. Therefore, in weak fields there are in iron 
crystals six and in nickel crystals eight different magnetic 
phases. The volume concentrations of these phases will be 
denoted by m1, ++ +i, ++ +%1, ++ *i- ++ where 7; denotes the 
concentration of the phase whose spins are anti-parallel to 
those of the phase with concentration 7;. 

The general formula obtained for the reversible sus- 


ceptibility is of the following form: 


Kr=Krypt+kri'**, la) 
: vy . : - . 
Krp =Kro ae | (Mit Ms) (e+e) (hi?+hi2) 
+k ‘ F 
‘ —2(ns—7i)(me— MAM +, ~b) 


Krill at . (“ ay h Nits )h?---, (1c) 
; Nye+rN; J 

where /; and fy are the cosines of the angles between the 
directions of the spins of phases 7 and k and the direction 
of the field; «1: and.«-2 are constant coefficients which 
depend on the mechanical properties of the crystals. When 
I=0, i.e. when n;=7; formula (1b) turns into the formula 
for xy of the previous letter. Here, however, it should be 
noticed that the expression for x» given in the previous 
letter represents the initial susceptibility while (1) for 
ny=n; represents only the initial reversible susceptibility 
kro; Which is defined as the limiting value of the reversible 
susceptibility at J=0 and /7=0. Even though for ;=%; 
(1) coincides with the expression obtained previously for 
ky except that the multiplying factor is different, experiment 
shows that in some cases ko >k,o may take place. 

Formula (1) be used to calculate the reversible 
susceptibility of iron crystals in different directions. 

(1) The magnetic field /7\|[100]. In this case h,=1, 
h,=h;=0. IT=(n,;—ni)Is, whence 


may 


The magnetization: 


Wm—hi=jJ; Ne—he=n3—h;=0---, (2) 
where j=1/T;. 
the limiting case of the 


most advantageous 


(a) In 
orientation when all the spins are parallel or antiparallel to 


the direction [100], 
nmy+n,=1, 


Substituting (2) and (3) in (1) we obtain 


Ne=Ne2=nN3;=3=0---. 3) 


kre =O, ket =Krni(l—j?)--- (4) 
whence if «-1>0 it follows that 
Kr kro=1—j?-->, (5) 
(b) If in the initial state when J=0 the concentrations 
of the magnetic phases are equal, strict calculation of #;, is 
very difficult if we exclude the particular case when 
7T||[111]. In the general case »; may be determined as a 
function of 7 only from measurements of the magneto- 
striction. One may attempt to calculate «x, for the case of 
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equal concentrations in the initial state by using Heisen- 
berg’s' formula for #;, obtained by statistical methods, The 
applicability of statistical methods to the calculation of » 
has as yet no theoretical basis, nevertheless experiment 
shows that values of 1; calculated by these methods are in 
quite good agreement with those determined from magneto- 
striction measurements in the case of not very soft crystals 
and polycrystals. When the magnetization is parallel to the 
axis [100] statistical calculation gives: 


m+, = — 3 +35(14+3)?)); 


Ne+R2=Ny+hs = 5 — 3(14+3/? © 
Substituting (2) and (6) in (1) we obtain: 
Krp =(2 Dura L—4— OP +51 +37]; | 
kr S+3j2—4(143)2)) Ls «CO 
ee 3 142 1+3j)) | 


From (4) and (7) it follows that when //, [100], the ratio 
ky kro depends on the relationship between «,; and «;: 

(2) The magnetic field /7|/[110). 
analogous to the above gives: (a) in the limiting case when 
+h =Net+he=3; n3=n3=0 the suscepti- 


bilities 


This calculation 
reversible 
Kr = 4kr2(1— 2?) ; 

Krit = 2kri(1—2;?) ; see, (S) 
whence, independently of the relationship between xp, and Ky. 


we have: - 
Kr/Kro=1—2j?->-. (9) 


(b) When the phase concentrations are equal in the 
initial state 
Kea =(2/9)weal — 272+ 4(1—37)47; 
5 —6j?—4(1 — 3j?)3 


Kl = 3kri*— 2—(1- 3 j7)4 
Zz— — §7*) 


(10) 


as in the analogous case for Z||[100] the ratio «,/Kyo 
depends on the relationship between «,; and «;2. It should, 
however, be noticed that when «;:<x,2 the expression for 
k,/kro Obtained from formula (10) is practically very little 
different from expression (9), 

(3) J7\/[111]. If in this case the concentrations are equal 
in the initial state calculation gives: 


(11) 


Hence again independently of the relationship between x, and 


key = (2/9)Kra(1—3j2); Kt = der(1—3;2). 


kro We obtain: 7 
- Kr/Kro= 1 — 3)”. (12) 


There will be divergences from (12) if the section of the 
specimen perpendicular to the direction of the magnetic 
field deviates considerably from the circular form and 
if the demagnetizing factor becomes either less or greater 
along one of the axes of easy magnetization than along the 
others. The expression for x, may be obtained in an 
analogous manner for crystals of nickel. 

Physical Institute of the State University, E. KONDORSKY 


Moscow, U. De De R 
April 8, 1938. 


1 E. Kondorsky, Phys. Rev. 53, 319 (1938), 

2 The deduction of formulae (1) will shortly be published in the 
Comptes Rendus Acad. Sci. U.S. S. R. Vol. 19, 4. 

>It should be noticed that if the concentrations of the magnetic 
phases are unequal in the initial state, the directions of easy magneti- 
zation can no longer be considered as physically equal, and Eqs. (1) 
become inapplicable. However, in the limiting cases of unequal concen- 
trations considered below, the concentrations of some of the phases 
are zero and those of the remainder are equal between themselves in 
the initial state. Therefore formulae (1) remain applicable to the 
limiting cases in question. 

4W. lleisenberg, Zeits. f. Physik 62, 287 (1931). 
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G. Singer—916(A) 

Radiography of thin microtome sections, E. Dershem 
202(A) 

Short exposure, high intensity x-rays, K. H. Kingdon, 
H. E. Tanis, Jr.—128 

X-ray generator for biological experiments, G. Failla 
916(A) 

Brownian movement 
Photographic measurements, C. Hire—205(A); 325(A) 


Collision diameter (see Atomic and molecular radii) 
Conductivity, electrical (see Electrical conductivity; Dis- 
charge of electricity in gases) 
Constants, physical, standards 
Discrepancy among experimental measurements, R. A. 
Beth—681 (A) 
Variability of constants, F. Zwicky—315(L) 
Contact potentials 
Of solid films on metals, W. D. Harkins, R. W. Mattoon 
—911(L) 
Cosmic radiation 
Absorption of penetrating component, W. M. Nielsen, 
K. Z. Morgan—915(A); V. C. Wilson—908(L) 
Analysis including hard component, L. W. Nordheim 
329(A); 694 
Artificial radioactivity, W. E. Ramsey, C. G. Mont- 
gomery, D. D. Montgomery—196(L) 
Bursts in equatorial zone, S. A. Korff—914(A) 
Change of mass, Ig. Tamm—1016(L) 
Cosmic-ray telescope, W. H. Pickering—335(A) 
East-west effect in India, H. V. Neher—328(A) 
Effect of earth’s magnetism, T. H. Johnson—329(A) 
Electrical conductivity in stratosphere, O. H. Gish, 
K. L. Sherman—434(L) 
Energy and intensity of primary, T. H. Johnson—499; 
682(A) 
Energy distribution, I. S. Bowen, R. A. Millikan, 
H. V. Neher—328(A); 915(A) 
Energy flux of corpuscular radiation, C. G. Mont- 
gomery, D. D. Montgomery—193; 329(A) 
Galactic rotation, I. A. Getting—914(A) 
Identification of primary rays, E. J. Schremp—915(A) 
Influence of earth’s field, W. Kolhérster—768(L) 
Influence of solar magnetic field, P. S. Epstein—862 
In stratosphere, L. F. Curtiss, A. V. Astin, L. L. Stock- 
mann, B. W. Brown, S. A. Korff—23; 330(A) 
Intensity and magnetic storms, S. E. Forbush—914(A) 
Intensities at great depths, V. C. Wilson—204(A); 337 
Ionization near a lead block, E. G. Schneider—615 
Large bursts, function of Pb shielding, W. P. Jesse, 
R. L. Doan—203(A); 691 
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Cosmic radiation (continued) 

Large bursts, theory of, C. G. Montgomery, D. D. 
Montgomery—682(A); 955 

Latitude effect, high altitudes, S. A. Korff, L. F. Curtiss, 
A. V. Astin—14 

Latitude effect, showers, H. V. Neher, W. H. Pickering— 
111; 329(A) 

Latitude effect of soft rays, B. Gross—607(L) 

Nature of penetrating component, R. M. Langer— 
494(L) 

Nature of radiation, I. S. Bowen, R. A. Millikan, H. V. 
Neher—214(A); 217; 855 

Nuclear disintegrations, R. B. Brode, M. A. Starr—3; 
214(A) 

Penetrating particles, S. H. Neddermeyer—102(L) 

Range of particles, A. Bramley—915(A) 

Scattering of particles, E. J. Williams—433(L) 

Showers at great depths, W. H. Pickering—215(A) 

Showers in lead, M. A. Starr—6 

Sidereal diurnal variation, S. E. Forbush—682(A) 

Specific ionization and mass of particles, D. R. Corson, 
R. B. Brode—215(A); 773 

Stopping power of heavy elements, J. H. Bartlett, Jr.— 
915(A) 

Theory of showers, D. Iwanenko, A. Sokolow—910(L) 

Transition effects in atmosphere, H. Snyder—960 

Cosmic rays 

Penetrating component, J. 
848(L) 

Crystal structure 

Form of acetylsalicylic acid in certain sugars, S. S. 
Sidhu—918(A) 

Of KH.(H;0).B;Ow, W. H. Zachariasen—202(A) 

Of N(CHs),Is, R. C. L. Mooney—851(A) 

Of NaHCOs, W. H. Zachariasen—917 (A) 

Space group of pyrrhotite, V. Hicks, S. S. Sidhu— 
207(A) 

Of N(CHs),4I Cle, R. C. L. Mooney—918(A) 

Of uranium, C. W. Jacob, B. E. Warren—108(A) 

Crystalline state 

Alkali halide-thallium phosphors, F. Seitz—685(A) 

Anisotropy of crystals, E. Kondorsky—1022(L) 

Anisotropy of vibrations in a Zn-Cu alloy, R. A. Howard 
—324(A); 966 

Asymmetry of atoms in Ag-Cd alloy, E. O. Wollan— 
203(A) 

Barium-copper-stearate films, properties of, C. Holley— 
534; I. Fankuchen—909(L); J. J. Bikerman, J. H. 
Schulman—909(L) 

Calculation of electronic energy bands, G. Wannier— 
671(L) 

Characteristic temperature of MgO, H. S. 
E. O. Wollan—332(A); 972 

Characteristic temperature of Zn crystals, G. E. M. 
Jauncey—203(A) 

Cold work and internal friction, C. Zener—582 

Critical cooling rate of iron-carbon alloys, T. G. Digges 
—925(A) 

Crystalline transitions and dielectric constant, H. v. R. 
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Zener—90; 100; 582; 
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Lattice vibrations in polar crystals, R. H. Lyddane, 
K. F. Herzfeld—948(A) 

Order versus disorder, lattice enlargement, F. C. Blake- 
333(A) 
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W. H. Zachariasen, 844(L); C. S. Barrett—1021(L) 


Rate of recrystallization, L. A. Young—686(A) 

Recalescence in uranium, W. L. Hole, R. Wright, 
H. B. Wahlin—768(L) 

Rotational transition in solid methane, J. A. Wheeler, 
C. V. Cannon—684(A) 

Segregation of Po in Bi crystals, A. B. Focke—324(A) 

Solubility at grain boundaries of solid solutions, G. R. 
Dean—202(A) 

Supercooling and freezing of water, N. E. 
925(A) 

Superstructure and magnetism of pyrrhotite, S. S. Sidhu, 
V. Hicks—207(A) 

Transition energies of polymorphous molecular crystals, 
M. Bruch-Willstatter, H. Sponer—324(A) 

Transitions in alkali-halides, R. B. Jacobs—930(A) 

Vitreous state in gelatin gels, B. J. Luyet—323(A) 
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Permanent volume polarization, R. W. Goranson, W. A. 
Zisman—668(L) 

Polar moment of cellosolves, W. H. Byers, A. A. Bless— 
852(A) 

Susceptibility of Rochelle salt, D. Bancroft—587 
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Clean-up effect in low pressure Hg vapor discharges, 
C. Kent y—936(A) 
Cu arc in air, A. B. White—935(A) 
Corona discharge, G. W. Trichel—214(A) 
Discharges under nearly complete ionization, F. L. 
Mohler—683(A); 923(A) 
Force on cathode of Cu arc, R. M. Robertson—578 
Glow-arc transition in Hg vapor, F. A. Maxfield 
204(A) 
Impulse breakdown in long tubes, J. R. Dietrich, L. B. 
Snoddy, J. W. Beams—923(A) 
Mechanism of positive column in Hg, E. Q. Adams, 
B. T. Barnes—556 
New mechanism in excited neon, E. W. Pike—920(A) 
Oxidation in arc cathode, C. G. Suits, J. P. Hocker 
670(L); J. D. Cobine—911(L) 
Positive-point to plane discharge, A. F. Kip—210(A) 
Recombination in pure Oz, M. E. Gardner—75 
Sparking potential, F. Ehrenkranz—210(A) 
Sparking voltage of sphere gap, effect of humidity, 
A. B. Lewis—850(A) 
States in high pressure discharge, C. G. Suits—609(L) 
Vaporization of Hg from anchored cathode spot, L. 
Tonks—936(A) 
Discharge of electricity in high vacuum (see Thermionic 
emission ; High voltage tubes and machines) 
Disintegration of nucleus (see also Radioactivity) 
Of Al by Het, H. L. Schultz—929(A) 
Of C; y-radiation, M. E. Rose—844(L) 
Angular distribution of disintegration products, R. D. 
Myers—325(A) 
Angular distribution of H! from H?+H?, R. O. Haxby, 
J. S. Allen, J. H. Williams—921(A) 
Of Be by He?, H. L. Schultz—681(A) 
Of B” and Li® by neutrons, M. S. Livingston, J. G. 
Hoffman—227 
B"+H! reactions, J. R. Oppenheimer, R. Serber—636 
Cloud chamber studies, L. W. Alvarez, W. M. Brobeck 
—213(A) 
Collision relations, A. E. Ruark, C. Jones—853(A) 
Of H? by He’, H. L. Schultz—622 
Deuteron reactions at high energy, L. W. Alvarez— 
326(A) 
Deviations from evaporation model of reactions, H. A. 
Bethe—675(A) 
Efficiency of Be+H!, L. S. Skaggs, G. T. Hatch, S. K. 
Allison—946(A) 
Energy of neutrons from C!+H?, T. W. Bonner— 
496(L) 


Energy-range relations, F. T. Rogers, Jr., M. M. Rogers 
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Energy spectra of particles, S. K. Allison, N. M. Smith, 
Jr.—946(A) 

Formation of Be’? and K electron capture, R. B. Roberts, 
N. P. Heydenburg—929(A) 

Formation of excited He*, T. W. Bonner—711 

Gamma-radiation from B+H!, W. A. Fowler, E. R. 
Gaerttner, C. C. Lauritsen—212(A); 628 

Loss of neutron and alpha-particle in Fe, relative proba- 
bility, J. M. Cork, B. R. Curtis—946(A) 
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Neutron-deuteron impacts, W. E. Lamb, Jr., R. Serber 
215(A) 

Neutron interaction with nuclei and neutron energy, 
P. N. Powers, H. H. Goldsmith, H. G. Bever, J. R 
Dunning—947 (A) 

Neutrons from A, Cl, Al and under Het bombardment, 
E. Pollard, H. L. Schultz, G. Brubaker—351 

Neutrons from B, H. Staub, W. E. Stephens—212(A) 

Neutrons from Li+H?, W. E. Stephens—212(A); 223; 
M. L. Pool—707 

Of N by neutrons, F. N. D. Kurie, M. Kamen—212(A) 

Nuclear excitation functions, M. S. Plesset, D. H. 
Ewing—943(A) 

Oppenheimer-Phillips process, H. A. Bethe—39; 205(A) 

Photodisintegration of deuteron, J. R. Richardson, 
L. Emo—234 

Production of N', R. B. Roberts, N. P. Heydenburg 
374 

Production of pairs, W. E. Lamb, Jr.—211(A) 

Scale for predicting nuclear transformations, V. Kara- 
petoff—675(A) 

Of Sc by Het, E. Pollard—679(A) 

Of Sc by neutrons, J. M. Cork, R. L. Thornton—866 

Thermonuclear reactions, G. Gamow, E. Teller—608(L) 

Of Ti by He*, W. L. Davidson, Jr., E. Pollard—-929(A) 

Width of resonance in B+H)!, J. S. Allen, R. O. Haxby, 
J. H. Williams—325(A) 


Dispersion 


Verification of the quantum mechanical theory, W. A. 
Bruce, G. E. M. Jauncey—332(A); W. A. Bruce—802 


Dissociation 


Of CCI,F2 and in CCl,, R. F. Baker, J. T. Tate—683(A) 
Of ethane by electron impact, J. A. Hipple, Jr.—530 
Of normal and isobutane, R. F. Baker, J. T. Tate 


944(A) 


Elasticity 


Determination of Young’s modulus, R. R. Jacobs, 
D. Bancroft—687 (A) 

Filtration of compressional waves in stratified medium, 
R. B. Lindsay—944(A) 

Fixed path acoustic interferometer, I. F. 
J. C. Hubbard—944(A) 

Internal friction in solids, C. Zener—90: 100; 582: 
686(A); 1010 


Zartman, 


-Thermoelastics and structure of rubber, F. E. Dart, 


E. Guth—327(A) 

Vibration modes of low decrement, K. S. Van Dyke 
945(A) 

Vibration of rotating propeller blades, W. Ramberg, 
S. Levy—945(A) 


Electrical conductivity and resistance 


Deviations from Ohm's law, E. Guth, J. Mayerhdéfer 
205(A) 

Drift of ions on filaments, R. P. Johnson—766(L) 

Effect of tension; Sn crystals, M. Allen—108(A) 

Electrical breakdown in alkali halides, R. J. Seeger, 
E. Teller—923(A) 

At frequencies of the order 300 mC/second, R. A. Chip- 
man—672(L) 
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Electrical conductivity and resistance (continued) 

Of metallic films condensed at low temperatures, H. R. 

Mitchell—250 
Electric moment (see Dielectrics and dielectric properties) 
Electrolytic cells 

Anodic properties of Al cell, W. B. Pietenpol, T. E. 
Devaney—325(A) 

Electrode capacitance of conductivity cells, S. F. Acree 
—923(A) 

Solubility of glass, behavior of glass electrodes, D. Hub- 
bard—935(A) 

Electromagnetic theory 

Electromagnetic and gravitational radiation, L. Infeld— 
836 

Light signals around a closed path, H. E. Ives—939(A) 

Mutual inductance between two coaxial helical wires, 
C. Snow—934(A) 

Oscillations of a prolate spheroid, L. Page, N. I. Adams, 
Jr.—819 

Electron diffraction (see also Electrons, scattering of) 

Molecules in unimolecular and multimolecular layers, 

L. H. Germer, K. H. Storks—322(A) 
Electrons (see also Positive electron and e/m) 

Energy loss in Pb, A. J. Ruhlig, H. R. Crane—618; 683(A) 

Production of heavy electrons, L. W. Nordheim, E. 
Teller—937 (A) 

Electrons, secondary 
From Na films on Ta, P. L. Copeland—328(A) 
From thorium-coated tungsten, E. A. Coomes—936(A) 
Electro-optical effects (see also Photoelectric effect and 
properties) 
Optical activity of NiSO., E. B. Nelson—849(A) 
E/m 

Analysis of classical experiments, C. T. Zahn, A. H. 
Spees—511 

Of B-rays from RaE, C. T. Zahn, A. H. Spees—365 

Change of rest mass during collisions, F. Zwicky— 
611(L); Ig. Tamm—1016(L) 

Determination for 8-particles, C. T. Zahn, A. H. Spees— 
357 

Energy losses of Pauli-Weisskopf particles, J. F. Carl- 
son—330(A) 

Heavy beta-rays, A. H. Compton—431(L); H. R. Crane 
—317(L);G. E. M. Jauncey—197(L); 265(L); 319(L); 
669(L); A. Ruark, C. C. Jones—264(L); 454; 496(L); 
679(A); C. T. Zahn—431(L); F. L. Arnot—1018(L)- 

Method for precision determination, A. E. Shaw—208(A) 

Nature of heavy electron, I. M. Freeman—606(L) 

Particle of intermediate mass, A. J. Ruhlig, H. R. Crane 
—266(L) 

From refraction of x-rays, J. A. Bearden—917(A) 

Energy states of atoms (see also Spectra, atomic) 

Self-consistent field for Fe, M. F. Manning, L. Gold- 

berg—662 
Energy states of molecules (see also Molecular structure 
and constants; Spectra, molecular) 

Energy levels of the asymmetrical rotator, E. E. Witmer 
—331(A) 

Energy states of nucleus (see also Nucleus) 
Of alpha-particles, excited states, H. Margenau—198(L) 
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Barytron theory of nuclear forces, H. A. Bethe —938(A) 

Binding energy of H?, H. A. Bethe—313 

Binding energy of O'*, G. Horvay, E. Feenberg—938(A) 

Compound nucleus, H. A. Bethe, E. J. Konopinski— 
688(A) 

Dynaton theory of nuclear forces, R. Serber—211(A) 

Energy of heavy nuclei, A. J. Dempster—869 

Excitation of nuclei, V. Weisskopf—1018(L) 

Excited He? in H?+H?, T. W. Bonner—928(A) 

Fine structure of ground level of Li’, G. Breit, J. R. 
Stehn—459; 684(A) 

Ground state of deuteron, R. G. Sachs, M. Goeppert- 
Mayer—991 

Hartree approximation, V. F. Weisskopf—680(A) 

Heavy particle interactions, D. R. Inglis—686(A) 

Of H* and He‘, S. S. Share—875 

Of Li*, K. G. Carroll, H. Margenau—921 (A) 

Multipole radiation probabilities, M. Phillips—681(A) 

Rotation of atomic nuclei, E. Teller, J. A. Wheeler— 


684(A); 778 

Saturation of nuclear forces, C. Critchfield, E. Teller— 
684(A); 812 

Saturation requirements for forces, G. Breit, E. Wigner 
—998 


Spacing of nuclear levels, J. Bardeen, E. Feenberg— 
938(A) 

Superposition of nuclear forces, H. Primakoff—938(A) 

Variational theory of alpha-particle, W. A. Tyrrell, Jr., 
H. Margenau—939(A) 

Errata 

Ferromagnetic anisotropy of nickel-iron crystals at vari- 

ous temperatures, J. D. Kleis—848(L) 


Geophysics 
Age of earth from isotope ratios, A. Bramley—323(A) 
Determination of geologic ages, C. Goodman, R. D. 
Evans—916(A) 
Substandards of feeble radioactivity, R. 
C. Goodman—916(A) 
Gravitation 
Electromagnetic and gravitational radiation, L. Infeld— 
836 


D. Evans, 


Hall effect 
Nature of the effect, E. Weber—935(A) 
Heat of combustion 
Of tetramethylmethane, J. W. Knowlton, F. D. Rossini 
—946(A) 
High voltage tubes and machines 
Accelerator using a quarter wave-length, P. L. Hartman, 
L. P. Smith—677 (A) 
Compact 750 kv Van de Graaff generator, J. B. Fisk, 
I. A. Getting—916(A) 
Cyclotron magnet-current stabilizer, T. Perry—943(A) 
Energy variations in a cyclotron ion beam, G. K. Green, 
P. G. Kruger—927(A) 
Focusing and energy of ions in cyclotron, M. E. Rose, 
H. A. Bethe—206(A); M. E. Rose—392; 675(A); 715: 
R. R. Wilson—213(A); 408 
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Fringing flux in large magnets, J. D. Howe, I. Waler- 
stein—334(A) 

Generator for production of neutrons, T. R. Folsom— 
687 (A) 

High air pressure electrostatic generator, D. B. Parkin- 
son, R. G. Herb, E. J. Bernet, J. L. McKibben—642 
High frequency systems for cyclotron, J. R. Dunning, 

H. L. Anderson—334(A) 
Ion optics of an accelerating tube, J. S. Allen, J. H. 
Williams—925(A) 
Minnesota pressure electrostatic generator, J. T. Tate, 
L. H. Rumbaugh, J. H. Williams—928(A) 
X-ray generator for biological experiments, G. Failla— 
916(A) 
Hydrodynamics 
Aerodynamics of golf club, S. J. Crooker—927(A) 
Theory of isotropic turbulence, H. P. Robertson— 
939(A) 
Hyperfine structure (see also Nuclear moments and spin) 
Of B, Y, Rh and Pd, D. T. Williams, L. P. Granath— 
685(A) 
Interval rule, D. H. Tomboulian, R. F. Bacher—678(A) 
Of Mnl, R. A. Fisher, E. R. Peck—685(A) 
Of spectral lines from atomic beams, K. W. Meissner— 
677(A) 


Instruments (see Methods and instruments) 
Ionization by particles (see Ionization potentials) 
Ionization potentials 
Of CCI.F, and in CCl,, R. F. Baker, J. T. Tate—683(A) 
And dissociation of ethane, J. A. Hipple, Jr.—530 
Of Hg, W. B. Nottingham—328(A) 
Of normal and isobutane, R. F. Baker, J. T. Tate—944(A) 
Use of isotopes in study of ionization processes, J. M. 
Delfosse, J. A. Hipple, Jr.—683(A) 
Ionosphere 
Ionization in E region, J. A. Pierce—942(A) 
Ionization and recombination, N. E. Bradbury—210(A) 
Lorentz polarization correction, H. G. Booker, L. V. 
Berkner—924(A) 
Photoelectric ionization, E. O. Hulburt—344 
Variation in F, ionization, E.O. Hulburt—670(L) 
Ions, mobility 
Liquids electrified by bubbling, S. Chapman—211(A) 
Isotopes 
Abundances of Pb isotopes, A. O. Nier—680(A) 
Of Ca, Ti, S and A, A. O. Nier—282 
Of C from CN spectrum, F. A. Jenkins, D. E. Wool- 
dridge—137 
Concentration of N“, H. G. Thode, J. E. Gorham, 
H. C. Urey—920(A) 
Of Gd, Dy, Er and Yb, A. J. Dempster—727 
Of heavy elements, A. J. Dempster—64 
Of uranium, A. O. Nier—922(A) 
Range of occurrence, H. Brown—846(L) 
Stable isotopes and atomic numbers, G. E. Gibson, 
H. Brown—846(L) 
Weight of He by direct comparison with O, K. T. Bain- 
bridge—922(A) 
Of Xe, Kr, As, Cs, I and Be, A. O. Nier—109(A) 
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Kinetic theory of gases 
Air-lubricated journal bearings, M. D. Hersey 
Motion of ions in gases, F. B. Pidduck—197(L) 


926(A) 


Persistence of molecular vibration in collisions, R. J. 
Dwyer—324(A) 


Liquids (see also Hydrodynamics) 

Atomic distribution in liquid K, C. D. Thomas, N. S. 
Gingrich—943 (A) 

Atomic distribution in liquid Na, F. H. Trimble, N. S. 
Gingrich—203(A); 278 

Atomic distribution in liquid white P, N. S. Gingrich, 
C. D. Thomas—943(A) 

Film lubrication of the journal bearings, S. M. Skinner— 
927(A) 


Supercooling and freezing of water, N. E. Dorsey 
925(A) 

Variation of ionic volumes, G. W. Stewart—323(A) 

Vitreous state in gelatin gels, B. J. Luyet—323(A) 


Luminescence 
Of pure crystals, D. H. Kabakjian—919(A) 
Spectra of cathodoluminescent materials, H. W. Lever- 
enz—919(A) 


Magnetic fields 
Demagnetization coefficient of oblate spheroids, L. P. 
Tarasov—105(L) 
Production of uniform fields, R. J. Stephenson, M. 
Ference, A. E. Shaw—208(A) 
Magnetic properties 
Anisotropy in ferromagnetic crystals, E. Kondorsky— 
319(L) 
Anisotropy of crystals, E. Kondorsky—1022(L) 
Diamagnetism of small regions, J. C. Slater—208(A) 
Ferromagnetic impurities, detection of, F. W. Constant, 
J. M. Formwalt—432(L); 852(A) 
Hysteresis testing, C. T. Razor—851(A) 
Interaction in Heusler alloy, L. W. McKeehan—199(L) 
Interaction in pyrrhotite and; in magnetite, L. W. 
McKeehan—307 
Interaction in strained crystals, L. W. McKeehan—301 
Magnetic viscosity in iron, C. W. Heaps—-334(A) 
Magnetostriction of polycrystalline material, W. F. 
Brown, Jr.—482 
Photomicrographs of magnetic colloid patterns, W. C. 
Elmore—933(A) 
Of pyrrhotite, S. S. Sidhu, V. Hicks—207(A) 
Structure of Co, W. C. Elmore—757 
Surface magnetization of Co, W. C. Elmore—108(A) 
Susceptibility of lead oxide, M. Petersen—933(A) 
Magneto-optical effects 
Magnetic double refraction; acetone, esters, J. B. 
Hawkes—84 
Verdet constants of CS, glycerine, isopropyl and amyl 
alcohols, G. C. Comstock—849(A) 
Mass defects (see Isotopes) 
Mass spectrography 
Ion source (Li), Hoff Lu—845(L) 
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Mechanics, quantum—general 

Lower bounds in quantum mechanics, A. F. Stevenson— 
199(L) 

Neutral particle, Theory of, W. H. Furry—679(A) 

Quadrupole-quadrupole interatomic forces, J. K. Knipp 
—734 

Theory of neutral particle, W. H. Furry—679(A) 

Mechanics, quantum—of solid bodies 
Calculation of electronic energy bands, G. Wannier— 

671(L) 
Compressibility of alkali metals, J. Bardeen—683(A) 
Meteorology 
Atmospheric absorption by water, W. M. Elsasser— 
768(L) 

Density and temperature of the atmosphere to 60 km, 
E. O. Hulburt—948(A) 

Determination of atmospheric ozone, B. O’Brien, H. S. 
Stewart, Jr.—949(A) 

Factors affecting pyrheliometers, L. F. Miller—948(A) 
Methods and instruments 
Accelerator using a quarter wave-length transmission 
line, P. L. Hartman, L. P. Smith—677(A) 

Atomic electron velocities; measuring of, A. L. Hughes, 
M. M. Mann, Jr.—50 

Barometer, for small volumes, W. Hurst—850(A) 

Capacity controlled magnetic suspension, C. S. Smith— 
933(A) 

Capture cross section, A. H. Spees, W. F. Colby, S. 
Goudsmit—326(A) 

Centrifuge, tubular, J. W. Beams—850(A) 

Circuit for G-M counters, E. W. Yetter—612(L) 
* Cloud chamber technique, F. N. D. Kurie—215(A) 
Coaxial tube counter, S. C. Brown, R. D. Evans—914(A) 
Coincidences in counter circuits, C. Eckart, F. R. Shonka 
—752 

Collimation of slow neutrons, J. G. Hoffman, M. S. 
Livingston—676(A) 

Collision relations, A. E. Ruark, C. Jones—853(A) 

Compact 750 kv Van de Graaff generator, J. B. Fisk, 
I, A. Getting—916(A) 

Counting losses of single scale recorders, H. Lifschutz, 
O. S. Duffendack—941 (A) 

Cyclotron magnet-current stabilizer, T. Perry—943(A) 

Deposition of Ra Cu on rotating cathode, J. Steigman— 
771(L) 

Detection of cocain-HCl in cerebrospinal fluid, L. A. 
Strait, R. B. Aird—213(A) 

Determination of e/m, A. E. Shaw—208(A) 

Determination of geologic ages, C. Goodman, R. D. 
Evans—916(A) 

Determination of specific charge of beta-particles, C. T. 
Zahn, A. H. Spees—357 

Displacement polarimeter, N. Underwood—853(A) 

Effective mutual inductance, alternating current method, 
H. L. Curtis, A. V. Astin, C. M. Sparks—934(A) 

Electron velocity analyzer, M. E. Bell, E. A. Coomes, 
W. B. Nottingham—328(A) 

Electronic stabilizer performance, F. V. Hunt, R. W. 
Hickman—913(A) 





SUBJECT INDEX 


Energy variations in a cyclotron ion beam, G. K. Green, 
P. G. Kruger—927(A) 

Errors in count from G-M tube, E. Rodgers—850(A) 

Factors affecting pyrheliometers, L. F. Miller—948(A) 

Ferromagnetic impurities, detection of, F. W. Constant, 
J. M. Formwalt—432(L); 933(A) 

Filament surfaces, effects of a.c. and d.c., R. P. Johnson 
—936(A) 

Fixed path acoustic interferometer, I. F. Zartman, J. C. 
Hubbard—944(A) 

Focusing and energy of ions in cyclotron, M. E. Rose, 
H. A. Bethe—206(A); M. E. Rose—392; 675(A); 715; 
R. R. Wilson—213(A); 408 

Focusing x-rays by rocksalt crystal, R. M. Bozorth, F. E. 
Haworth—538; I. Fankuchen—910(L) 

Fringing flux in large magnets, J. D. Howe, I. Walerstein 


—334(A) 

Gas free surfaces, H. E. Farnsworth, R. P. Winch 
935(A) 

Geiger counter, T. G. How, K. Lark-Horovitz—-334(A) 


Geiger counter circuit, I. A. Getting—103(L); H. V. 
Neher, W. H. Pickering—316(L); T. H. Johnson 
914(A) 

Geiger counters, use of, H. V. Neher, W. H. Pickering 
214(A) 

G-M tube design, J. L. Lawson, A. W. Tyler—605(L) 

Generator for production of neutrons, T. R. Folsom 
687 (A) 

Grating, mounting for, H. Sponer, W. Hurst—851(A) 

High frequency systems for cyclotron, J. R. Dunning, 
H. L. Anderson—334(A) 

High pressure capillary lamps, C. Bol—214(A) 

High speed electrostatic motor, C. S. Smith—688(A) 

Hollow cathode discharge source for Zeeman effect, 
R. A. Fisher, A. S. Fry, J. R. Platt—934(A) 

Infra-red wave-lengths, for calibration, P. E. 
E. K. Plyler—851(A) 

Integrated intensity from x-ray powder photographs, 
E. O. Wollan—208(A) 

Intensifying screens in x-ray diffraction work, C. 
Gamertsfelder, N. S. Gingrich—332(A) 

Intensities of spectral lines, D. Foster, H. W. Highriter— 
673(L) 

Ion source (Li), Hoff Lu—845(L) 

Ion sources using focused beams, L. P. Smith, G. W. 
Scott, Jr.—677(A) 

Ionization current meter, G. Failla—687(A) 

Light source for cloud chamber, W. M. Powell—327(A) 

Lissajous figures, F. H. Mitchell—850(A) 

Magnetic fields and susceptibilities in solenoids, A. R. 
Kaufman—933(A) 

Magnetic laboratory at M.I.T., F. Bitter—933(A) 

Mass spectrometer, new, W. Bleakney, J. A. Hipple, Jr. 
—521 

Measurement of the speed of a rotating shaft, F. Wenner, 
I. L. Cooter—927(A) 

Measurement of super-voltage radiation, L. S. Taylor, 
G. Singer—916(A) 

Measuring nuclear magnetic moment, I. I. Rabi, J. R. 
Zacharias, S. Millman, P. Kusch—318(L) 


Shearin, 
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Mechanical harmonic synthesizer-analyzer, S. L. Brown 
333(A) 

Multivibrator counter circuit, A. Ruark—316(L) 

Multivibrators from 200,000 to 1 cycle second, W. D. 
George—924(A) 

New ether-drift interferometer, W. B. Cartmel—108(A) 

New mass spectrometer, W. Bleakney, J. A. Hipple, Jr. 
—521 

Nonideal Lloyd's mirror, F. W. Doermann—420 

Nuclear recoil in 8-ray emission, J. Halpern, H. R. Crane 
—676(A) 

Numerical integration, M. F. Manning, J. Millman— 
673(L) 

Numerical solution of Laplace’s equation, G. H. Short- 
ley, R. Weller—207(A) 

With oscillograph, C. T. Razor—851(A) 

Permanent electrets, W. M. Good—323(A) 

Photoelectric counters with Cs cathodes, G. L. Locher— 
333(A) 

Production of uniform magnetic fields, R. J. Stephenson, 
M. Ference, A. E. Shaw—208(A) 

Radiography of thin microtome sections, E. Dershem— 
202(A) 

Reduction of barium oxide, J. P. Blewett—935(A) 

Scale of eight circuit, H. Lifschutz, J. L. Lawson—204(A) 

Scale for predicting nuclear transformations, V. Karape- 
toff—675(A) 

Sensitivity of balanced space charge method, A. Rostagni 
—729; R. N. Varney—732 

Short exposure; high intensity x-rays, K. H. Kingdon, 
H. E. Tanis, Jr.—128 

Single crystal diffraction patterns, F. Miller, Jr., L. G. 
Parratt—200(L) 

Solar spectrum photography, H. S. Stewart, Jr., B. 
O’Brien—949(A) 

Surface viscosimeter, W. D. Harkins, J. G. Kirkwood— 
322(A) 

Temperature control for absorption cells, D. Williams, 
R. Taschek—851(A) 

Theory of scaling circuit, L. Alaoglu, N. M. Smith, Jr.— 
832 

Thermionic cathode, new type, A. W. Hull—936(A) 

Three-component accelerometer, J. E. Shrader—333(A) 

Ultramicrometer, D. L. Arenberg, P. M. Roope—687 (A) 

Unusual demonstrations with piezoelectric resonators, 
K. S. Van Dyke—686(A) 

Vacuum furnace for production of single crystals, F. C. 
Nix—917(A) 

Vacuum tube frequency meter, H. Lifschutz—950(A) 

Vacuum-type centrifuge, S. A. Black, J. W. Beams, 
L. B. Snoddy—924(A); J. W. Beams, C. Skarstrom— 
327(A) 

Velocity of sound in air, R. C. Colwell, A. W. Friend, 
D. A. McGraw—686(A) 

Viewing and measuring instrument for cloud chamber, 
C. Jones, A. Ruark—677(A) 

Viscometer for lubricating oils, B. W. Thomas, R. B. 
Dow, W. R. Ham—926(A) 

Voltage multiplying circuit, S. K. Allison, G. T. Hatch, 
L.S. Skaggs, N. M. Smith, Jr.—204(A) 


Wilson cloud chamber control, I. I. Hopkins, R. T. 
Dickerson, W. M. Nielsen—850(A) 

X-ray camera for high temperatures, C. Nusbaum 
918(A) 

X-ray generator for biological experiments, G. Failla 
916(A) 

Young’s modulus, determination of, R. R. Jacobs, 
D. Bancroft—687 (A) 


Mobility of ions (see Ions, mobility) 
Molecular structure and constants (see also Spectra, 


molecular and Raman spectra) 

Deuterium bonds, W. Gordy—851(A) 

Potential barrier in ethylene chloride, J. Y. Beach, K. J. 
Palmer—685(A) 

Potential functions of the methyl halides, D. M. Denni- 
son, Z. I. Slawsky—930(A) 


Neutrino (see also Radioactivity) 


Evidence for neutrino, H. R. Crane, J. Halpern—789 
Heavy neutrino, F. L. Arnot—1018(L) 
Heavy particles; existence of, G. E. M. Jauncey—106(L) 


Neutrons (see also Scattering of electrons, neutrons and 


ions) 

Absorption in boron, H. C. Torrey—266(L) 

Absorption by Cl, E. O. Salant, W. J. Horvath—943(A) 

Albedo, slow neutrons, E. Fermi, E. Amaldi, G. C. Wick 
—493(L) 

Angular distribution, M. L. Pool—334(A); J. G. Hoff- 
man, M. S. Livingston—1020(L) 

From B+H?2, H. Staub, W. E. Stephens—212(A) 

From C"®+H?2, T. W. Bonner—496(L) 

Collimation of slow neutrons, J. G. Hoffman, M. S. 
Livingston—676(A) 

Distribution from Li+H?, M. L. Pool—707 

Emission from A, Cl, Al and under He* bombardment, 
E. Pollard, H. L. Schultz, G. Brubaker—351 

Interaction with normal and parahydrogen, J. R. Dun- 
ning, H. J. Hoge, J. H. Manley, F. G. Brickwedde 
205(A) 

From Li+H?, W. E. Stephens—212(A); 223 

Magnetic properties of free neutrons, O. Frisch, H. von 
Halban, Jr., J. Koch—719 

Magnetic scattering and magnetization of iron, H. G. 
Beyer, H. Carroll, C. Witcher, J. R. Dunning—947(A) 

Minimum energy for disintegration, R. Sagane—492(L) 

Multiple ejection from excited nuclei, R. L. Thornton, 
J. M. Cork—922(A) 

Passage through lead, E. Hudspeth, T. W. Bonner- 
928(A) 

Properties of slow neutrons, J. G. Hoffman, M. S. 
Livingston—929(A) 

Scattering of fast neutrons, R. F. Bacher, D. C. Swanson 
—676(A); 922(A) 

Scattering by gases, H. Carroll, H. G. Beyer, K. Wilhelm 
J. R. Dunning—947(A) 

Spin, H. H. Goldsmith, L. Motz—947(A) 


Nuclear moments and spin (see also Hyperfine structure) 


Electron pair theory of particle interaction, R. E. Mar- 
shak, H. A. Bethe—677(A); Erratum 950(A) 
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Nuclear moments and spin (continued) 
Hyperfine structure of spectral lines from atomic beams, 
K. W. Meissner—-677 (A) 
Of In™5, S. Millman, I. I. Rabi, J. R. Zacharias—331(A); 
384 
Of Li and K", J. E. Gorham—331(A) 
Of Li’; magnetic, H. A. Bethe—842(L) 
Magnetic moment of Li’, perturbation theory, D. R. 
Inglis—880; 939(A) 
Magnetic moments of Li®, Li? and F9, I. I. 
Millman, P. Kusch, J. R. Zacharias—495(L) 
Magnetic moments of Li® and K", signs of, J. E. Gorham 
—563 
Magnetic moments; theory, D. R. Inglis—470 
Of Mn, R. A. Fisher, E. R. Peck—685(A) 
Measuring nuclear magnetic moment, I. I. Rabi, J. R. 
Zacharias, S. Millman, P. Kusch—318(L) 
Of Ne and A, J. M. B. Kellogg, N. F. Ramsey—331(A) 
Quadrupole and magnetic moments, K. Way—685(A) 
Relativistic spin-orbit coupling, S. M. Dancoff, P. 
Morrison—211(A) 
Spin of neutron, H. H. Goldsmith, L. Motz—947(A) 
Nucleus (see also Disintegration of nucleus; Energy states 
of nucleus) 
Barytron theory of nuclear forces, H. A. Bethe—938(A) 
Dynaton theory of nuclear forces, R. Serber—211(A) 
Electromagnetic properties of nuclear systems, W. E. 
Lamb, Jr., L. I. Schiff—651 
Excitation of nuclei, V. Weisskopf—1018(L) 
Ground state of deuteron, R. G. Sachs, M. Goeppert- 
Mayer—991 
Penetrability of potential barrier, F. C. Hoyt—673(L) 
Relativistic equations, nuclear particles, G. Breit—153 
Rotation of atomic nuclei, E. Teller, J. A. Wheeler— 
684(A); 778 
Saturation of nuclear forces, C. Critchfield, E. Teller— 
684(A); 812 
Saturation requirements for forces, G. Breit, E. Wigner 
—998 


Rabi, S. 


Optical constants and properties 
Of fluorite in the 6u region, D. Williams—930(A) 
Intermetallic compounds, A. P. Friesen—688(A) 
Of liquid thallium, A. P. Friesen—678(A) 
Lorenz-Lorentz constant for paraffin oil, A. P. Friesen— 
688(A) 
Refractive index and dispersion of distilled water, L. W. 
Tilton, J. K. Taylor—918(A) 
Refractive index of air, D. Bender—852(A) 
Optical instruments (see also Methods and instruments) 
Oblique ether-drift interferometer, W. B. Cartmel— 
949(A) 
Optical theory 
Ether-drift interferometer, theory of, W. B. Cartmel— 
949(A) 
Nonideal Lloyd’s mirror, F. W. Doermann—420 


Packing fraction (see also Isotopes) 
Energy of heavy nuclei, A. J. Dempster—64; 869 
Of Ti, Fe, Cu, Ag, Au and Pt, A. J. Dempster—202(A) 


SUBJECT 


INDEX 


Photoelectric effect and properties; cells 
Energy distribution of photoelectrons from Na, A. G. 
Hill—184 
Sensitization of Al, L. P. Thein—287 
Spectral sensitivity of Sr and Mg, R. J. Cashman, E. 
Bassoe—919(A) 
Thoriated tungsten surfaces, A. King—570 
Threshold and intensity from Bi film, A. H. Weber—895 
Photography 
Camera error and photogrammetric mapping, I. C. 
Gardner—945 (A) 
Variation of contrast with wave-length, W. R. Koch— 
325(A) 
Photovoltaic effect 
In polarized light, R. G. Wilson—264(L) 
Piezoelectric effect 
Unusual demonstrations with piezoelectric resonators, 
K. S. Van Dyke—686(A) 
Vibration modes of low decrement, K. S. Van Dyke— 
945(A) 
Positive electron 
Energy loss in Al, B. R. Curtis—986 


Quantum mechanics (see Mechanics, quantum) 


Radiation (see also Thermal radiation) 

Asymmetric radiation produced by 6-rays, G. B. Collins, 
V. G. Reiling—205(A); 948(A) 

Discharges under nearly complete ionization, F. L. 
Mohler—683(A); 923(A) 

XT relation for black-body radiation, F. Benford, A. G. 
Worthing—948(A) 

Luminosity of yellow and green Hg lines, W. E. Forsythe, 
B. T. Barnes, A. L. Shrider-—678(A) 

Radioactivity (see also Disintegration of nucleus) 

Of air, water, cave and soil gases, O. G. Fryer—325(A) 

Alpha-particle mean range standards, M. S. Livingston, 
M. G. Holloway—922(A) 

Of As, B. R. Curtis, J. M. Cork—681(A) 

8-emission, nuclear K-electron capture, E. C. Campbell 
—206(A) 

8-spectrum of RaE, A. I. Alichanian, A. I. Alichanow, 
B. S. Dzelepow—766(L) 

B-spectrum of ThC, A. I. Alichanian, S. J. Nikitin— 
767(L) 

Of Be’, R. B. Roberts, N. P. Heydenburg, G. L. Locher— 
1016(L) 

Of Co, J. J. Livingood, G. T. Seaborg—847(L) 

In Cu, long period, S. W. Barnes, G. Valley—946(A) 

Of Cu, Zn, Ga and Ge isotopes, R. Sagane—212(A) 

Deposition of Ra Cu on rotating cathode, J. Steigman— 
771(L) 

Disintegration constant of thorium, A. F. Kovarik, N. I. 
Adams, Jr.—928(A) 

Energy of 8-rays, G. E. M. Jauncey—106(L) 

y-radiation from B+H2?, W. A. Fowler, E. R. Gaerttner, 
C. C. Lauritsen—942(A) 

y-radiation from In, B. R. Curtis, J. R. Richardson— 
942(A) 
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y-radiations from Eu’? and Au’, J. R. Richardson— 
942(A) 

Gamma-radiation from N!3, J. R. Richardson—610(L) 

Gamma-rays from I, A. Roberts, J. W. Irvine, Jr.— 
609(L) 

Gamma-rays from V*, N'3, Cu, Mg?? and Na®, J. R. 
Richardson—124 

Of In, by slow neutrons, A. C. G. Mitchell—269 

Induced by high energy protons, L. A. DuBridge, S. W. 
Barnes, E. O. Wiig, J. H. Buck, C. V. Strain—326(A); 
C. V. Strain, J. H. Buck—943(A); L. A. DuBridge, 
S. W. Barnes, J. H. Buck, C. V. Strain—447 

Induced in Ni, Cu, Mg and Ag by H!, L. N. Ridenour, 
L. A. Delsasso, M. G. White, R. Sherr—770(L) 

Induced in the rare earths by neutrons, M. L. Pool, 
L. L. Quill—326(A); 437 

Of I, J. J. Livingood, G. T. Seaborg—1015(L) 

Isomerism in Ra Ag and K electron capture, M. L. Pool, 
E. C. Campbell—272 

Isomers of In™*, A. C. G. Mitchell, L. M. Langer—505; 
675(A); M. L. Pool—611(L) 

Half-life of actino-uranium, A. O. Nier—922(A) 

Mean life of Ru 87, E. J. Konopinski, H. A. Bethe— 
679(A) 

Of Ni®? and Ni®, J. J. Livingood, G. T. Seaborg—847(L) 

Of K, A. Bramley, A. K. Brewer—502 

K electron capture in gallium, L. W. Alvarez—606(L) 

Produced by cosmic-rays, W. E. Ramsey, C. G. Mont- 
gomery, D. D. Montgomery—196(L) 

Proton bombardment of various elements, R. Sherr, 
M. C. Henderson, M. G. White, L. A. Delsasso, L. N. 
Ridenour—946(A) 

Range and specific ionization of a-particles, M. G. 
Holloway, M. S. Livingston—676(A) 

Of Sc, J. M. Cork, R..L. Thornton—866 

Short-lived radioelements, L. W. Alvarez—215(A) 

In Ag, M. L. Pool—116; 206(A) 

Substandards of feeble radioactivity, R. D. Evans, C. 
Goodman—916(A) 

Of Ta, O. Oldenberg—35; 109(A) 

Of Te, G. F. Tape, J. M. Cork—676(A) 

In Ur, P. Abelson—211(A) 

Of Zn, R. L. Thornton—326(A) 

Of Zn by a-particle bombardment, W. B. Mann—212(A) 

Of Zn and Co, C. Perrier, M. Santangelo, E. Segré— 
104(L) 

Raman spectra 

Of disubstituted acetylenes, F. F. Cleveland, M. J. 
Murray—330(A); 931(A) 

Of S, S. Bhagavantam—1015(L) 

Reflecting power (see also Absorption and reflection of 
light) 

Infra-red spectrometer with fluorite prisms, E. S. Barr— 
852(A) 

Relativity 

Eddington’s treatment of Dirac’s equation, B. Podolsky 
—591 

New ether-drift interferometer, W. B. Cartmel—108(A) 

Oblique ether-drift interferometer, W. B. Cartmel— 
949(A) 


Relativistic equations; nuclear particles, G. Breit—153 
Trouton-Nobel experiment, P. B. Carwile—924(A) 
Resistance, electrical (see Electrical conductivity and 
resistance) 


Scattering of atoms and molecules (see also Atomic and 
molecular beams) 
a-particles by Ne and H?, J. T. McCarthy—30 
Scattering of H! by H!, L. R. Hafstad, N. P. Heyden- 
burg, M. A. Tuve—239 
Scattering of electrons, neutrons, and ions (see also Elec- 
trons, scattering of; Electron diffraction) 
Angular distribution, J. G. Hoffman, M. S. Livingston 
1020(L) 
Atomic electron velocities, A. L. Hughes, M. M. Mann, 
Jr.—S0 
Effect of valence electrons on intensities, V. A. Johnson, 
H. M. James—327(A) 
Energy loss of positrons, B. R. Curtis, J. M. Cork 
206(A) 
Of fast electrons in gases, J. B. H. Kuper—993 
Of fast electrons and positrons in Pb, J. Oppenheimer, 
W. A. Fowler—928(A) 
Of fast neutrons, elastic and inelastic, D. C. Grahame, 
G. T. Seaborg—795 
Of neutrons by deuterons, P. G. Kruger, W. E. Shoupp, 
R. E. Watson, F. W. Stallmann—1014(L) 
Of He* by A, O, and Ne, G. Brubaker—921(A) 
Inelastic scattering by Ag, J. C. Turnbull, H. E. Farns- 
worth—944(A) 
Magnetic scattering and magnetization of iron, H. G. 
Beyer, H. Carroll, C. Witcher, J. R. Dunning—947(A) 
Multiple scattering of electrons, M. E. Rose—937(A) 
Multiple scattering of neutrons, O. Halpern, R. Luene- 
burg, O. Clark—173 
Of neutrons by gases, H. Carroll, H. G. Beyer, K. 
Wilhelm, J. R. Dunning—947(A) 
Of neutrons by gases, H. Carroll, P. N. Powers, H. G. 
Beyer, J. R. Dunning—680(A) 
Of neutrons by paramagnetic media, J. H. Van Vleck— 
680(A) 
Of neutrons by protons, theory of, R. G. Sachs, M. 
Goeppert-Mayer—991 
Of neutrons by various elements, W. H. Zinn, S. Seely, 
V. W. Cohen—921(A) 
Of H! and H? by deuterium and by helium, N. P. Hey- 
denburg, R. B. Roberts—922(A) 
Profile of Compton band determined by electron scatter- 
ing, A. L. Hughes, M. M. Mann, Jr.—204(A) 
Resonance scattering of charged particles, O. Halpern, 
M. H. Johnson, Jr.—939(A) 
Scattering of fast neutrons, R. F. Bacher, D. C. Swanson 
—676(A); 922(A) 
Of slow positive ions, A. G. Emslie—680(A) 
Scattering of light (see also Raman spectra) 
Diffraction of electromagnetic waves, E. Guth—688(A) 
Secondary electrons (see Electrons, secondary) 
Spectra, absorption (see also Absorption of light) 
Atmospheric absorption by water, W. M. Elsasser— 
768(L) 
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Spectra, absorption (continued) 

Of benzene at 2700-2200A, H. Sponer, G. Nordheim 
A. Sklar—932(A) 

Of C1302! at 4.3754, A. H. Nielsen—983 

Of CsH, G. M. Almy, M. Rassweiler—890 

Of diatomic fluorides in carbon-tube furnace, F. A. 
Jenkins, G. D. Rochester—213(A) 

4.3u band in CO:, D. M. Cameron, H. H. Nielsen—246 
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